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ABSTRACT 


This  effort  was  expended  to  investigate  the  feasibility  of  using  an  optical 
scanning  system  to  inspect  printed  circuit  boards  for  defects.  To  imple¬ 
ment  this  effort  a  Data  General  C/300  Eclipse  was  used  interfacing  with  a 
Colorada  Video  digitizer  and  a  Panasonic  TV  Camera /Monitor  system.  It 
was  found  that  the  system  could  operate  well  on  the  bare  P.  C.  board  without 
"pre-schooling and  would  reveal  imperfections  in  circuit  lands,  holes,  and 
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SUMMARY 


The  expended  effort  was  intended  to  determine  how  optical  techniques 
might  be  applied  to  an  automated  fault  detection  board  scanning  system. 
The  specific  faults  shown  below  were  analyzed  and  a  summary  of  the 
results  are  depicted  in  the  Table  below. 


1.  Base  electrical 

a.  Delamination 

b.  Measling 

2.  Component  Alignment 

3.  Conductor  Flaw 

4.  Lifted  Pad 

5.  Scratches 

6.  Surface  Quality 

7.  Hole  Registration 

8.  Legibility 

These  analyses  will  be  discussed  in  detail  in  the  following  section.  In 
summation,  all  areas  demonstrated  a  positive  result  with  the  one  exception 
of  the  component  alignment  check.  The  approach  selected  eminated  a  flair 
for  two-dimensional  examination  but  had  difficulty  with  any  "3-D"  require¬ 
ments.  In  order  to  improve  upon  this,  a  higher  resolution  system  (per 
unit  time)  would  have  to  be  used,  such  as  a  solid-state  video  psuedo  stare- 
state  mode  monitor  that  would  register  a  vertical  line  of  images  at  a  time, 
and  choose  these  data  points  that  are  significant.  This  would  involve  a 
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hardware  pre-processing  to  prevent  excessive  CPU  time.  Direct  benefits 
are  obvious  on  a  production  scale  with  the  currently  developed  software  and 
a  solid  state  scanner. 

"Pre-schooling"  would  enhance  the  inspection  process,  however,  the  major 
thrust  of  the  effort  was  directed  toward  a  general  workmanship  specifica¬ 
tion  approach.  This  was  achieved  with  the  one  exception  of  component 
alignment  checks.  "Pre-schooling"  has  been  shown  to  greatly  improve  this 
form  of  analysis. 


II. 


INTRODUCTION 


This  document  constitutes  the  final  submittal  for  Contract  DAAK40-77-C- 
0105. 

Negotiations  and  final  signing  of  the  contract  were  completed  on  May  3,  1977. 
Since  cost  data  is  available  at  SCI  using  months  of  4,  4,  and  5  weeks  per 
month  during  each  quarter,  the  monthly  progress  and  monthly  financial 
reports  were  based  on  this  same  schedule. 

Per  the  agreements  prior  to  negotiation  of  the  contract,  SCI  provided  a  com¬ 
puter  system  for  use  while  this  work  was  being  conducted  at  SCI.  A  Data 
General  C/300  Eclipse  computer  system  was  moved  into  the  work  area  and 
checked  for  proper  operation.  Additional  electrical  service  was  provided  to 
the  lab  in  order  to  accommodate  the  computer  system  plus  other  normal 
laboratory  equipment.  The  elements  of  the  system  are  depicted  in  Appendix 
A. 

An  evaluation  of  available  video  digitizers  resulted  in  a  decision  to  purchase 
an  image  digitizer,  Model  270A  by  Colorado  Video,  which  was  easier  to 
interface  and  operate  than  others. 

Other  equipment  in  the  system  included  a  graphics  monitor /terminal,  a  TV 
camera,  and  a  video  monitor.  A  paper-tape  reader  was  also  interfaced  to 
the  C300  computer.  It  was  used  to  allow  access  of  a  wider  variety  of  general 
support  software  than  was  available  on  magnetic  tape. 

The  task  of  locating  sources  for  board  defects  was  conducted  along  with  a 
review  of  board  and  solder  specifications. 


III. 


DISCUSSION 


Implementation  of  Hardware 

A  standard  unmapped  real  time  disk  operating  system  (RDOS  Rev.  5)  was 
installed  in  the  computer  and  various  service  routines  checked  out.  This 
was  later  replaced  by  INFOS.  An  illumination  system  was  designed,  which 
included  a  PC  board  mounting  technique. 

An  illumination  box  was  designed  to  allow  proper  illumination  of  the  printed 
circuit  board  under  both  diffuse  and  spot  lighting  conditions  with  a  manually 
moveable  mount  for  the  board  located  on  the  illumination  box.  A  printed 
wiring  board  was  mounted  on  the  outside  of  the  box  so  that  the  camera,  which 
is  mounted  on  the  opposite  side,  could  see  the  board  through  holes  in  the  box. 
Lamps  were  also  mounted  on  the  box  such  that  either  a  spot  can  be  placed  on 
the  board  or  scattered  light  from  inside  the  box  impinges  on  the  board.  The 
inside  of  the  box  is  both  diffuse  and  white  to  get  maximum  scattering  of  the 
light.  A  television  camera  was  selected  to  serve  as  the  primary  sensor. 

With  the  camera  on,  tests  were  made  on  the  preliminary  software.  These 
tests  indicated  that  more  resolution  would  be  required  than  was  available 
using  the  standard  lens  which  came  with  the  camera.  An  adaptor  was  made 
to  allow  the  use  of  a  photographic  camera  lens  and  extension  tube.  This 
increased  the  resolution  significantly. 

The  camera  and  lens  and  two  lamp  sources  were  mounted  to  a  wooden  base, 
and  the  diffusion  box  was  fastened  to  this  base.  With  this  system,  the  diffuse 
illumination  was  satisfactory.  The  specular  illumination  scheme  did  not 
produce  the  anticipated  results.  This  was  due  to  lack  of  reflectance  of  the 
board  material  and  to  a  dullness  of  some  of  the  solder  on  the  sample  board. 
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This  condition  is  normally  encountered  in  industry  standards  and  while  it  is 
a  negative  result,  it  should  'ce  noted  as  an  approach  not  to  be  attempted. 

For  convenience  and  clarity,  the  technical  requirements  that  guided  this 
effort  are  included  in  the  following  section. 
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TECHNICAL  REQUIREMENT  NO.  6148 


1.  SCOPE 

1.1  GENERAL 

The  purpose  of  this  project  was  to  determine  methods  and  technology  of 
inspecting  printed  wiring  boards  at  the  bare  board  or  populated /soldered 
level  by  an  automatic  video  scanning  method  and  to  demonstrate  the  feasibility 
of  such  a  process  as  a  cost  reducing  method  for  the  determination  of  an 
optimum  methodology  for  the  location/evaluation  of  specific  classes  of  ' 
defects. 


2.  APPLICABLE  DOCUMENTS 


2.  1 

MlL-STD-454  -  Standard  General  Requirements  for  Electronic  Equipment. 

2.2 

MlL-STD-275  -  Printed  Wiring  for  Electronic  Equipment. 

2.3 

MIL-STD-1495  -  Multilayer  Printed  Wiring  Boards  for  Electronic  Equipment. 


2.4 

MIL-S-45743  -  Soldering,  Manual  Type,  High  Reliability. 
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2.  5 

MIL-S-46344  -  Solder  Bath  Soldering  of  Printed  Wiring  Assemblies,  Automatic 
Machine  Type. 

3.  REQUIREMENTS 

3.  1  TECHNICAL  EFFORT 

The  contractor  provided  the  personnel,  materials,  equipment  and  expertise 
to  perform  the  tasks  as  hereinafter  described. 

3.2  PLANNING 

Technical  effort  conducted  under  these  technical  requi rements  were  planned 
in  terms  of  technical  milestones  to  be  accomplished  in  pursuit  of  the  program. 
Each  milestone  indicated  completion  of  a  significant  portion  of  the  effort.  A 
milestone  chart  showing  the  entire  program  was  prepared  by  the  contractor 
and  submitted  in  an  oral  presentation  within  thirty  (30)  days  after  the  contract 
award.  The  government  reserved  seven  (7,  days  after  receipt  to  review  and 
concur/nonconcur  in  the  detailed  proposed  plan;  if  the  contractor  did  not 
receive  notification  within  seven  (7)  days,  he  was  to  assume  concurrence. 

The  milestone  chart  and  the  government  concurrences/nonconcurrences 
were  documented. 

3.3  DETAILS  OF  EFFORT  TO  BE  PERFORMED 

3.  3.  1 

The  contractor  demonstrated  an  image  processing  system  consisting  of  a 
contractor  supplied  computer  and  the  following  items  obtained  and  assembled 
as  a  part  of  this  contract: 
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iHt  lMkUEST  fHfrt  Yn.lt.  IS  4"  iHfc  SullNu  fncTu* 


F  1  LENrti-ib  if. 


,'TuUP  THRESHOLDS  BET  WEEN  42  AND  101 

HOI  t  iN  HFEH  1 

'MINIMUM  TH I UK NESS  HRuUND  HOLE  iS  .  12  21 

HOLE  L  EN  l  ER—  27..  7 

AVERAGE  RHD  1 1  If  =  5.  41 

RADIUS  RANGE  FROM-  4  72  1 U  6.  5b 
CIRCUMFERENCE  OF  HOLE  If  15.  91 
AREA  OF  HOLE  IS  117.  50 
4+ PI +HREA/ < C ♦ *2 > =  1.  14 

LOCH  7  x  ON  -  IMN*  22  I  MX*  32  ,.MN-  72  JMX-  82 

HR  EH  1  BRIGHTNESS  BETWEEN  0  AND  10b 

PER i METER*  249.  54  INTERNAL  AREA*  1790  00 

MINIMUM  INTERNHL  THICKNESS-  8.  75 
I M I N—  10  I MAX*  43  JMIN*  1  JHRX*  97 
RREH  CON TH I NED  1  HOLES  RND  0  SPOTS 

HOLE  IN  RREH  2 

MINIMUM  THICKNESS  HRUUND  HOLE  IS  12 
HOLE  CENTER—  SO*  31 
AVERAGE  RADIUS*  5.  75 
RHDIUS  RANGE  FROM-  4.  25  TO  7.  76 
CIRCUMFERENCE  OF  HOLE  IS  37.  91 
RREH  OF  HOLE  IS  136  25 
4  *F  I  +  HREH.'  ■■.  C++ 2  >  —  1.  19 

LOCHTION  -  I MN=  56  I  MX*  66  JMN=  26 

RREH  :  2  BRIGHTNESS  BETWEEN  0  RND  106 

PERIMETER*  152.  08  INTERNAL  AREA-  1329.  28 

MINIMUM  INTERNAL  THICKNESS*  10.  OO 
IMIN*  45  IMHX*  77  JMIN*  1  J  MAX*  51 

RREH  CONTAINED  1  HOLES  RND  0  SPOTS. 

AREA  THICKNESS  INDE1 ERMINHTE 

AREA  3  BRIGHTNESS  BETWEEN  0  AND  106 

PERIMETER*  72.  86  INTERNAL  AREA*  186.  88 

MINIMUM  INTERNAL  THICKNESS*  *20000 
IMIN*  1  I MAX*  7  JMIN*  16  JMAX*  49 

AREA  CONTAINED  O  HOLES  AND  O  SPOTS. 

THERE  WERE  3  AREAS  LOCATED. 
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FIGURE  4  PRINT  OF  CHARACTERISTICS  OF  BOARD  IN  FIGURE  2 
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FILENAME:  IS  XXX  t 

CONTOUR  THRESHOLDS  BETWEEN  27  AND  100 
H<Jl_E  IN  HREA  1 

MINIMUM  THICKNESS  HKUUND  HOLE  Is  .  9.  62 

HOLE  LENT  EE!—  2L-  ■  71 

AVER'HGE  RAD  iiJS—  0.  _-U 

RADIUS  E'HNijE  FROM—  S.  SI.  t  U  6. 

c i eilumfee.enl e  uf  hole  is.  ..t.^ui 

heeh  of  hole  is  iol.  >5 

4  t-F  1  ♦  HR  E  AX  •'  i_!  ♦  +  2  >  —  1  1 

1.  ULHT  ION  -  iMN-  IS  I  MX-  2a  JMN*  67  JHX-  76 

hE'EH  .1.  BE:  I  OH  I  NESS  BETWEEN  O  HND  lo6 

E'ER  i  ME  TER-  244.  IS  INTERNAL  HR  Eft-  1546.68 

MINIMUM  INTERNAL  THICKNESS*  8.  7b 
1M1N*  11  I MAX*  40  JMIN-  1  JMAX*  S4 

HREH  I.  ON  I  H 1  NED  1  HOLES  HND  0  SPOTS. 


HOLE  1 N  HREH  2 

MINIMUM  THICKNESS  HROUNO  HOLE  is  10.  00 

HOLE  CEN  I  ER—  58.*  25 

AVERAGE  RADIUS*  5  25 

RHDIUS  RANGE  FROM—  2.  Si  T  O  €•.  S_- 

CIRCUMFERENCE  OF  HOLE  IS  24.  61 

HREH  OF  HOLE  IS  102.  ?5 

4+P I  +-AKEAX  C++- 2 )  —  1.  12 

LOCATION  -  IMN*  54  I MX*  oi  JMN*  il  OMX*  20 

AREA  :  2  BRIGHTNESS  BETWEEN  0  AND  106 

PERIMETER*  146.  72  INTERNAL  HREH*  1218.  75 

MINIMUM  INTERNAL  THICKNESS*  10.  U0 
I M I N*  46  I MAX*  78  JMIN*  1  JMAX*  48 
AREA  CONTft I NED  1  HOLES  HND  0  SPOTS. 

AREA  THICKNESS  INDETERMINATE 


AREA  :  2  BRIGHTNESS  BETWEEN  0  AND  106 

PERIMETER*  79.  06  INTERNAL  AREA*  225.  62 
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FIGURE  3  PRINT  OF  CHARACTERISTICS  OF  BOARD  IN  FIGURE  1 


A  specific  sequence  of  operation  will  appear  as  follows: 


1.  The  DOTM AT  Routine  is  called  and  responds  with  the  question 
"File  Name?"  which  directs  the  entry  of  the  desired  data  file. 
This  data  file  is  acquired  by  use  of  live  scan  using  routine  SCAN 

2.  The  name  of  the  file  is  arbitrary  but  must  not  exceed  10 
characters  in  length. 

2.  The  DOTM  AT  Routine  provides  the  information  to  define  the 
total  number  of  detected  points.  This  information  can  be  used 
to  compare  with  the  expected  norm,  or  to  establish  a  norm. 

Also  provided  by  DOTMAT  is  visual  image  of  the  area  being 
scanned.  PRTM AT  is  a  companion  routine  that  outputs  on  the 
line  printer  when  a  hard  copy  is  required. 

3.  HISTOGRAM  is  run  on  the  data  file  created  by  SCAN  2  and  dis¬ 
plays  the  information  needed  to  select  upper  and  lower  thresholds 
to  be  used  when/if  PCCARD  is  used.  If  Step  1  is  used  on  a  test 
card  and  the  printout  "Probable  Error  on  Board"  appears,  then 
Steps  2  and  3  apply.  Otherwise,  the  test  series  would  complete 
with  Step  1.  The  CRT  output  of  HISTOGRAM  is  shown  in  Figure  7. 

4.  PCCARD  is  activated  by  a  series  of  inputs: 

a.  "Print  or  Type"  Enter  P  or  T 

b.  "Live  Scan?"  Enter  Y  or  N 

c.  "Enter  Filename" 

* 

d.  File  is  opened  and  question  displayed  "Smooth  Data?" 

Y  or  N 

e.  Thresholds  are  defined  and  are  optional;  override  entry 
is  provided. 

The  output  information  provided  includes  the  minimum  distance  from  exterior 
and  interior  surfaces  of  the  solder  lands,  and  the  data  displayed  by  Figure  3. 
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System  Oce  ration 


Operating  the  system  under  INFOS  the  operator  first  calls  for  the  PCBRD 
directory.  Then  he  has  a  choice  of  several  analytical  routines  to  use  de¬ 
pending  on  the  nature  of  the  inspection.  First  a  "standard"  is  recorded  of  a 
normal  board  and  stored.  The  test  boards  are  scanned  and  compared.  The 
recommended  procedure  is  to  use  the  HISTOGRAM  Routine  to  determine  the 
upper  and  lower  threshold  limits,  and  the  DOTMAT  Routine  to  establish  the 
number  of  points  detected.  This  can  then  be  followed  by  PCCARP  to  analyze 
the  nature  of  the  card  in  question  in  regard  to  a  given  spot.  This  above 
procedure  is  presented  in  more  detail  next. 

A  printed  circuit  board  should  be  tested  in  the  following  manner.  A  sample 
card  is  pulled  from  a  production  run.  The  board  is  scanned  with  TV  monitor 
system  and  the  digitized  data  is  compared  with  the  standard.  These  two 
scans  can  be  given  arbitrary  file  names.  In  this  example,  the  scans  are 
given  the  names  XXX2  for  the  standard,  and  XXXI  for  the  sample  board. 

If  more  than  a  1%  difference  in  the  number  of  detected  points  occurs,  an 
error  message  will  appear  as  "probable  Error  on  Board".  If  not  "Board 
OK"  will  appear.  Following  this  if  further  checks  are  needed,  then  the 
HISTOGRAM  or  PRTHIST  routines  are  employed  to  determine  the  threshold 
settings  to  use.  Next,  the  PCCARP  Routine  is  implemented  and  will  call  in 
the  data  file  and  will  determine  the  minimum  distance  across  any  PC  land 
and  will  display  it  on  the  Tektronix  memory  terminal.  (See  Figures  1  &  2). 
Also,  Figures  3  &  4  depict  information  that  can  be  observed  in  hard  copy  or 
CRT  display.  In  this  copy,  the  error  shown  is  an  off-center  hole  in  a  PC  land. 
This  is  visible  in  Figure  1.  The  units  displayed  are  mils.  PRTHIST  depicts 
a  plot  of  the  relative  rate  of  occur  ranee  of  light  levels  on  an  arbitrary  scale 
of  0  to  256.  The  routine  HISTOG RAM  displays  the  same  plot  on  the  CRT 
terminal,  but  does  not  display  the  threshold  card  scale  factor  data.  The 
output  of  PRTHIS  is  displayed  in  Figures  5  &  6  for  the  boards  in  Figures  1  &  2 


This  final  report  in  narrative  form  was  delivered.  This  report  summarizes 
the  work  accomplished  under  this  contract.  This  report  fully  documents 
the  equipment,  its  interfacing,  and  programming  manuals  required  to  duplicate 
and  utilize  the  demonstrator  scanner  for  inspecting  printed  wiring  boards 
within  a  production  environment. 


3.3.3 

Printed  circuit  boards  with  dimples  of  several  levels  of  the  defects  under 
consideration  were  obtained.  A  rating  system  was  devised  to  establish  the 
severity  of  each  defect  and  to  allow  a  decision  as  to  acceptable /non-acceptable 
for  each  type  defect. 

3.  3.4 

Methods  of  detecting  and  evaluating  printed-ci rcuit  board  defects  were 
devised  and  tested  using  the  system  of  3.  3.  1.  Implementation  was  in 
software  wherever  possible. 

3.  3.  5 

Based  on  the  above,  this  report  was  prepared  containing: 

(A)  Discussion  of  the  defects  and  related  scoring  system  developed 
in  3.  3.  3. 

(B)  Discussion  and  explanation  of  the  detection  methods  developed  in 
3.3.4. 

4.  DOCUMENTATION 

4.  1 

Data  delivered,  schedule  of  delivery,  and  distribution  requirements  are 
specified  on  DD  Form  1423. 

4.2 

In  addition,  the  contractor  reviewed  the  progress  of  the  work  to  the  government 
at  mutually  agreed  times  and  places.  The  number  of  such  reviews  held  at 
places  other  than  the  contractor's  facilities  did  not  exceed  three  (3). 
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(A)  Input  Imaging  Device  with  computer  interface  containing  vidicon 
TV  camera,  scan  control,  analog /digital  conve  rte r,  compute  r  bus 
interface  electronics. 

(B)  Graphics  display  terminal  capable  of  accepting  images  from  the 
computer,  storing  the  image,  and  displaying  the  image  with 
several  shades  of  gray. 

(C)  Mechanical  system  capable  of  holding  and  manipulating  printed 
circuit  boards  under  manual  control  for  placement  under  the 
imaging  and  illumination  system. 

(D)  Image  processing  software  capable  of  controlling  and  operating 

the  Input  Imaging  Device  in  order  to  bring  images  into  the  computer 
memory  along  with  special  processing  of  these  images  to  include 
thresholding  and  variation  of  number  of  resolution  elements. 

3.3.2 

Board  Parameters  evaluated  include: 

(A)  Rough  or  ragged  edges 

(B)  Rough,  burred  holes 

(C)  Haloing  extending  between  conductors 

(D)  Measling 

(E)  Delamination 

(F)  Localized  blistering 

(G)  Exposed  glass  weave 

(H)  Board  crazing 

(I)  Legibility  of  screened,  printed  or  etched  marking 

(J)  Plating  modulation  in  holes 

(K)  Holes  with  voids 

(L)  Bridged  circuit  lines 


FIGURE  7  CRT  DISPLAY  BY  HISTOGRAM 
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FIGURE  8  EXAMPLE  OF  MEASUNG 


The  type  of  system  faults  that  are  common  on  printed  circuit  cards  are  as 
follows: 

1.  Base  Electrical 

A.  Delemination 
5.  Measling 

2.  Component  Alignment 

3.  Conductor  Flaw 

4.  Lifted  Pad 

5.  Scratches 

6.  Surface  Quality 

7.  Hole  Registration 

8.  Legibility 

The  system  can  readily  spot  board  elamination  or  measling  but  has  difficulty 
in  differentiating  between  the  two.  Component  alignment  checks  have  to  be 
done  on  a  looser  tolerance.  Conductor  flaws  on  a  bare  board  are  the  most 
consistently  recognized  of  all.  An  example  of  measling  is  shown  in  Figure 
8  with  characteristics  print  on  Figure  9. 

Scratches  or  poor  surface  quality  can  be  detected  by  change  in  the  number  of 
detected  points  using  DOTMAT  or  PRTMAT.  The  use  of  PCCARD  can  then  be 
used  to  isolate  and  measure  the  problem  areas.  This  comparison  technique 
requires  that  the  light  be  maintained  constant  during  a  series  of  board  to 
minimize  the  number  of  false  alarms. 

Legibility  of  letters  and  numbers  can  be  examined  by  use  of  PCCARD  to 
determine  the  minimum  width  on  each  letter  or  number  and  displaying  the 
features  on  the  Tektronix  terminal.  A  scan  of  letters,  in  this  case  a  logo, 
is  shown  on  Figure  10.  ' 


r  I  ?  ENAMF  iS  XXX  J 

CONTOUR  THR  F'-  HOu/i  BE  i  WEEN  45  AM  •  62 

## ;LN  iUUi SMUMili. 

CONTOUR  VHRESHGi  DS  6£  rWEEN  45  AND  82 
HOLE  IN  ARr.n  1 

M I N I HUM  THICKNESS  AROUND  HOLE  IS  .24.  62 
HOLE  CENTER—  40,  LI 
AVERAGE  r  HD  .<  > .io~  2.  BO 

EHNi.it  F  R  ON—  1.  HO  Tm  0  EL- 
C I R OUMFEF ENCt  OF  HOLE  IS  .25.  16 
BFFH  OF  Hui.E  IS  57  50 
4*Rl  +  ARFA.4  L  13**2,'  =  1.  14 

L.QCATION  -  IMN*  So  I  MX*  44  .IMN  •  17  .JMX-  25 
HGi.F  I  N  fir  EH  1 

MINIMUM  THICKNESS  AROUND  HOLE  IS  : 

HGL  K  CENT  cp>  27,  3-4 
AVERAGE  RADIUS—  14  89 
RADIUS  RANGE  FROM—  2  00  TO  52. 

C IRCiJMEEr ENCB  OF  HOLE  IS  156  84 
AREA  OF  HOLE  I S  1358.  00 
4*  FI •+  ARE  AXv I,.+:;tL' >  —  O.  52 

LOCATION  -  I NN—  10  I  MX*  44  J'HN 

HOLE  IN  AREA  I. 

MINIMUM  THICKNESS  AROUND  HOLE  IS  5.  3:9 

HGL  t  CENTER-  1 7,  60 

AVERAGE  RADIUS-  7.  26 

RADIUS  RANGE  FROM—  2.  26  TG  12.  25 

C I RCUMFERENCE  OF  HOLE  IS  69.  42' 

AREA  OF  HOLE  IS  23:5  GO 
4  +  F I  * ARE AX v  C* +2)*  0  61 

LOCATION  -  IMN*  7  I MX-  26  JMN*  54  JMX*  68 

AREA  :l  BRIGHTNESS  BETWEEN  O  AND  104 

PERIMETER*  220  93:  INTERNAL  AREA-  4300.  62 

MINIMUM  INTERNAL  THICKNESS-  55  23. 

7  MIN-  l  1  MAX—  69  JMIN-  1  JilAX-  83 
AREA  CONTAINED  2  HOLES  AND  4  SHOTS. 

ONL V  ONE  AREA  SO  NO  DISTANCES  AVAILABLE. 


10.  61 


*  19  JMX-  52 


FIGURE  9  PCCARD  OUTPUT  FOR  MEASLING  EXAMPLE 
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Other  detectable  defects  are  shown  in  other  figures  included.  Figures  11 
and  12  illustrate  a  concentricity  defect.  Figures  13  and  14  illustrate  an 
example  of  nodule  defect.  Figures  15  and  16  depict  the  same  board  without 
a  nodule.  Figure  18  displays  the  results  from  PRTMAT  of  a  scratched 
board.  Figure  18  is  a  display  of  the  output  of  a  routine  named  CALGRID 
which  provides  a  method  of  calibration  of  absolute  distances.  Figures  19 
A  and  B  present  the  image  of  a  broken  land  and  Figures  20  A  and  B  present 
an  image  of  the  same  board  without  any  broken  land.  Figure  21  illustrates 
the  general  system  configuration. 
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X  .X  X  1 


FILENAME  IS 

contour  thresholds  between  37  and 

HOLE  IN  AREA  1 

MINIMUM  THICKNESS  mPOUND  HOLE  IS  .  5.  62 

HOLE  CENTER-  23,  71 

AVERAGE  RADIUS-  5.  3v 

RADIUS  RANGE  FROM-  3  91  TO  6.  93 

CIRCUMFERENCE  OF  HGi_E  IS  34  01 

AREA  OF  HOLE  IS  103  75 

4*PI*AREA/<C**2>*  1.  13 

LOCATION  -  IMN=  19  IMX*  28  UMN=  67  JMX= 

AREA  1  BRIGHTNESS  BETWEEN  0  AND  106 

PERIMETER*  244  19  INTERNAL  AREA*  1846  S3 

MINIMUM  INTERNAL  THICKNESS*  8.  75 
I MIN*  11  IMAX*  45  JMIN*  1  UMAX*  94 
AREA  CONTAINED  1  HOLES  AND  0  SPOTS. 

HOLE  IN  AREA  2 

MINIMUM  THICKNESS  AROUND  HOLE  IS  :  10  00 

HOLE  CENTER-  53,  25 

AVERAGE  RADIUS*  5.  25 

RADIUS  RANGE  FROM-  3.  91  TO  6.  93 

CIRCUMFERENCE  OF  HOLE  IS  34.  01 

AREA  OF  HOLE  IS  103.  75 

4*FT*AREA./(C**2>  =  1.  13 

LOCATION  -  IMN*  54  IMX*  63  JMN*  21  -JMX* 

AREA  :  2  BRIGHTNESS  BETWEEN  0  AND  106 

PERIMETER*  146.  73  INTERNAL  AREA*  1318.  75 

MINIMUM  INTERNAL  THICKNESS*  10.  00 
I MIN*  46  IMAX*  78  JMIN*  1  UMAX*  48 

AREA  CONTAINED  1  HOLES  AND  0  SPOTS. 

AREA  THICKNESS  INDETERMINATE 

AREA  .  3  BRIGHTNESS  BETWEEN  0  AND  106 

PERIMETER*  79.  06  INTERNAL  AREA*  225.  63 

MINIMUM  INTERNAL  THICKNESS*  *20000 
IMIN*  1  IMAX*  8  ,JMIN=  12  UMAX*  47 

AREA  CONTAINED  0  HOLES  AND  0  SPOTS. 

THERE  WERE  3  AREAS  LOCATED. 
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FIGURE  12A 


DATA  FOR  BOARD  IN  FIGURE  11 
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FIGURE  12  B  CONCENTRIC  HOLE 
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rgbaele  error  on  board 
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«««»  «»•««*•« 

»«««*  •«««••• 

«•#«« 

•tit###  *#«##• 

»*»##  #•####* 

•  •«»»  •  #•#•#•# 

It  It  It  tilt  II  #IMHHMHt«tt 

»««««  IMMtIHMttMtWW 

JMIItfttt  til  It  ##«#*••## 

•  •••  »IHMI«tMHMHHH>« 

WMffNM*  •IMMtIMMMMMMMMMt 

»«  It  *  *•*#  •*#«#••«»## 


•  •  •••••» 


#•## 

ttttt# 

•• 

tttt 
•••• 
•tttt 
It  It  ttltlt 

•  •tit 

•  ••»  It 

««•« 


FIGURE  13  EXAMPLE  OF  NODULE 
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FILENAME  IS  :  PCNCD 


CONTOUR  THRESHOLDS  BETWEEN  77  AND  120 

HOLE  IN  AREA  1 

MINIMUM  THICKNESS  AROUND  HOLE  IS  :  3.  00 

HOLE  CENTER-  75,  35 

AVERAGE  RADIUS*  14.  09 

RADIUS  RANGE  FROM-  11.94  TO  16.00 

CIRCUMFERENCE  OF  HOLE  IS  94.  23 

AREA  OF  HOLE  IS  753.  75 

4*PI*AREA/(C**2>=  1.  07 

LOCATION  -  IMN=  65  IMX=  87  JMN*  21  JMX*  51 

AREA  :  1  BRIGHTNESS  BETWEEN  0  AND  124 

PERIMETER*  268.  26  INTERNAL  AREA*  1368.  75 

MINIMUM  INTERNAL  THICKNESS*  9.  38 
IMIN*  44  IMAX*  92  JMIN*  13  UMAX*  100 
AREA  CONTAINED  1  HOLES  AND  0  SPOTS. 

HOLE  IN  AREA  2 

MINIMUM  THICKNESS  AROUND  HOLE  IS  :  3.  20 

HOLE  CENTER-  19,  37 

AVERAGE  RADIUS*  13.  82 

RADIUS  RANGE  FROM-  12.  00  TO  16.  01 

CIRCUMFERENCE  OF  HOLE  IS  92.  67 

AREA  OF  HOLE  IS  718.  75 

4*PI*AREA/<C**2)=  1.  05 

LOCATION  -  IMN*  8  IMX*  31  JMN*  24  JMX*  51 

AREA  :  2  BRIGHTNESS  BETWEEN  0  AND  124 

PERIMETER*  163.  64  INTERNAL  AREA*  809.  38 

MINIMUM  INTERNAL  THICKNESS*  9.  38 
IMIN*  1  IMAX*  35  JMIN*  14  UMAX*  64 
AREA  CONTAINED  1  HOLES  AND  0  SPOTS. 

THERE  WERE  2  AREAS  LOCATED. 

MINIMUM  DISTANCE  FROM  AREA  1  TO  2  IS  20.  16 
20.  16  47  56  33  46 

20.  16  47  56  33  46 


FIGURE  14  DATA  FOR  BOARD  IN  FIGURE  13 
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FIGURE  15  BOARD  WITHOUT  A  NODULE 
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FILENAME  IS  PCNONQD 

CONTOUR  THRESHOLDS  BETWEEN  68  AND  114 

AREA  .  1  BRIGHTNESS  BETWEEN  0  AND  118 

PERIMETERS  163.  34  INTERNAL  AREA=  531.  25 

MINIMUM  INTERNAL  THICKNESS*  5  00 
IMIN=  1  IMAX=  24  JM I N=  11  UMAX*  53 

AREA  CONTAINED  0  HOLES  AND  0  SPOTS. 

HOLE  IN  AREA  2 

MINIMUM  THICKNESS  AROUND  HuLE  IS  :  0.  00 

HOLE  CENTER-  65,  29 

AVERAGE  RADIUS*  15.  15 

RADIUS  RANGE  FROM-  12.85  TO  13.  11 

CIRCUMFERENCE  OF  HOLE  IS  104.  23 

AREA  OF  HOLE  IS  831.  88 

4*PI*AREA/(C**2>=  0.  96 

LOCATION  -  IMN*  54  I  MX*  78  JMN*  14  JMX*  45 

AREA  :  2  BRIGHTNESS  BETWEEN  0  AND  118 
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MINIMUM  INTERNAL  THICKNESS*  10.  00 
IMIN*  31  IMAX*  80  JMIN*  10  UMAX*  100 

AREA  CONTAINED  1  HOLES  AND  0  SPOTS. 

AREA  :  3  BRIGHTNESS  BETWEEN  0  AND  118 
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IMIN*  90  IMAX*  100  JMIN*  53  UMAX*  100 

AREA  CONTAINED  0  HOLES  AND  0  SPOTS 

THERE  WERE  3  AREAS  LOCATED. 
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FIGURE  16  DATA  FOR  BOARD  IN  FIGURE  15 
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F  ig  u  v  e  A  - 1  Data 


APPENDIX  A 

ELEMENTS  OF  THE  SYSTEM 


more  rapid  inspection  and  would  reduce  the  amount  of  memory- 
required.  The  hardware  would  not  access  thru  except  when  the 
desired  attribution  in  the  data  are  found. 

5.  A  relatively  small  computer,  possibly  a  microprocessor.  To 
produce  the  necessary  report  to  go  with  each  board.  With  the 
hardware  preprocessing,  a  relatively  small  amount  of  computa¬ 
tion  would  be  necessary. 
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holes  can  be  detected  accurately  enough.  The  problem  then  reduces  to  find¬ 
ing  a  device  with  a  large  number  of  resolution  elements,  and  a  fast  scan  rate. 
The  best  choice  for  this  is  probably  one  of  the  solid  state  line  arrays.  This 
type  of  array  could  be  combined  with  a  board  moving  mechanism  (like  a  con¬ 
veyor  belt)  to  give  a  rapid  high  resolution  scan.  The  software  would  need  to 
be  changed  to  accumulate  the  features  of  the  areas  in  real  time  so  that  storage 
of  the  complete  image  would  not  be  necessary.  The  storage  of  complete 
images  must  be  eliminated  if  high  resolution  scans  of  moderate  sized  boards 
are  to  be  accomplished.  In  order  to  get  a  greater  processing  speed,  it  may 
be  necessary  to  do  a  portion  of  the  computation  in  hardware. 

Additional  work  can  be  done  in  the  thresholding.  The  present  algorithms 
work  well  if  the  image  illumination  is  flat,  but  if  there  is  shading  of  the 
field,  then  the  segmentation  is  poorly  done.  Some  image  processing  systems 
have  eliminated  this  shading  problem  with  hardware  filters.  This  would 
certainly  be  a  plausable  approach. 

The  final  recommendations  for  inspection  of  drilled  and  etched  boards  would 
be  to  use: 

1.  A  solid  state  line  array  for  rapid  scanning. 

2.  A  continuously  moving  PC  board  ideal  for  production. 

3.  Hardware  Dynamic  Thresholding  -  In  order  to  do  better  scene 
segmentation,  a  thresholding  algorithm  based  on  rate  of  change 
of  contrast  will  probably  need  to  be  implemented  in  hardware. 

The  task  of  thresholding  and  scene  segmentation  is  crucial  to  all 
of  the  image  processing  and  should  receive  maximum  attention. 

4.  Hardware  Processing  -  Hardware  Processing  for  features  such 
as  area,  perimeter,  minimum  thickness,  etc.,  would  permit 
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While  the  current  system  is  a  very  excellent  analytical  system,  the  routines 
would  of  course  benefit  from  more  on-line  refinement  of  the  diagnostic  and 
detection  methods  employed.  This  would  provide  a  greater  attractiveness' for 
production  line  applications  of  the  system.  One  outstanding  ability  of  the 
system  is  searching  for  dimensional  variations  and  for  minimum  distances. 
Lack  of  concentricity  in  hole  locations  is  an  example,  vis-a-vis: 


Using  this  ability  an  inspection  could  be  made  of  a  sample  lot  and  the  result¬ 
ing  distribution  of  dimensional  variances  could  be  used  to  characterize  the 
quality  control  for  that  particular  production  run. 

The  next  effort  to  be  made  in  the  complete  board  automatic  inspection  should 
probably  be  in  the  training  portion.  That  is,  a  technique  for  storing  the  data 
on  a  "good"  board  in  some  compacted  form  should  be  developed.  One  suggested 
technique  is  to  store  the  printout  data  for  each  portion  of  the  board  scanned. 
This  is  considerably  more  compact  than  storing  image  files  and  has  the  advan¬ 
tage  that  it  is  less  sensitive  to  alignment.  The  particular  portion  of  the  print¬ 
out  data  that  would  be  useful  is  the  number  of  area's,  their  perimeters,  their 
areas,  and  the  number  of  holes.  The  distances  do  not  need  to  be  saved  since 
they  are  covered  by  specification. 

A  limiting  factor  on  the  speed  of  inspection  is  the  scanning  device. 

The  scanning  resolution  of  approximately  0.004"  needs  to  be  maintained  so 
that  the  widths  can  be  resolved  and  the  hole  roundness  on  0.  020  inch  di  ^eter 
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The  fact  that  the  computer  is  tireless  permits  100%  inspection  of  boards  to 
the  spacing  and  line  width  specifications.  It  is  also  simple  to  100%  inspect 
for  hole  roundness. 

The  automated  system  can  detect  delamination  and  measling.  if  the  defect  is 
severe  enough  to  be  high  enough  contrast,  and  if  the  defect  is  in  an  area  that 
causes  the  width  of  spacing  to  be  affected.  Conductor  flaws  can  also  be 
easily  detected,  if  they  affect  line  widths.  Conductor  flaws  that  appear  as 
normal  runs  with  breaks  in  them  can  be  missed  by  the  automated  system 
unless  the  system  uses  a  training  board. 

The  defects  on  soldered  boards  were  not  studies  using  this  system  since  the 
illumination  scheme  did  not  give  satisfactory  differentiation  between  soldered 
areas  and  board  material.  Therefore,  defects  such  as  lifted  pads  were  not 
detected,  nor  were  any  of  the  defects  that  appear  only  on  soldered  boards. 

Component  alignment  was  not  readily  detected  since  the  contrast  on  a  com¬ 
ponent  generally  varies  too  greatly  for  the  scene  segmentation  schemes  which 
were  used  to  operate  satisfactorily. 

Scratches,  surface  quality  (within  limits),  and  hole  registration  were  readily 
detected. 

Legibility  of  legends  is  a  nebulous  defect,  but  if  the  legends  are  designed 
to  the  same  specifications  in  line  width  and  spacing  (a  very  reasonable 
approach!)  then  their  inspection  would  be  identical  to  that  of  the  rest  of  the 
board. 
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IV. 


CONCLUSIONS  &  RECOMMENDATIONS 


The  system  as  configured  has  shown  exceptional  resolution  and  good  sensitivity. 
The  use  of  the  automatic  threshold  setting  feature  proved  to  be  less  than  ideal 
except  for  a  given  set  of  lighting  conditions.  To  gain  greater  control  over 
threshold  settings,  a  manual  entry  option  was  added  to  the  PCCARD  Routine 
allowing  entry  of  both  lower  and  upper  thresholds. 

Sufficient  repeatability  was  attained  by  use  of  mechanical  indexing  of  the 
PC  card. 

One  goal  of  this  effort  was  to  determine  how  much  of  the  PC  board  inspection 

task  could  be  done  using  the  computer  as  an  inspector,  so  that  it  could 

inspect  a  board  without  having  seen  a  "training"  board.  The  word  on  this 

particular  technique  was  limited  to  dimensional  analysis  of  a  bare  (unloaded- 

unsoldered)  board.  We  found  that  the  computer  could  readily  inspect  for 

minimum  line  widths  and  spacings,  and  that,  at  our  high  resolution,  it  could 

determine  the  quality  of  the  holes  in  the  etch.  By  measuring  the  ratio  of  the 

square  of  the  perimeter  to  the  area  of  the  holes  out  of  round  holes  or  holes 

with  modules  could  be  detected.  By  pre-selecting  line  spacing  and  line  width 

2 

limitations  and  by  pre-selecting  the  limitation  on  the  "roundness"  (P  /A) 
criterion,  fully  automatic  inspection  is  feasible  for  production  rise.  The 
scanning  of  a  complete  board  and  it's  inspection  for  dimensional  tolerances 
can  be  done  by  the  computer  with  no  need  for  any  training.  This  is  useful 
for  areas  of  the  board  Limited  to  normal  point  to  point  runs  with  round  com¬ 
ponent  holes.  One  interesting  fact  we  discerned  was  that  if  the  lettering  in 
the  legend  area's  of  the  board  was  not  designed  to  the  same  line  width  and 
spacing  as  the  rest  of  the  board,  it  would  be  flagged  as  an  out  of  specification 
area. 
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FIGURE  21  SYSTEM  CONFIGURATION 
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FIGURE  20-B  MEASUREMENTS  FOR  BOARD  SHOWN  IN  FIGURE  20-A 
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FILENAME  IS  :  PCGAP 


CONTOUR  THRESHOLDS  BETWEEN  60  AND  99 

AREA  :  1  BRIGHTNESS  BETWEEN  0  AND  102 

PERIMETER3  219.  19  INTERNAL  AREA3  584.  38 

MINIMUM  INTERNAL  THICKNESS3  4.  19 
IMIN=  1  IMAX=  79  JMIN3  4  JMAX=  26 

AREA  CONTAINED  0  HOLES  AND  0  SPOTS. 

AREA  :  2  BRIGHTNESS  BETWEEN  0  AND  102 

PERIMETER3  46.  51  INTERNAL  AREA3  77.  50 

MINIMUM  INTERNAL  THICKNESS3  4.  00 
I MIN3  86  I MAX3  100  JMIN3  13  UMAX3  22 

AREA  CONTAINED  0  HOLES  AND  0  SPOTS. 

HOLE  IN  AREA  3 

MINIMUM  THICKNESS  AROUND  HOLE  IS  :  1.  00 

HOLE  CENTER-  47,  40 

AVERAGE  RADIUS3  9.  47 
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FIGURE  19B  DATA  FOR  FIGURE  19A 
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FIGURE  17  SCRATCHED  BOARD 


i 


31 


Figure  A-4  Illumination  Specular  &  Diffuse 


THU-  I;  A  DESCRIP  f  ION  OF  THE  SOFTWARE  ON  THE 
PRINTED  CIRCUIT  BOARD  INSPECTION  PR Ob RAM.  (W0-4502) 


THE  F 01  LOW I NO  PROGRAM  ii  USED  TO  GATHER  IMAGES  FROM  THE 
TELEVISION  DIGITIZER  SYSTEM. 

SC  AN  2 

THIS  PROGRAM  SChNS  A  100  EY  100  BLOCK  OF  POINTS 
CENTERED  IN  THE  TELEVISION  PICTURE  THE  PROGRAM  REQUESTS  A 
FILENAME  FOR  STORAGE  ON  THE  DISK, AND  STORES  A  HISTOGRAM  OF 
THE  IMAGE  IN  THE  FIRST  BLOCK  OF  THE  FILE.  THE  DATA  IS  PACKED 
TWO  BYTES  PER  WORD  IN  THE  FILE 

THE  FOLLOWING  PROGRAMS  ARE  USED  TO  DISPLAY  THE  DATA  ON  THE 
TEKTRONIX  TERMINAL  OR  THE  LINE  PRINTER 

H I STOGRAM 

THIS  PROGRAM  DISPLAYS  THE  HISTOGRAM  ON  THE  SCREEN  WITH 
A  SCALE  AT  THE  BOTTOM  WHEN  CALLED,  IT  REQUESTS  A  FILENAME. 

DOTMAT 

THIS  PROGRAM  DISPLAYS  THE  DATA  ON  THE  TEKTRONIX 
TERMINAL  AS  A  CARTOON  THE  PROGRAM  REQUESTS  A  FILENAME  AND 
THEN  REQUESTS  THE  LOWER  AND  UPPER  LIMITS  OF  THE  VALUES  TO  BE 
DISPLAYED. 

PRTMAT 

THIS  PROGRAM  IS  SIMILAR  TO  THE  ABOVE  DOTMAT,  EXCEPT 
IT  DISPLAYS  THE  DATA  ON  THE  LINE  PRINTER. 

PRTHIST 

THIS  PRO OR AM  IS  SIMILAR  TO  THE  ABOVE  "HISTOGRAM",  EXCEPT 
IT  PRINTS  THE  DATA  ON  THE  LINE  PRINTER. 

PRTSUM 

THIS  PROGRAM  PRINTS  THE  SAME  DATA  AS  " PRTHIST" ,  EXCEPT 
THE  DATA  IS  IN  TABULAR  FORM.  THIS  FORMAT  USES  ONLY  ONE  PAGE  ON  THE 
PRINTER 

THE  ABOVE  PROGRAMS  REQUIRE  IMAGE  FILES  THAT  WERE  SCANNED  USING  "SCAN2". 
[THESE  PROGRAMS  WERE  WRITTEN  BY  BOB  .JONES,  AND  USE  THE  FOLLOWING 
SUBROUTINES.  1 

THE  FOLLOWING  SUBROUTINES  ARE  USED  IN  THE  DETAILED  DATA  HANDLING. 
SPOT 

THIS  PROGRAM  IS  USED  TO  INPUT  DATA  FROM  THE  CAMERA. 

LOG  I C 

THIS  PROGRAM  IS  USED  TO  PACT-.  AND  UNPACK  THE  DATA  IN  ORDER 
TO  CONSERVE  STORAGE  SPACE. 

[THE  ABOVE  MACHINE  LANGUAGE  SUBROUTINES  (AND  MANY  EARLIER  EXPERIMENTAL 
PROGRAMS)  BY  TOMMY  REYNOLDS.  J 


B-l 


INC 


THE  FALLOWING  PROGRAMS  WERE  PREPARED  FOR  SCI  SYSTEMS, 

BY  BILL  POPE,  TELCOM  DATA  CORPORATION. 

THERE  ARE  3  ENHANCED  PCEOARD  INSPECTION  PROGRAMS,  BOARD,  CARD,  AND  PCCARB. 

ALL  PROGRAMS  ARE  SIMILAR  IN  FUNCTION.  THEY  EXAMINE  A  100X100  POINT 
SCAN  OF  A  PCEOARD  AND  DETERMINE  VARIOUS  ITEMS  CONCERNING  THE  GEOMETRY 
OF  THE  BOARD  THESE  INCLUDE  OUTLINING  ALL  RUNS  ON  THE  BOARD,  DETERMINING 
THE  SIZE  AND  LOCATION  OF  THE  RIJNS,  PROVIDING  CALCULATIONS  FOR  ANY  HOLES 
FOUND.  AND  EXAMINING  THE  MINIMUM  DISTANCES  BETWEEN  RUNS. 


BOARD 

BOARD  WILL  AUTOMATICALLY  TAKE  A  SERIES  OF  LIVE  SCANS  OF  A  PC  BOARD. 
ALTERNATELY  IT  WILL  TAKE  A  SINGLE  PRE-SCANNED  FILE.  FOR  EACH  SCAN. 
IMAGE  THRESHOLDS  ARE  DETERMINED  BY  LOOKING  AT  THE  SCAN  DATA  AS  IF  IT 
WERE  A  CONTOUR  MAP.  AS  THE  IMAGE  MOVES  BETWEEN  LEVELS,  THE  AVERAGE 
VALUES  BETWEEN  LEVELS  IS  DETERMINED  AS  THE  THRESHOLD  SEPERATING 
THE  BOARD  FROM  THE  RUN.  THESE  THRESHOLDS  ARE  USED  FOR  RUN 
DETERMINATION  AND  HOLE  DEFINITION.  FINALLY  CRITICAL  DISTANCES 
BETWEEN  RUNS  ARE  EXAMINED  BY  TRACING  A  CRITICAL  DISTANCE  BORDER 
AROUND  EACH  RUN. 

SUBROUTINES  USED  BY  BOARD  ARE. 

BOARD,  AEDGE,  TSPOT,  LOGIC,  UNPACK,  SMOOTH,  TRACE,  FILLIN, 

INTERNAL,  DISTANCE,  TRACK,  CKPT,  CONVAL,  SPDOT,  SPLIN,  SPMOV, 

HOLE, SPDQTLIN 


PCCARD 

PCCARD  IS  ALMOST  THE  SAME  AS  BOARD  WITH  THESE  EXCEPTIONS.  FOR  LIVE 
SCANS,  PCCARD  PROVIDES  ONLY  A  SINGLE  SCAN.  THE  MINIMUM  DISTANCE 
FROM  EACH  RUN  TO  ANY  OTHER  RUN  IS  DETERMINED  AS  OPPOSED  TO  JUST 
LOOKING  AT  CRITICAL  DISTANCES. 

SUBROUTINES  USED  ARE. 

PCCARD,  AEDGE,  TRACE,  SEGMENT.  LSCAN,  SPOT,  UNPACK,  SMOOTH,  FILLIN, 
HOLE,  INTERNAL,  SVDISTANCE,  GUESS,  CONHIST,  LOGIC,  SPDOT,  SPLIN. 
SPMOV,  SPDOTL IN 

THE  FOLLOWING  SUBROUTINES  ARE  USED  TO  DISPLAY  GRAPHICS  ON 
THE  TEKTRONIX  4006- 1  TERMINAL 

SPLIN  (IXPOS,  IYPOS,  IXEND,  IYEND)  THIS  SUBROUTINE  IS  USED 
TO  PLOT  LINES  FROM  LOCATION  ( IXPOS, IYPOS)  TO  ( I XEND,  I YEND ) . 

SPMOV  ( IXPOS,  IYPOS)  THIS  SUBROUTINE  IS  USED  TO  MOVE 
THE  POSITION  OF  THE  NEXT  CHARACTER  TO  BE  PRINTED  TO  ( I XPOS,  IYPOS ) . 

SPDOT  ( IXPOS, IYPOS)  THIS  SUBROUTINE  IS  USED  TO  PLOT 
A  DOT  AT  POSITION  ( IXPOS, IYPOS). 


C THESE  GRAPHICS  PROGRAMS  WERE  WRITTEN  BY  PETER  P.  PRYOR  JR.  ] 


CARD 


CARD  WAS  THE  INITIAL  DEVELOPMENT  PROGRAM.  IT  ASLO  IS  SIMILAR  TO  THE 
OTHER  TWO  PROGRAMS  IT  USES  ONLY  PRESCANNED  FILES  CREATED  BY  PROGRAM 
SCANS'  THE  METHOD  OF  DETERMINING  THRESHOLDS  IS  INSPECTION  OF  THE 
HISTOGRAM  DATA  THRESHOLDS  MUST  EE  ENTERED  MANUALLY  AND  PLOTTING 
OF  THE  HISTOGRAM  DATA  IS  OPTIONAL. 

SUBROUTINES  USED  ARE. 

CARD,  TRACE,  SEGMENT,  UNPACK,  PEAKS,  SMOOTH,  PLHIST,  EDGE,  FILLIN, 
INTERNAL,  SVDISTANCE,  GUESS,  CONHIST,  LOGIC,  SPDOT,  SPLIN,  SPMOV, 
HOLE, SPDOTLIN 


TWO  OTHER  PROGRAMS,  BOTH  INCONSEQUENTIAL,  WERE  LEFT  ON  THE  FILE.  CONTOUR  AND 
BVIEW  PROVIDE  A  3D  CONTOUR  PLOT  AND  BACK-VIEW  OF  PRESCANNED  DATA  FILES. 
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PROGRAM  PC CARD 


PREPARED  FOR  SCI  SYSTEMS,  INC. 

BY  BILL  POPE,  TELCOM  DATA  CORPORATION. 

THIS  SUBROUTINE  IS  THE  MAIN  DRIVER  PROGRAM  WHICH  EXAMINES  SCANS 
OF  PC  BOARDS. 


COMMON  ARRAY, SUM, LP,  IMIN,  IMAX,  JMIN, JMAX, NCNT , NAREA, AMAX 
I NTEGER  ARRAY ( 1 00,  1 00 ) ,  SUM ( 256 ) ,  AMAX  <  5, 20 ) 

EXTERNAL  OVERO, 0VER1, OVERS, 0VER6, 0VER7 
C 

C  INITIALIZE  OVERLAYS 

C 

CALL  OVOPN  ( 5,  "PCCARD.  OL",  I  ERR) 

IF  (  I  ERR.  EQ.  1  )  GO  TO  380 

TYPE  "  OVERLAY  OPENING  ERROR  ",  I ERR 

STOP 

SO  CONTINUE 


LP=10 

TYPE  "  PRINT  OR  TYPE  RESULTS?  P, T  " 
READ! 1 1 , 100) NANS 
IF  (NANS.  NE.  "P"  >G0  TO  7 
LP=  1 2 

CALL  F0PEN( 12,  "*LPT"  ) 

CONTINUE 


FIND  OUT  ABOUT  LIVE  SCAN 

0  TYPE  "  LIVE  SCAN?  Y.  N" 

READ< 11, 100) NANS 
00  FORMAT (SI) 

IF  (NANS.  EQ.  "N"  )G0  TO  20 
IF  (NANS.  NE.  "Y"  >G0  TO  10 

LOAD  SCAN  OVERLAY 

CALL  0VL0D(5, 0VER7, 0, IERR ) 
IF  ( I  ERR.  NE.  1  )  GO  TO  980 
C 

CALL  LSCAN 
C 


B-4 


go  r  o 


CONTINUE 

INITIALIZE  AND  UNPACK  DATA 

CALL  OVLOD  (5,  OVERO.  O.  I  ERR  ) 

I F  i  I EF:R.  NE  1  'GO  TO  980 

CALL  UNPACK 

CONTINUE 

CALL  OVLOD < 5 . OVER 1 . 0, I ERR ) 

IF(  I  ERR.  NE.  1  )  GO  TO  980 

SMOOTH  MAKES  EACH  POINT  EQUAL  THE  AVERAGE  OF  ITS  NEIGHBORS 


CALL  SMOOTH 

PEAK  TAKES  HISTOGRAM  IN  SUM  AND  DETERMINES  PEAKS  FIRST  MAJOR 
PEAK  FOUND  LIMITS  ARE  RETURNED  IN  LTH  AND  UTH. 

CALL  PEAKS 

TEST  THRESHOLD 

CALL  OVLOD <  5, OVER*. 0-  I ERR ) 

IF  (  I  ERR.  NE.  1)00  TO  980 

CALL  CONHIST ( ARRAY, LP. LTL, LTU ) 


CALL  EDGE  TO  BREAK  ARRAY  INTO  AREAS  WITHIN  LTH,  UTH 
CALL  AEDGE 1 LTL, LTU) 
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980 


CALL  0VL0D(5- OVERS, 0, I ERR) 
IF  <  I  ERR  NE.  1  'GO  TO  930 


NOW  TRY  TO  FIND  DISTANCES  BETWEEN  AREA 


CALL  DISTANCE 


CALL 

CALL 

STOP 

TYPE 

STOP 

END 


CLOSE (5, I ERR) 

RESET 

"  ERROR  LOADING  OVERLAY 


",  I  ERR 


T I TL  LOGIC 
DOC  FORTRAN  LOGIC 
FUNCTIONS 

ENT  I AND,  I  OR.  I LEFT 


ENT  IRIGHT 

.  EXTN  FRET 
.  EXTD  CPYLS 

VAL  =  -167 

ARG1  =  VAL+1 

ARG2  =  ARG 1+1 

SIZE  =  2 

NREL 


I  AND  (  I.  J> 
SIZE 


I  AND: 

■JSR 

CPYLS 

LDA 

0  SARG1  3 

LDA 

1  <2  ARG  2  3 

AND 

1  0 

STA 

0  @VAL  3 

FRET 

i 

I OR i I. J> 

.  SIZE 

I  OR: 

JSR 

&  CPYLS 

LDA 

0  @ARG 1  3 

LDA 

1  0ARG2  3 

COM 

1  1 

AND 

1  0 

ADC 

1  0 

STA 

0  <2VAL  3 

FRET 

i 

I  LEFT  U. 

NUMBER )  , 

i 

IRIGHT ( I, NUMBER) 

.  SIZE 

ILEFT 

SUB 

2  2 

JMP 

COMMN 

SIZE 

IRIGHT. 

ADC 

2  2 

COMMN : 

JSR 

<1  CPYLS 

LDA 

0  @APG 1  3 

LDA 

1  SARG2  3 

NEC 

1  1  SNR 

JMP 

OUT 

FUNCTION  VALUE  RETURNED 
PARAMETER  1 
PARAMETER  2 
FRAME  STACK  SIZE 
NORMAL  MEMORY 


COPY  ARGUMENT  ADDRESSES 
GET  PARAMETER  1 
GET  PARAMETER  2 
PERFORM  LOGICAL  AND 
SET  FUNCTION  VALUE 
RETURN  TO  CALLER 


COPY  ARGUMENT  ADDRESSES 

GET  VARIABLE 

GET  PARAMETER  2 

FORM  1@S  COML 

MASK  OUT  BITS 

INSERT  BITS 

SET  RETURN  VALUE 

RETURN  TO  CALLER 


;  INDICATE  LEFT 
-  GO  TO  COMMON  ROUTINE 

i 

•  INDICATE  RIGHT 
.  COPY  ARGUMENT  ADDRESSES 
i  GET  PARAMETER  1 
.  GET  SHIFT  COUNT 
;  NEGATE  COUNT 
,  RETURN  IF  ZERO 
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L  OOF- 


MOV 

2  2  SNP 

MOV  1 L 

0  C>  SKF‘ 

MOVZR 

i  <  .I 

INC 

1  1  SIR 

..IMF' 

LOOP 

STA 

0  eVAL 

FRET 

END 


E-KIP  IF  RIGHT 
E.HIFT  LEFT 
-;HIFT  RIGHT 
E-KIP  IF  DONE 
z#  I  p  SOME  MORE 
E.ET  OUTPUT  VALUE 


__  o 


E 


L  Tufr  -L  T  H+NB  I  m‘£- 
•JTHF  1 U  r  H  ^  N E*  I 

rn  C'  J,  :t  t  i  i  t 

jmini=jnin+i 
I  F  .  JM  I  i  i  1  ijT.  iJNm x  >  RETURN 
JO  15u  J=.JMIN1  ■  JMA* 

INITIALIZE  FLmG.  AO  WE  SEARCH  THIS:  PORTION  OF  ARRAY ,  POINTS  MAY 
FALL  INSIDE  THRESHOLD  WHILE  BEING  EXTERNAL  TO  SEGMENT  WE  ARE 

examining  will  set  flag  when  we  find  first  marked  edge  to 

INDICATE  MOVING  INTERNAL  TO  SEGMENT. 

NFL*'.* 

DO  149  I - I M I N • I MAX 
IF  (NFL.  GT  O'  GO  TO  145 

TO  THIS  LEADING  EDGE  OF  SEGMENT. 

I F  <  A  (I ,  .J  ).  GE  NE I  AS  >  NFL*  1 
GO  TO  147 

DETERMINE  IF  POINT  IS  WITHIN  THRESHOLDS. 

IF  (A*  I  -  J)  LT  LTH.  OR.  A  \  I ,  •  J ) .  GT,  UTH)GO  TO  146 


FOUND  A  POINT  LOOK  AT  ADJACENT  POINTS  TO  SEE  IF  ANY 
ARE  OUT  OF  THRESHOLD  LIMITS.  IF  SO.  CONSIDER  AS 
START  OF  HOLE  IN  AREA  GO  TO  TRACE  TO  FIND  EDGE  OF  HOLE. 
LOOKING  AHEAD  TO  NEXT  ROW  AND  COLUMN  TO  "ANTICIPATE"  HOLES. 

N=J 
M*  I  - 1 

IFvM.  LT.  I  GO  TO  12 

IF  (AIM,  Ni.  GE.  LTH.  AND  A  <  M.  N  ) .  LE.  UTH  >  GO  TO  12 
IF l  A(M,  N).  GE.  LTHFlGO  TO  12 
GO  TO  25 

M=I  +  1 

IF(M  GT  1  •.»*>>  GO  TO  14 

IF  (  A  (  M>  N  ) .  GE  LTH.  AND  A(M,  N>.  LE.  UTH)  GO  TO  14 
IF  (  A  <  M,  N ) .  GE.  LTHF  >  GO  TO  14 
GO  TO  25 

M*  I 
N=  J-l 

I F ' N  LT  1)  GO  TO  If 

I F  v  A  •:  M ■  N  >  GE.  LTH  AND.  A  <  M.  N  ;  LE.  ijTH'GO  TO  If 
I F ( A  I M ,  N  >  GE  L  THF ) GO  TO  1 o 
GO  TO  25 


V'-ERI.  A'  u'-'ERh 

SUBROUTINE  F  I LL 1 1-1  •.  L  Th  ■  UTH,  NBIAS,  RaREA.  I  DEBUG,  B,  NEP,  NHOLE,  NSPOT  ) 

r  re  fared  f  op  s c  i  s  i  s  r e m ?■  •  i  n c 
BY  BILL  PORE.  TELCOM  DATA  CORPORATION. 

THIS  SUBROUTINE  WORKS  WITH  AN  AREA  WHOSE  EDGES  HAVE  BEEN  DEFINED  IN 
SUBROUTINE  EDGE.  GIVEN  THAT  THE  AREA  FALLS  BETWEEN  THE  GIVEN 
I  AND  J  MINS  AND  MAX  VALUES.  AND  THAT  THE  EDGE  POINTS  HAVE  BEEN 
GIVEN  A  BIAS  OF  NBIAS,  THIS  SUBROUTINE  FILLS  IN  ALL  THE  POINTS  INTERNA 
TO  THE  AREA  (FALLING,  OF  COURSE,  BETWEEN  THE  SEARCH  LIMITS). 

IT  ALSO  CALCULATES  THE  AREA  FOUND  BY  SUMMING  THE  NUMBER 
OF  POINTS 

WHILE  PERFORMING  THE  INTERNAL  SEARCH  FOR  POINTS  WITHIN  THE  SEGMENT, 

IF  A  POINT  IS  FOUND  WHICH  IS  OUTSIDE  THE  SEARCH  LIMITS  THEN  IT  IS 
ASSUMED  TO  EE  THE  EDGE  OF  AN  INTENAL.  "HOLE"  THE  SUBROUTINE  HOLE 
IS  CALLED  TO  TRACE  THE  HOLE  BOUNDARIES  AND  DETERMINE  HOLE  RELATED 
PARAMETERS. 


VARIABLES  USED  INCLUDE: 

LTH,  IJTH  LOWER  AND  UPPER  THRESHOLDS  FOR  INTERNAL  SEARCH 

NBIAS  THE  VALUE  USED  TO  "MARK"  THE  EDGES  OF  THIS  SEGMENT 

(MARK  BY  ADDING  NBIAS  TO  ARRAY ( I ,  J ) ) 

RAREA  THE  CALCULATED  INTERNAL  AREA  OF  THE  SEGMENT.  VALUE 

IS  REAL  NUMBER  AND  DOES  NOT  INCLUDE  HOLE  AREAS. 

B ( 2, NPT )  INTEGER  ARRAY  OF  2,  NPT  POINT  CONTAINING  THE  I ,  J 

VALUES  OF  THIS  SEGMENTS  EDGE. 

NHOLE  NUMBER  OF  HOLES  FOUND 

NSPOT  NUMBER  OF  SPOTS  FOUND  (HOLE  LESS  THAN  20  PERIMETER). 

IMIN, IMAX  THE  MINIMUM  AND  MAXIMUM  COORDINATES  OF 

JMIN, UMAX  THE  CURRENT  SEGMENT. 


COMMON  ARRAY,  SUM,  LP,  IMIN,  IMAX,  .JMIN,  UMAX,  NCNT ,  NaREA,  AMAX 
COMMON  /PLT/  ISX, ISY, MIDX, MIDY, NXS, NYS 
INTEGER  ARRAY < 100,  100), SUM (256),  AMAX (5, 20) 

INTEGER  Av  100,  luO) ,  UTH,  UTHF 
INTEGER  B <  2,  1500; 

ECU I VALENCE  ( A ( 1 , 1 ) , ARRAY (1,1)) 


NH0LE=0 
NSP0T=0 
I  A=0 

SET  PSUEDO  THRESHOLDS  FOR  POINTS  "MARKED"  BY  THIS  NBIAS 
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SAVE  MINIMUM  AND  MAXIMUM  VALUES  FOR  SUBROUTINE  DISTANCE. 


AMAX' 1, NON T ) = I M I N 
AM AX « 2. NCNT > = IMAX 
AMAX  v 3.  NCNT)=JMIN 
AMAX  (4,  NCNT >  =JMAX 
AMAx (5, NCNT ) =NHO 
CONTINUE 
CONTINUE 
MCNT =NCNT 
RETURN 

TYPE  "  SECOND  LEVEL  OVERLAY  LOAD  ERROR" 

STOP 

END 


o  o 


NOW  CALL  'SUBF.OUT  INE  TRACE  TO  TRACE  OUT  THE  EDGE 


C 

c 

r 
C 
0 

NSPL®25c- 

C  MARK  A'E  GUTTER  EDGE  OF  APEh 

C 

CALL  TRACE ( I ROW. JCOL. 10. JO, LTH, UTH, PERIM, AREA, NSPL, 

1  I  MIN,  IMAX.  .JMIN.  UMAX.  B-  NPT  > 

IF  (PERIM.  NE.  U.  0  i  GO  TO  77 
C 

C  ONLY  A  SINGLE  POINT  SO  BACK  OFF  AREA. 

AREA=AREA-NAREA 
GO  TO  73 
77  CONTINUE 

C 

c 

C  NOW  WE  HAVE  THE  BOUNDARY,  FILL  IN  THE  REST 

C 

C 

c 

c 

NCNT=NCNT+1 

CALL  FILL  IN  ( LTH,  UTH,  AREA,  RAREA,  I  DEBUG.  B,  NPT,  NHO,  NSF'T ) 

C 

C 

C  GO  LOOK  AT  AREAS  INTERNAL  DIMENSIONS 

C 


CALL  0VL0D(5, OVERB, 0, IER ) 
IF (  IER.  NE.  1  >G0  TO  330 
CALL  INTERNAL (B. NPT, XDIM) 


WRITE(LP,  260  )NCNT ,  LTH,  UTH 

260  FORMAT ( "OAREA  .  ",  13,  "  BRIGHTNESS  BETWEEN  ",I4,  "  AND  ",I4) 

WRITE  < LP- 261 > PERIM, RAREA 

261  FORMAT (  "  PERIMETER®  ",  F3.  2,  5X,  "  INTERNAL  AREA®  ",F8.  2) 

WRITE ( LP' 162) XDIM 

162  FORMAT < "  MINIMUM  INTERNAL  THICKNESS®  ",Ffc.  2) 

WRITEvLP, 262) IMIN. IMAX, JMIN. UMAX 

262  FORMAT ( "  IMIN=",I4, "  IMAX®" , 14, "  JMIN®", 14, "  JMAX®",I4) 

WR I TE • LR , 263 ) NHO ■ NSPT 

263  FORMAT ( "  AREA  CONTAINED  ",  13,  M  HOLES  AND  ",  13,  "  SPOTS.  " > 
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T»  w 


T 


tr 


C  PEAKS  FOUND  FROM  EXAMINATION  OF  THE  "HISTOGRAM".  WHEN 

C  SEARCHING  FOR  THE  EDGES  OF  AREAS,  A  FUDGE  FACTOR  WILL  BE 

C  USED  TO  EXTEND  THE  SEARCH  RANGE.  THE  POINTS  WITHIN  THE 

C  EXTENDED  RANGE  WILL  ONLY  BE  ACCEPTED  IF  THEY  ARE  ADJACENT 

C  TO  A  POINT  WITHIN  THE  ORIGINAL  SEARCH  RANGE. 

C 

c 

ND=(UTH1-LTH1 )/10 

LTH=LTH1-ND 

LTH=0 

IF  ( LTH.  LT.  0 )  LTH— 0 

IJTH=UTH  1 +ND 

IF < UTH.  GT.  255)  UTH-255 


SIZE  FACTORS  TO  ALLOW  FOR  THE  FACT  SCAN  IS  NOT  SQUARE. 

NXSI ZE=5 
NYSI ZE=4 
MIDX-512 
MIDY-370 

ISTARTX=MIDX- ( ( NXSI ZE* 100 ) /2 ) -NXSI ZE 
I STARTY=M I DY+ ( ( NYS I ZE* 1 00 ) /2 ) +NYS I ZE 
ISTARTX=200 
I ST ARTY-620 

LOOP  THROUGH  ENTIRE  ARRAY 

DO  79  JC= 1 , 100 
DO  78  IR  =  1, 100 
JCOL-JC 
IROW=IR 

IS  DATA  POINT  WITHIN  ORIGINAL  THRESHOLDS 

IF<A<  IROW,  JCOL).  LT.  LTH1.  OR.  A(  IROW,  JCOL).  GT.  UTH1  )  GO  TO  78 

YES,  FOUND  FIRST  POINT  OF  NEW  AREA 
BACK  UP  ON  I  FOR  EXTENDED  SEARCH  RANGE 

IF <  IROW.  EQ.  1  )G0  TO  10 
I - I ROW— 1 

IF ( A (  I.  JCOL).  LT.  LTH.  OR.  At  I,  JCOL).  GT.  UTH)GO  TO  10 
IROW=I 
GO  TO  9 
10  CONTINUE 

I 0=1 ROW- 1 
JO— JCOL 

AREA=AREA+NAREA 

C 

CALL  0VL0D(5, OVER A, O, IER) 

IFtlER.  NE.  1)00  TO  880 
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SUBROUTINE  AEDGEiLTHl,  UTH1 > 

PREPARED  FOR  SCI  SYSTEMS-  INC. 

BY  BILL  POPE,  TELCOM  DATA  CORPORATION. 

THIS  SUBROUTINE  EXAMINES  THE  SCAN  DATA  FOUND  IN  ARRAY  A. 

ALL  POINTS  WHOSE  VALUE  IS  BETWEEN  THE  LOWER  THRESHOLD  LTH 
AND  THE  UPPER  THRESHOLD  UTH  ARE  DIVIDED  INTO  SEGMENTS. 

A  SEGMENT  IS  AN  AREA  OF  POINTS  CONTIGUOUS  TO  EACH  OTHER  BUT 
SEPERATED  FROM  ALL  OTHER  POINTS  WITHIN  THE  THRESHOLD  LIMITS. 


A  IS  AN  ARRAY  DIMENSIONED  100  BY  100  CONTAINING  SCAN  DATA  POINTS. 
LTH1  IS  THE  LOWER  BRIGHTNESS  THRESHOLD 
IJTH1  IS  THE  UPPER  BRIGHTNESS  THRESHOLD 
NCNT  IS  THE  NUMBER  OF  AREAS  FOUND. 


COMMON  ARRAY,  SUM,  LP,  IMIN,  IMAX,  -JMIN,  JMAx,  NCNT,  NAREA,  AMAX 
COMMON  ,/PLT  / 1  ST  ART  X ,  ISTARTY,  MIDX,  MIDY,  NXSIZE,  NYSIZE 
INTEGER  ARRAY ( 100, 100 ) ,  SUM ( 256 > , AMAx (5, 20 ) 

I NTEGER  B  <  2,  1 500 ) , SEG  <  3,  20 ) 

INTEGER  UTHF , UTH 1 

INTEGER  A(100> 100),  AREA,  UTH 

EQUIVALENCE  ( Av 1. 1 ) , ARRAY (1,1)) 

EXTERNAL  OVERA, OVERB 


OTHER  VARIABLES  USED  ARE. 

MI  -  MAXIMUM  I  'ROW)  WITHIN  CURRENT  COLUMN  WHICH  IS  PART  OF  AREA 
LI  -  LOWEST  I  (  ROW)  WITHIN  CURRENT  COLUMN 
NLI  -  LOWEST  I  IN  LAST  COLUMN 

IS  -  START  I  VALUE  FOR  EXAMINING  CURRENT  COLUMN 

INITIALIZE  DATA.  NAREA  IS  SET  TO  VALUES  THAT  WILL 
UTILIZE  THE  HIGH  ORDER  BITS  OF  A  SCAN  POINT  AS  A  FLAG 
INDICATING  THE  AREA  NUMBER  TO  WHICH  IT  BELONGS. 

AREA=0 
NCNT =0 
NAREA* 1024 
I DEBUG* 1 
ISP=2S8 
IUS«287 


IF ( (LTH1+UTH1 ).  EQ  0) RETURN 


LTH1  AND  UTH 1  REPRESENT  THE  NOMINAL  DIVISION  VALUES  BETWEEN 
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o  o 


4" 


415 

C 


NBC  T  =  - 1 
IB=" 

IT-0 

DO  14  1  =  1'  LRENO 

FI  ROT  i_.  HECK  NEW  POINT  FOR  MAXIMUM  BRIGHTNESS 

IF ■ ARRAY > I, J> .  LT.  MIN  > MIN=ARRAY (I, J) 

NS=ARRA  Y  ■  1 ,  J )  -ARF'A  y  (  I +STP .  J ) 

•JUMP  DEPENDENT  ON  WHICH  REGION  WE  HAVE  SEEN  IN. 

IF <  SLOPE  >30.  40,  50 

HAVE  BEEN  SEARCHING  ALONG  LEVEL  (SLOPE  =  •.*). 

IF  THE  ABSOLUTE  VALUE  OF  THE  SLOPE  EXCEEDS  THE  LIMIT  NSTP, 
THEN  CONSIDER  STARTING  MOVE  TO  NEXT  LEVEL. 


SAVE  THE  PRESENT  VALUE  IN  NTQP  OR  NBOT  DEPENDENT  ON  WHICH 
DIRECTION  WE  ARE  STARTING  TO  MOVE. 

IF( I ABS i NS ’ .  LT  NSTP  > GO  TO  24 
I F  (  NS.  GT.  0 )  GO  TO  45 
NTOP=ARRAY (  I ,  J ) 

GO  TO  20 

NBOT  ® ARRAY « I , J  > 

GO  TO  20 


C  AT  LABELS  30  AND  50,  WE  HAVE  BEEN  MOVING  BETWEEN  LEVELS.  WHEN  THE 

C  VALUE  OF  THE  SLOPE  FALLS  BELOW  THE  LIMIT  NSTP,  THINGS  ARE  "LEVELING 

OUT.  IF  THERE  WAS  A  PREVIOUS  LEVEL  (NOT  STARTING  A  SIDE  OF  SCAN), 
THEM  GET  AVERAGE  BRIGHTNESS  AND  SAVE. 

C 

30  IF  (NS.  LT.  -NSTP)  GO  TO  24 
IF < NTOP.  LT  0 > GO  TO  31 

AVG=AVG+ ( FLOAT ( NTOP+ ARRAY l I , J ) > ) /2. 

CNT=CNT+1 

31  NTOP=- 1 
NE0T=-1 
GO  TO  20 

5*!>  IF  (NS.  GT  NSTP)  GO  TO  24 

IF ( NBOT  LT  0 ) GO  TO  31 
AVG=AVG+ <  FLOAT ( APRA  Y  c I , J ) +NBOT )  ) /2. 

CNT =CNT  + 1 
GO  TO  31 
20  SLOPE=NS 

I F ( I ABS ( NS )  LT.  NSTP > SL0PE=0 
24  CONTINUE 
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OVERLAY  OVER* 

SUBROUTINE  CONHIST  •;  ARRAY  .  LP,  MIN,  LT ) 

PREPARED  FOR  SCI  SYSTEMS,  INC. 

BY  BILL  POPE,  TELCOM  DATA  CORPORATION. 


THIS  SUBROUTINE  TRIES  A  DIFFERENT  METHOD  FOR  THRESHOLDS. 

ASSUMMING  THAT  THE  DATA  ARRAY  RECEIVED  FROM  A  SCAN  BASICALLY 
FORMS  A  CONTOUR  OF  THE  BRIGHTNESSES,  CONHIST  SEARCHES  ONE  ROW 
AT  A  TIME  LOOKING  FOR  THE  CHANGES  IN  CONTOUR.  THE  CHANGING  CONTOUR 
CAN  THEN  BEE  SEEN  AS  LEVELS  REPRESENTING  THE  BOARD,  THE  RUNS,  AND 
THE  "BOTTOM"  OF  HOLES.  BY  FINDING  THE  AVERAGE  MIDPOINT  BETWEEN  LEVELS 
EACH  TIME  THE  CONTOUR  MOVES,  THESE  MIDPOINTS  CAN  THEN  BE 
AVERAGED  TO  FIND  THE  "UPPER",  OR  DARKER,  THRESHOLD.  AT  THE 
SAME  TIME,  THE  BRIGHTEST  SPOT  ON  THE  BOARD  CAN  BE  FOUND  TO  USE 
AS  THE  LOWER  THRESHOLD. 


INTEGER  ARRAY ( 100,  100),  STP,  SLOPE 


THE  METHOD  FOR  DETECTING  MOVEMENT  BETWEEN  LEVELS  IS  EXAMINATION  OF  THE 
SLOPE  OF  THE  CONTOUR.  SINCE  THE  "X"  COORDINATE  WILL  BE  CONSTANT  FOR 
EACH  COMPARISON.  ONLY  THE  "Y"  COORDINATE  (DIFFERENCE  IN  BRIGHTNESS) 
WILL  BE  EXAMINED. 


INITIALIZE  VALUES 

NSTF-10 

STP=1 

LPEND= 1 00— STP 
MIN=255 
CNT =0. 

AVG=0. 

DO  25  J=5,  100,  5 

AS  WE  EXAMINE  EACH  ROW,  NS  REPRESENTS  THE  NEW  SLOPE  FOR  THE  NEXT 
POINT  TO  EXAMINE  WHILE  SLOPE  IS  THE  GENERAL  SLOPE  IN  THE  REGION 
WE  ARE  SEARCHING.  A  SLOPE  OF  ZERO  INDICATES  WE  ARE  ON  THE  BOARD, 

A  NEGATIVE  SLOPE  INDICATES  MOVEMENT  TOWARDS  A  BRIGHTER  LEVEL,  AND 
A  POSITIVE  SLOPE  TOWARD  A  DARKER  LEVEL. 


NS=ARRAY  < 1 , J  > -ARRAY ( STP+ 1 ,  J ) 
SLOPE=NS 

IF ( IABS(NS).  LT.  NSTP)SL0PE=0 
NT0P=-1 


OVERLAY’  OVER1 
SUBROUTINE  SMOOTH 
C 

C  THIS  SUBROUTINE  TAKES  THE  1UO  By  100  SCAN  DATA  IN  ARRAY 

C  AND  SMOOTHS  IT  BY  REVALUING  EACH  POINT  TO  THE  AVERAGE 

C  OF  IT'S  EIGHT  ADJACENT  NEIGHBORS.  A  NEW  256  VALUE 

C  "HISTOGRAM"  IS  THEM  CONSTRUCTED  IN  SUM. 


C 


C 


100 

C 

C 

c 

9 

C 

c 

c 

c 


17 

27 


13 

28 

C 

c 

c 


COMMON  ARRAY, SUM, LP, IMIN, IMAX, JMIN, UMAX ■ NCNT ■ NAREA, AMAX 
INTEGER  AMAX (5, 20) 

I NTEGER  ARRAY ( 1 00#  1 00 ) ,  SUM  t  256 ) 

TYPE  "  SMOOTH  DATA?  Y.  N“ 

READ  < 11,  100) NANS 
FORMAT (SI > 

IF(NANS.  EQ.  "N"  ) RETURN 

CLEAR  OLD  "HISTOGRAM"  AS  NEW  ONE  WILL  BE  FORMED. 

DO  9  1  =  1,256 
SUM ( I ) =0 


AVERAGE  ARRAY  POINTS. 

DO  11  J= 1,100 
DO  10  1=2,99 
I SUM=0 

I SUM=ARRAY \ 1  +  1, J ) + ARRAY ( I - 1 ,  J ) 

NN= J- 1 
NS=I-1 
NE=I+1 

I F  ( J.  EQ.  1  )G0  TO  27 
DO  17  K=NS> NE 
ISUM=ISUM+ARRAY<K,  NN) 

CONTINUE 
CONTINUE 
NN=.J+ 1 

IF  ( J.  EQ.  100) GO  TO  28 
DO  18  K=NS> NE 
I SUM= I SUM+ARRA Y < K , NN > 

CONTINUE 

ND=3 

I F  ( J.  EQ.  1.  OR.  J.  EQ.  1 00 )  ND=5 
ARRAY ( I , J ) = I SUM/ND 

NOW  USE  NEW  VALUE  FOR  "HISTOGRAM". 

N=ARRAY < I , J ) + 1 
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LOW  +  1 


LOW  = 

M  =  0 

IF  <  DELTAY.  LT.  STEF'2)  GO  TO  11 

DELTAY  =  0 

LOW  =  LOWERY 

M  =  0 

CONTINUE 

END  OF  SCAN.  MAKE  SURE  ALL  POINTS  PUT  IN  ARRAY 


IF<  I  A.  GT.  100) RETURN 
I F  <  K.  EQ.  0  >  RETURN 


DO  50  1  =  1,  K 

ARRAY •  I A, JA  ) =POINTS t 3,  I > 
JA= JA+2 

IF ( JA.  LE.  100) GO  TO  50 
JA=JA— 99 

I F  ( -JA.  EQ.  2)  GO  TO  50 

•JA=  1 

IA=IA+1 

CONTINUE 


RETURN 


L  OOP 

Tu  CPEm 

IE  5<-’0  POINTS  rO  SPOT 

r  a=  t 

“ 

JA-- 1 

2' 

DO  10 

COLUMN 

=  LEFTX,  RIGHTX, STEP 

1 1 

DO  20 

ROW  = 

LOW.  UPPER Y.  STEP2 

K  =  k 

1 

POINT 

5  <  1 .  K  < 

=  COLUMN 

POINT 

S  v  2  •  K  > 

=  ROW 

M  =  M 

+  1 

IF  f  K 

EQ.  500 

)  GO  TO  30 

GO  TO 

20 

r 

c 

CALL 

SPOT 

30 

CALL 

SPOT'K. 

POINTS) 

c 

c 

c 

c 

c 

c 

e  put  data  into  core  array 

c 

c 

DO  45  1=1. 500 

ARRA i < I A , UA >  =P0 1 NTS (3.  I ) 

JA=JA+2 

IF(.JA  LE.  1 00  >  GO  TO  45 
UA= JA-99 

IF  OJA.  EQ.  2  >00  TO  45 

•JA=  1 

IA=IA+1 

45  CONTINUE 
C 
C 
C 

c 

c 

K  =  0 
GO  TO  21 

20  CONTINUE 
C 

c  note  that  there  is  a  double  do  loop  using  STATEMENT  NUMBERS  10  AND  20 

C  CARE  SHOULD  BE  TAKEN  IN  FOLLOWING  THE  CORRECT  PATH 

c 

21  DELTAY  =  DELTAY  +  STEP 
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OVERLAY  0“ER7 
SUBROUT  I WE  LSCAN 

PREPARED  FOP  SCI  SYSTEMS.  INC. 

BY  BILL  POPE.  TELCQM  DATA  CORPORATION. 


THIS  PROGRAM  IS  NAMED  LSCAN 

BASICALLY  AN  ARRAY  CALLED  POINTS  IS  CREATED  THRU  A  DO  LOOP 
AND  PASSED  TO  AN  ASSEMBLY  LANGUAGE  SUBROUTINE  CALLED  SPOT 
WHICH  DOES  THE  ACTUAL  SCANNING.  THE  THIRD  ARGUMENT  OF  THE 
THREE  DIMENSIONAL  ARRAV  CONTAINS  THE  RETURNED  DIGITIZED 
VALUE  WHICH  CAN  VARY  FROM  0  TO  255.  THE  500  VALUES 
RETURNED  ARE  PLACED  IN  THE  DATA  ARRAY (1 00, 100 ).  WHEN 
IT  IS  FULL,  CONTROL  IS  RETURNED  TO  THE  DRIVER. 


LEFTX  IS  THE  X  OF  THE  BEGINNING  <  X , Y )  POSITION 
RIGHTX  IS  THE  X  OF  THE  LAST  (X,Y>  POSITION 

UPPERY  IS  THE  X  OF  THE  BEGINNING  <X,Y>  POSITION 

LOWER Y  IS  THE  Y  OF  THE  LAST  (X,Y)  POSITION 

ROW  IS  VARIABLE  CONTAINING  CURRENT  ROW  BEING  PROCESSED 
COLUMN  IS  VARIABLE  CONTAINING  CURRENT  COLUMN  BEING  PROCESSED 
STEP  IS  VARIABLE  ALLOWING  DIFFERENT  SCAN  METHODS 

POINTS  IS  AN  ARRAY  WHICH  IS  PASSED  TO  AN  ASSEMBLY  LANGUAGE 
SUBROUTINE  WHICH  DOES  THE  ACTUAL  SCANNING. 

COMMON  ARRAY,  SUM,  LP-  I  MIN,  I  MAX,  .JMIN,  UMAX,  NCNT,  NAREA,  AMAX, 
LEFTX, RIGHTX, UPPERY, LOWERY, STEP 
INTEGER  ARRAY ( 100,  100 > , AMAX  v 5, 20 ) 

INTEGER  LEFTX,  RIGHTX,  UPPERY,  LOWERY 
INTEGER  YD IS,  XDIS,  AREA,  ROW,  COLUMN 
INTEGER  LOW,  DELTAY,  STEP,  STEP2 
INTEGER  POINTS (3,  500),  SUM (256) 


THIS  PROGRAM  CAN  RUN  ON  EITHER  TERMINAL  FOR  SCAN  PURPOSES 


SET  UP  BEGINNING  LOCATION-  TO  SCAN 


STEP2=2*STEP 

K  =  0 
M  =  0 

DELTAY  =  0 
LOW  =  LOWERY 
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N- J+ 1 

IF  t  M  GT  1  )  GO  TO  '30 

IF  A iM.N  :•  GE  L  TH  AND.  AtM- Ni  LE.  UTH'GO  TO  30 
IF  O  '  i'1,  H  >  GE  L  T  HF  *  GU  T U  30 
=<  1 0=N 

.JO--N 

FOUND  BEGINNING  OF  HOLE.  GO  TRACE  EDGE  ETC 
CALL  HOLE  i  I.  J,  10-  JO.  LTH.  IJTH-  PERIM,  NBIAS,  B,  NEP ) 

CHECK  TO  SEE  IF  POINTS  FOUND.  DEPENDENT  ON  PERIMETER  SIZE  CLASSIFY 
AS  A  HOLE  OR  SPOT. 

IF  (PERIM.  LE.  0.  O'  GO  TO  30 
IF  ( PERIM  LE.  20.  >NSP0T=NSP0T+1 
IF ( PERIM  GT  20.  ' NHOLE=NHOL£+ 1 
0  I F  ( A  (  I ,  J  > .  GT.  NE:  I  AS- '  GO  T 0  1 4t- 

FOUND  INTERNAL  POINT.  MARK  IT  WITH  NEIaS. 

All-  J)*NBIAS+A*  I--J) 

I A= I A+ 1 
I XP=ISX+I#NXS 
I  YP  = I S Y — J*N YS 

IF  l 1  DEBUG.  EQ.  O'CALL  SFDOT>  IxP,  I  YP  ) 

SET  FLAG  THAT  WE  ARE  NOW  INTERNAL  TO  SEGMENT. 

NFL=2 
GO  TO  147 

IF  WE  WERE  ON  EDGE  OF  SEGMENT  AND  CURRENT  POINT  HAS 
NOT  BEEN  MARKED,  THEN  WE  MUST  BE  MOVING  COMPLETELY 
OUTSIDE  SEGMENT.  SET  FLAG. 

46  IF  (At  I.  J).  LT.  NBIAS.  AND.  NFL  EQ.  1  >NFL=0 
I F  <  A  ( I ,  -J ) .  GE.  NB I  AS  •  NFL=  1 

CALCULATE  AREA  OF  SEGMENT.  METHOD  USED  INVOLVES  LOOKING  AT 
CURRENT  POINT  AND  THE  THREE  POINTS  ADJACENT  IN  THE  I - 1 ,  J-l 
DIRECTION  THESE  3  POINTS  HAVE  ALREADY  BEEN  EXAMINED  FOR 
BEING  PART  OF  THE  SEGMENT  AND  MARKED  IF  SO.  BY  COUNTING  THE 
NUMBER  OF  CORNER  POINTS  THAT  WERE  IN  THE  SEGMENT,  THEN 
THE  PORTION  OF  THE  BOX  THEY  FORM  THAT  IS  CONTAINED  IN  THE 
SEGMENT  CAN  BE  DETERMINED  IF  ALL  4  POINTS  IN  SEGMENT,  THEN 
WHOLE  BOX  CONTAINED  IF  3  POINTS,  THEN  HALF  BOX.  LESS  THAN  3 
POINTS,  THEN  THERE  IS  NO  AREA  ENCLOSED.  NOTE  THAT  THE 
VALUES  ADDED  TO  THE  AREA  COMPENSATE  FOR  THE  NON-SQUARE  SCAN. 

47  NM=0 


IF<  I  Eli.  1  i GO  TO  149 
IF  i  J.  EG.  1  >  00  TO  149 
IF  \  A I .  J ) .  OE  ME  I  AS  )  NM=NM+ 1 
I F  t  A  (  I  - 1 .  J  > .  OE.  MB  I  AS  )  NM=NM+ 1 
IFf A< 1-1, J-l )  OE.  NE I AS ) NM=NM+ 1 
IF(  A!  I,  J-l  >.  OE.  NBIAS)NM=NM+1 
IF  <  MM.  EG.  3  > RAREA=RAREA+0.  625 
I F  <  MM.  EG.  4  )  RAREA-RAREA+ 1 .  25 

149  CONTINUE 

1 50  CONTINUE 


RETURN 

END 
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THIS  PROGRAM  PRODUCES  A  CARTOON  OF  THE  IMAGE  SCANNED  USING 
AN  UPPER  AND  LOWER  THRESHOLD  VALUE  ENTERED  THROUGH 
THE  CONSOLE  BY  SELECTING  THE  CORRECT  THRESHOLDS  ONE  CAN 
OBTAIN  A  RESONABLE  PICTURE  OF  WHAT  WAS  SCANNED. 

THE  PROGRAM  ALLOWS  MULTIPLE  PLOTS  USING  DIFFERENT  THRESHOLDS 
THUS  ONE  CAN  FILL  IN  AREAS  BY  DETERMINING  WHICH  THRESHOLDS 
WILL  GENERATE  WHICH  CARTuON.  TO  EXIT  THE  PROGRAM  THRESHOLD 
VALUES  OF  0  FOR  BOTH  UPPER  AND  LOWER  THRESHOLDS  WILL  CAUSE  THE 
PROGRAM  TO  GO  TO  STOP.  CHECKS  ARE  MADE  DURING  PROCESSING  TO 
MAKE  SURE  THE  FILE  IS  VALID  AND  THE  THRESHOLDS  ARE  CORRECT. 

THE  NUMBER  OF  POINTS  WHICH  WILL  BE  PLOTTED  IS  GIVEN  ON  THE 
CONSOLE  AFTER  THE  THRESHOLDS  ARE  TYPED  IN.  THIS 
TOTAL  IS  THEN  PROGRAMAT I CALL Y  CALCULATED  AND  LATER  CHECKED 
AGAINST  THE  FIRST  NUMBER  TO  VALIDATE  PROCESSING. 

SUBROUTINE  LOGIC  IS  USED  TO  UNPACK  DATA  FROM  BLOCKS  2-21. 
SUBROUTINE  SPDOT  IS  USED  TO  PLOT  ".  "  ON  THE  SCREEN  OF  THE 
TEKTRONIX  4006-1  TERMINAL. 


ARRAY  WILL  CONTAIN  THE  SCANNED  VALUES  OF  EACH  BLOCK  WHICH 
WILL  EE  UNPACKED  BY  SUBROUTINE  LOGIC. 

PACK  WILL  CONTAIN  THE  RELATIVE  BLOCK  1  THRU  RELATIVE 
BLOCK  20  DATA  WHICH  WILL  BE  READ  FROM  THE  FILENAME. 

I HOLD  WILL  BE  USED  IN  THE  UNPACKING  DO  LOOP 

SUM  WILL  HOLD  THE  HISTOGRAM  DATA  AND  IS  USED  TO  VALIDATE 

THE  FILENAME 

TOT  WILL  BE  USED  TO  ACCUMULATE  TOTALS  TO  CHECK  AGAINST 
THE  VALUES  IN  SUM  TO  MAKE  SURE  EVERYTHING  IS  OK 

SIZE  IS  THE  MULTIPLE  OF  A  100  X  100  SCAN  FILE 
WHICH  WILL  BE  DISPLAYED  ON  THE  SCREEN.  THUS  A 
SIZE  OF  4  WILL  GIVE  A  400  X  400  CARTOON  ON  THE  SCREEN. 

UTH  IS  THE  UPPER  THRESHOLD  TAKEN  FROM  THE  CONSOLE 
STEP  IS  THE  WAY  THE  IMAGE  WAS  ORIGINALLY  SCANNED  BY  SCAN2 
AND  IS  NEEDED  IN  ORDER  TO  UNPACK  THE  DATA  INTO  ITS  RELATIVE 
SCANNED  POSITION 

I NAME  IS  THE  NAME  OF  THE  FILE  TO  BE  PROCESSED. 


I NTEGER  ARRAY ( 500 ) ,  PACK  <  256 ) ,  I HOLD  1 256 ) ,  SUM ( 256 ) ,  TOT ( 256 ) 
INTEGER  SIZE,  UTH,  UTHLOC,  STEP 
DIMENSION  I NAME <  6 ) 

INITIALIZE  ARRAY 

DO  11=1,  500 
ARRAY ( I J  =  0 


CONTINUE 


ZERO  OUT  SEVERAL  ARRArS  FOR  INITIALIZATION  PURPOSES 

DO  2  1  =  1.-  25c. 

PACK* I  =  0 
I  HOLD*  I  )  ~  0 
SUM ; I )  -  0 
TOT '  I  1  =  O 
CONTINUE 

I TOT  SHOULD  BE  THE  TOTAL  NUMBER  OF  POINTS  IN  SUM 
I COUNT  IS  COUNTER  INCREMENTED  EACH  TIME  A  DOT  IS  DISPLAYED 
I BLOCK  IS  CURRENT  BLOCK  BEING  PROCESSED  AND  IS  COMPARED 
AGAINST  IBLK  WHICH  IS  EQUAL  TO  PACK ( 1 >  OF  THE  BLOCK  JUST 
READ.  THIS  COMPARISON  IS  DONE  FOR  VALIDATION  PURPOSES. 

I END  IS  FLAG  SET  AFTER  BLOCK  20  IS  READ  TO  PREVENT  READING 
PAST  VALID  DATA.  I END  IS  SET  TO  1  AFTER  BLOCK  20  IS  READ. 

I TOT  =  0 
I COUNT  =  0 
I BLOCK  =  0 
I END  =  0 

THIS  PROGRAM  SHOULD  BE  RUN  FROM  THE  DGC  TERMINAL. 

CALL  OPENU,  "*TT01"  .  0,  IER.  12S) 

TYPE  "  ENTER  FILENAME  UP  TO  10  CHARACTERS  : " 

READ; 11,  100)  I  NAME  (  1  > 

FORMAT  (S10) 

WRITE  110,  101)  I NAME (1) 

FORMAT  <"  FILENAME  IS  :  ",  S10) 

CALL  FOPEN<  2,  I NAME,  512) 

TYPE  "  FILE  IS  OPEN" 

LOCATE  TO  RELATIVE  BLOCK  0  FOR  HISTOGRAM  INFO 
CALL  FSEEK  v  2,  0) 

READ  RELATIVE  BLOCK  v  INTO  SUM 
READ  BINARY; 2)  SUM 

CHECK.  TO  SEE  IF  loouu  POINTS  IN  BLOCK  0 

DO  3  1  =  1,  25c. 

I  TOT  =  I  TOT  +  SUM< I) 

CONTINUE 

IF  ( I  TOT.  NE.  1  )  GO  TO  99u 


INITIALIZE  VARIABLES  FOR  MULTIPLE  PASSES 
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O  O  O  -fc>  O  O  O  u  O  O  O  O  O  O  0)  O  O  O  ■' 


THESE  VARIABLES  WILL  EE  REINITIALIZED  FOR  EACH  2  NEW 
UPPER  AND  LOWER  THRESHOLDS  ENTERED  THROUGH  THE  CONSOLE. 

LTH  IS  LOWER  THRESHOLD 

ITOTPTS  IS' TOTAL  NUMBER  OF  POINTS  FROM  LOWER  TO  UPPER 
THRESHOLD  WHICH  SHOULD  BE  PLOTTED  BASED  ON  THE  VALUES 
FROM  THE  ARRAY  CALLED  SUM.  THIS  TOTAL  (ITOTPTS)  IS  CHECKED 
AGAINST  I COUNT  TO  VERIFY  IF  EVERYTHING  IS  GFK. 

1  UTH  =  0 

LTH  =  0 
UTHLOC  =  0 
LTHLOC  =  0 
ITOTPTS  =  0 
SIZE  ■  0 
STEP  =  2 


ACCEPT  "  ENTER  LOWER  THRESHOLD  ;  ",  LTH 
ACCEPT  "  ENTER  UPPER  THRESHOLD  :  ",  UTH 

VERIFY  THE  THRESHOLDS  ACCEPTED  FROM  CONSOLE 

IF  (UTH.  LT.  LTH)  GO  TO  991 

IF  (LTH.  EG.  0.  AND.  UTH.  EG.  0)  GO  TO  999 

IF  (LTH.  LT.  0  OR.  UTH.  GT.  255)  GO  TO  992 

CALCULATE  NUMBER  OF  POINTS  THAT  SHOULD  PRINT 
USING  BLOCK  0  MATRIX 

LTHLOC  =  LTH  +  1 
UTHLOC  =  UTH  +  1 
DO  4  M  =  LTHLOC,  UTHLOC 
ITOTPTS  «  ITOTPTS  +  SUM(M) 

CONTINUE 
TOTL= ITOTPTS 

TYPE  "  TOTAL  NO  OF  POINTS  TO  PLOT  IS  :  ",  TOT 

ACCEPT  "ENTER  STANDARD  NUMBER  OF  POINTS  FOR  CARD  BEING  TESTED" , STNRD 
CALL  TEST  < TOTL, STNRD  > 

IF  (  ITOTPTS.  EQ.  0)  GO  TO  993 

INITIALIZE  DISK  FILE  TO  READ  BLOCKS 

SIZE  =  4 
IMIDX  =  512 
IMIDY  =  370 

ISTARTX  =  IMIDX  -  ((SIZE  *  100)  /  2) 

ISTARTY  =  IMIDY  +  ((SIZE  *  100)  /  2) 

IULX  =  ISTARTX 
IULY  =  ISTARTY 
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I ELOCh  =  0 
I  END  = 

I  BLK  »  0 
I  COUNT  = 

READ  BLOCK  INTO  ARRhy  CALLED  PACK 

READ  BINARY (2)  PACK 
I BLOCK  =  I BLOCK  +  1 
I BLK  =  PACK < 1 ) 

IF  < I BLOCK.  NE.  I BLK)  TYPE  "  BLOCK  COUNTS  NOT  EQUAL  -  BUG" 
IF  (  I  BLK.  EQ.  20 )  I  END  =  1 

UNPACK  DATA 

.J  =  6 

DO  6  1  =  1  »  500 <  2 
J  =  J  +  1 

ARRAY ( I >  =  I R I OHT ( PACK <  J  > , 8 ) 

I HOLD <  J )  =  I LEFT < ARRAY (I > , S ) 

ARRAY  <  I  +  1  )  =  PACK  <  .J )  -  I HOLD < J ) 

CONTINUE 

SEARCH  ARRAY  TO  DETERMINE  WHICH  POINTS  TO  PLOT  BASED  ON 
WHETHER  THE  VALUE  LIES  BETWEEN  THE  LOWER  AND  UPPER 
THRESHOLDS. 

TWO  PASSES  WILL  BE  MADE  DOWN  EACH  COLUMN  OF  THE  SCREEN 
BECAUSE  THE  SCAN2  PROGRAM  SCANNED  IN  THIS  MANNER. 


J  =  0 
K  =  0 

DO  10  I  =  1 1  500 

IF  (ARRAY!  I  >.  GE.  LTH.  AND.  ARRAY  (I  ).  LE.  UTH  )  GO  TO  11 
J  =  J  +  1 

IF  (J.  EQ.  50)  GO  TO  13 
GO  TO  1 0 
K  =  K  +  1 

GO  TO  <14,  15,  14,  15,  14,  15,  14,  15,  14,  15),  K 

•J  =  0 

IULY  =  ISTARTY  -  SIZE 
GO  TO  10 
J  =  0 

IULY  =  ISTARTY 
IULX  =  IULX  +  SIZE 
GO  TO  10 


PRINT  DOT  FOR  THRESHOLD  VALUE  FOUND 


o  o  o 


IYPQS  =  I UL i  -  tSTEP*SIZE»*J 
IXPOS  =  I UL  X 

CALL  S  P  D  0  T  (  I  <  P  0  S  >  I < P 0 S ) 

I COUNT  =  I COUNT  +  1 
GO  TO  12 


CONTINUE 


IULY  =  I START Y 

IF  ( I  END.  EQ.  1  >  GO  TO  99 

GO  TO  5 


pf  that 


THIS  PROGRAM  PRINTS  THE  IMAGE  ON  THE  LINE  PRINTER. 

THE  THRESHOLDING  AND  OTHER  DATA  MANIPULATION  IS  SIMILAR  TO 

THAT  FOUND  IN  THE  PROGRAM  DOTMAT  ■  WHICH  IS  COMMENTED  IN  GREATER 

DETAIL.  THE  PRIMARY  DIFFERENCE  IS  THAT  THIS  PROGRAM  PRINTS 

THE  IMAGE  SIDEWAYS,  TO  PREVENT  IMAGE  REVERSAL.  OF  COURSE  A 

IS  PRINTED  INSTEAD  OF  A  DOT.  AND  THE  BLANK  SPACES  ARE  CHARACTER 

SPACES. 

THIS  PROGRAM  PRODUCES  A  CARTOON  OF  THE  IMAGE  SCANNED  USING 
AN  UPPER  AND  LOWER  THRESHOLD  VALUE  ENTERED  THROUGH 
THE  CONSOLE.  Br  SELECTING  THE  CORRECT  THRESHOLDS  ONE  CAN 
OBTAIN  A  RESONABLE  PICTURE  OF  WHAT  WAS  SCANNED. 

TO  E' IT  THE  PROGRAM  THRESHOLD  VALUES  OF  0  FOR  BOTH  UPPER 
AND  LOWER  THRESHOLDS  WILL  CAUSE  THE  PROGRAM  TO  STOP. 

CHECKS  APE  MADE  DURING  PROCESSING  TO  MAKE  SURE  THE  FILE 
IS  VALID  AND  THE  THRESHOLDS  ARE  CORRECT. 

THE  NUMBER  OF  POINTS  WHICH  WILL  BE  PLOTTED  IS  GIVEN  ON  THE 
CONSOLE  AFTER  THE  THRESHOLDS  ARE  TYPED  IN.  THIS  TOTAL 
IS  THEN  PROGR AMAT I CALL V  CALCULATED  AND  LATER  CHECKED 
AGAINST  THE  FIRST  NUMBER  TO  VALIDATE  PROCESSING. 

SUBROUTINE  LOGIC  IS  USED  TO  UNPACK  DATA  FROM  BLOCKS  2-21. 

ARRAY  WILL  CONTAIN  THE  SCANNED  VALUES  OF  EACH  BLOCK  WHICH 
WILL  EE  UNPACKED  BY  SUBROUTINE  LOGIC. 

PACK  WILL  CONTAIN  THE  RELATIVE  BLOCK  1  THRU  RELATIVE 
BLOCK  20  DATA  WHICH  WILL  BE  READ  FROM  THE  FILENAME. 

I HOLD  WILL  BE  USED  IN  THE  UNPACKING  DO  LOOP 

SUM  WILL  HOLD  THE  HISTOGRAM  DATA  AND  IS  USED  TO  VALIDATE 

THE  FILENAME. 

TOT  WILL  BE  USED  TO  ACCUMULATE  TOTALS  TO  CHECK  AGAINST 
THE  VALUES  IN  SUM  TO  MAKE  SURE  EVERYTHING  IS  OK 

UTH  IS  THE  UPPER  THRESHOLD  TAKEN  FROM  THE  CONSOLE. 

STEP  IS  THE  WAY  THE  IMAGE  WAS  ORIGINALLY  SCANNED  BY  SCAN2 
AND  IS  NEEDED  IN  ORDER  TO  UNPACK  THE  DATA  INTO  ITS  RELATIVE 
SCANNED  POSITION 

I NAME  IS  THE  NAME  OF  THE  FILE  TO  BE  PROCESSED. 


I NTEGER  ARRAY  <  500  > •  PACK ( 256 > ,  I HOLD < 256  > ,  SUM  v 256 ) ,  TOT ( 256 ) 
INTEGER  SIZE,  UTH,  UTHLOC,  STEP 
I NTEGER  PRTL I NE  U  <X> ) 

DIMENSION  I NAME <  & t 

INITIALIZE  ARRAY 


DO  1  I 


1 ,  500 


O  O  O  O  O  O  O  O  O  *•*  O  O  O  O  O  O  O  O  O  O  O  O  Ki 


ARRAY ( I )  =  0 

1  CONTINUE 

c 

C  ZERO  OUT  SEVERAL  ARRAY S  FOR  INITIALIZATION  PURPOSES 

C 

DO  2  I  =  1.  256 
PACK  (  I )  =0 
I HOLD ( I )  =  0 
SUM ( I )  =0 
TOT ( I >  =  0 

2  CONTINUE 

DO  21  I  =  1.  10O 

PRTLINE ( I )  =  "  " 

1  CONTINUE 

I TOT  SHOULD  EE  THE  TOTAL  NUMBER  OF  POINTS  IN  SUM 
I COUNT  IS  COUNTER  INCREMENTED  EACH  TIME  A  POINT  IS  PLOTTED. 
I BLOCK  IS  CURRENT  BLOCK  BEING  PROCESSED  AND  IS  COMPARED 
AGAINST  I  ELK  WHICH  IS  EQUAL  TO  PACK(  1?  OF  THE  BLOCK  -JUST 
READ.  THIS  COMPARISON  IS  DONE  FOR  VALIDATION  PURPOSES. 

I END  IS  FLAG  SET  AFTER  BLOCK  20  IS  READ  TO  PREVENT  READING 
PAST  VALID  DATA.  I END  IS  SET  TO  1  AFTER  BLOCK  20  IS  READ. 

I TOT  =  0 
I COUNT  =  0 
I BLOCK  =  0 
I END  =  0 

THIS  PROGRAM  CAN  RUN  ON  EITHER  TERMINAL 
CALL  FOPEN< 12,  "$LPT" ) 

TYPE  "  ENTER  FILENAME  UP  TO  10  CHARACTERS  : " 

READ( 1 1 ,  100)  I NAME « 1 ) 

OO  FORMAT  (S10) 

WRITE  (10,  101)  I NAME < 1 ) 

01  FORMAT  ("  FILENAME  IS  :  ",  S10) 

CALL  FOPEN( 2,  INAME,  512) 

TYPE  "  FILE  IS  OPEN" 

LOCATE  TO  RELATIVE  BLOCK  0  FOR  HISTOGRAM  INFO 
CALL  FSEEK( 2,  0) 

READ  BLOCK  0  INTO  ARRAT  CALLED  SUM 
READ  BINARY (2)  SUM 


CHECK  TO  SEE  IF  10000  POINTS  IN  BLOCK  0 


DO  3  I 


1,  256 
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r ro r  =  hot  +  sum.-  i  , 

CONTINUE 

IF  •  I TOT  NE.  1  Ov'.".’ 1  00  TO  99 0 

INITIALIZE  VARIABLES  FOR  MULTIPLE  PASSES 

THESE  VARIABLES  WILL  EE  REINITIALIZED  FOR  EACH  2  NEW 

UPPER  AND  LOWER  THRESHOLDS  ENTERED  THROUGH  THE  CONSOLE. 

LTH  IS  LOWER  THRESHOLD. 

ITOTPTS  IS  TOTAL  NUMBER  OF  POINTS  FROM  LOWER  TO  UPPER 
THRESHOLD  WHICH  SHOULD  BE  PLOTTED  BASED  ON  THE  VALUES 
FROM  THE  ARRAr  CALLED  SUM  THIS  TOTAL  C ITOTPTS)  IS  CHECKED 
AGAINST  I COUNT  TO  VERIFY  IS  EVERYTHING  IS  OK. 


1  UTH  =  0 

LTH  =  0 
UTHLOC  =  0 
LTHLOC  =  0 
ITOTPTS  =  0 
SIZE  =  0 
STEP  ®  2 


ACCEPT  "  ENTER  LOWER  THRESHOLD 
ACCEPT  -  ENTER  UPPER  THRESHOLD 


",  LTH 
" ,  UTH 


VERIFY  THE  THRESHOLDS  ACCEPTED  FROM  CONSOLE 

IF  (  UTH.  LT.  LTH)  GO  TO  99 1 

IF  (LTH.  EQ.  0.  AND.  UTH.  EQ.  0)  GO  TO  999 

IF  (LTH.  LT.  0.  OR.  UTH.  GT.  255)  GO  TO  992 


CALCULATE  NUMBER  OF  POINTS  THAT  SHOULD  PRINT 
USING  BLOCK  0  MATRIX 


LTHLOC  =  LTH  +  1 
UTHLOC  =  UTH  +  1 
DO  4  M  =  LTHLOC,  UTHLOC 
ITOTPTS  =  ITOTPTS  +  SUM(M) 

CONTINUE 
TOTL= ITOTPTS 

TYPE  "  TOTAL  NO  OF  POINTS  TO  PLOT  IS  :  ",  ITOTPTS 

ACCEPT "ENTER  STANDARD  NUMBER  OF  POINTS  TO  PLOT  IS>",STNRD 
CALL  TEST ( TOTL , STNRD ) 

IF  (  ITOTPTS.  EQ.  V)  GO  TO  993 


SET  UP  VARIABLES  TO  PRINT 
L  =  100 
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onoo'  *.  ,  o  o  cn  o  o  o 


DO  32  1  =  1.  1  «X> 

PRTLINEd;  =  "  " 

32  CONTINUE 

C 

C  WRITE  MESSAGE  TO  GIVE  FILENAME  AND  OTHER  INFO 

f 

WRITE  <12,  102  )  I NAME  < 1 ) ,  LTH,  UTH,  ITOTPTS 
102  FORMAT  <,,<14>",  "  FILENAME  IS  :  ",  S10,  "  LOWER  THRESHOLD  IS  : 

213,  "  UPPER  THRESHOLD  IS  :  ",  13,  “  TOTAL  POINTS  TO  BE  PLOTTED  IS 

315  > 

INITIALIZE  DISK  FILE  TO  READ  BLOCKS 

SIZE  =  4 
IMIDX  =  512 
IMIDY  =  370 

I START X  =  IMIDX  -  < (SIZE  *  100)  /  2) 

ISTARTY  =  IMIDY  +  ((SIZE  *  100)  /  2) 

I IJL  X  =  I  START  X 
IULY  =  ISTARTY 
I BLOCK  =  0 
I END  =  0 
I BLK  *  0 
I COUNT  =  0 

READ  BLOCK  INTO  ARRAY  CALLED  PACK 

READ  BINARY < 2)  PACK 
I BLOCK  =  I BLOCK  +  1 
I BLK  =  PACK ( 1 ) 

IF  ( I BLOCK.  NE.  I BLK)  TYPE  "  BLOCK  COUNTS  NOT  EQUAL  -  BUG" 

IF  ( I  BLK.  EQ.  20)  I  END  =  1 

UNPACK  DATA 

J  =  6 

DO  6  I  ■  1,  500,  2 
•J  =  U  +  1 

ARRAY ( I )  *  IRIGHT ( PACK ( J ) , 8 ) 

I HOLD ( J )  =  I LEFT (ARRAY ( I ) , 8) 

ARRAY  (I  +  1)  =  PACK  (  J )  -  I  HOLD  (-J) 

CONTINUE 


L  =  100 
•J  =  0 
K  =  0 

DO  10  I  =  1,  500 

IF  ( ARRAY ( I  ).  GE.  LTH.  AND.  ARRAY <  I  ).  LE.  UTH)  GO  TO  11 
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12  J  =  -J  +  1 

L  =  L  -  2 

I F  '  .J  EC1  50  '■  00  T 0  1  3 
GO  TO  1 

I  ;;  k  -  K  +  1 

00  TO  (14.  15..  14,  15,  14,  15,  14,  15,  14,  15* ,  K 

14  J  =  O 

IULt  =  I START V  -  SIZE 
L  =  99 
00  TO  10 

15  J  =  0 

IULt  =  IS TARTy 
I UL  X  =  IULX  +  SIZE 

WRITE  (12,  20*0  <  (  PRTLINE  (N  1 ,  N  =  1 ,  1 OO  > 

200  FORMAT  (  5  X  •  100A1 » 

DO  16  N  =  1,  1 00 

PRTLINE ( N  >  =  M  " 

16  CONTINUE 
L  =  100 
00  TO  1 C* 

0 

C 

II  IYPOS  =  IULY  -  <STEP*SIZE)*U 
IxPOS  =  IULX 

IXEND  *  IULX 
IYEND  =  IYPOS 
I COUNT  =  I COUNT  +  1 
PRTLINE  (  L  )  =  "#'• 

00  TO  12 
C 
C 

c 

10  CC  JTINUE 

C 

c 

c 

IULY  =  ISTARTY 

IF  ( I  END.  EQ.  1  )  00  TO  99 

GO  TO  5 


99  IF  (I COUNT.  NE.  ITOTPTS)  00  TO  994 

C 

C  LOCATE  BACK  TO  RELATIVE  BLOCK  1  FOR  MULTIPLE  PASSES 

C 

CALL  FSEEK ( 2,  1) 

00  TO  31 
C 

C  ERROR  MESSAGES 
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990 
99  1 
992 
99:3: 


994 

r 


C 

C 

999 


TYPE  "  HOT  NE  loOoO  - 
GO  TO  999 

TYPE  "  UPPER  THRESHOLD 
GO  TO  SI 

TYPE  "  THRESHOLD  VALUE: 
GO  T  0  3 1 

TYPE  "  NUMBER  OF  POINT 


BAD  DISK  " 

LT  LOWER  THRESHOLD  - ' TRY  AGAIN' 
;  OUT  OF  RANGE  0-255  TRY  AGAIN" 

>  TO  PLOT  IS  ZERO  -  TRY  AGAIN” 


GO  TO  31 

TYPE  ”  I COUNT  NOT  EQUAL  TO  ITOTPTS  -  PROGRAM  BUG" 


CALL  RESET 

STOP 

END 


! 


i 


i 


SUBROUT I NE  TEST  < PU I N Ts •  STNRD ) 

IF  < HBS <  PO I NTS-STNRD  > /STNRD-  005)1,1 ,2 

1  WRITE <12,  50) 

GO  TO  S3 

2  WRITE <12,  51) 

50  FORMA T< IX,  3HB0ARD  OK) 

51  FORMAT < IX,  2SHPR0BABLE  ERROR  ON  BOARD) 
SS  RETURN 

END 


i 


i 

j 


PROGRAM  HISTOGRAM 


THIS  PROGRAM  DISPLAYS  A  HISTOGRAM  OF  THE  SCAN  FILES 
CREATED  BY  PROGRAM  SCAN2.  THE  TEKTRONIX  4006-1 
TERMINAL  IS  USED  FOR  DISPLAY.  BOTH  TERMINALS  ARE  NEEDED  TO 
RUN  THIS  PROGRAM.  THE  PROGRAM  SHOULD  BE  RUN  ON  THE  DGC 
TERMINAL  WITH  OUTPUT  OF  THE  HISTOGRAM  GOING  TO  THE  TEKTRONIX 
TERMINAL  THE  NAME  OF  THE  FILE  WHOSE  HISTOGRAM  IS  WANTED  IS 
ENTERED.  THE  PROGRAM  ONLY  HAS  TO  ACCESS  THE 
FIRST  BLOCK  OF  THE  FILE  TO  OBTAIN  THE  HISTOGRAM  DATA.  THIS 
DATA  IS  READ  INTO  AN  INTEGER  ARRAY  CALLED  SUM  WHICH  IS 
256  WORDS  IN  LENGTH.  A  DO  LOOP  THROUGH  THIS  ARRAY  CREATES 
A  SERIES  OF  LINES  WHOSE  LENGTH  EQUALS  THE  VALUE  IN  THE 
WORD  BEING  PROCESSED.  A  SUBROUTINE  CALL  TO  SPLIN  PLOTS  A 
HORIZONTAL  LINE  ON  THE  SCREEN.  256  LINES  WILL  BE  PLOTTED. 
AFTER  THE  LINES  ARE  PLOTTED,  THE  SUBROUTINES  SPLIN  AND  SPMOV 
ARE  USED  TO  PRODUCE  SCALE  VALUES  UNDER  THE  HISTOGRAM  USEFUL 
IN  OBTAINING  THRESHOLD  VALUES  FOR  VALLEYS  AND  PEAKS. 


SUM  CONTAINS  THE  FIRST  BLOCK  OF  DATA  READ  FROM  THE  INPUT  FILE 
I NAME  IS  NK  ME  OF  FILE  TO  BE  READ 

INTEGER  SUM (256; 

DIMENSION  I NAME ( 6 ; 

INITIALIZE  SUM 

DO  11=1,  256 
SUM ( I )  =0 
CONTINUE 

TYPE  "  ENTER  FILENAME  UP  TO  10  CHARACTERS  :  " 

READ  (11,  1 00 )  I  NAME  <  1  ) 

•O  FORMAT  <S10) 

WRITE  (10,  101;  INAME(1? 

>1  FORMAT  ("  FILENAME  IS  :  ”,  S10) 

CALL  FOPEN<  2,  I NAME,  512) 

TYPE  "  FILE  IS  OPEN  " 

CALL  OPEN  (1,  "*TT01 ”,  0,  IER,  128) 

SET  UP  INITIAL  STARTING  POSITION  FOR  PLOTTING 

IXPOS  =  2 
IYPOS  =  60 
IXEND  =  2 
I  VEND  =  60- 

READ  FIRST  BLOCK  OF  DATA  FROM  FILE  OPENED  INTO  ARRAY  SUM 
READ  BINARY ( 2  >  SUM 
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2/2 


UNCLASSIFIED 


F/G  9/5 


NL 


o  o  o  o  o  o  >  i  o  o  to  ►*  o  o  o  o i  o 


CHECK  DISK  FILE  OUT  TO  SEE  IF  CREATED  BY  SCAN2 
I  TOT  =  0 

DO  1 1  J  =  1 ,  25* 

I  TOT  =  I  TOT  +  SUM  ( .J  > 

CONTINUE 

IF  I  TOT  IS  NOT  EQUAL  TO  10000  THE  DISK  FILE  IS  NOT  GOOD 
IF  (I TOT.  EQ.  1 000*  >)  GO  TO  13 

TYPE  "  FILENAME  WAS  NOT  CREATED  BY  SCAN2  JOB  ENDS" 

GO  TO  99 


PLOT  25*  HORIZONTAL  LINES 

DO  2  I  =  1,  25* 

IYEND  =  IYEND  +  SUM  v I ) 

CALL  SPL IN  (IX POS •  I YPOS .  I X END ,  I YEND ) 

IXEND  =  IXEND  -*■  4 

IxPOS  *  IXEND 

IYEND  =»  *0 

CONTINUE 

SET  UP  TO  PRINT  NUMBERS  UNDER  HISTOGRAM 


GET  BACK  TO  STARTING  POINT 

IX POS  =  2 
I YPOS  =  60 
IXEND  =  2 
IYEND  =  60 

CALL  SPLIN  (IXPOS,  IYF'OS,  IXEND,  IYEND) 

DO  THREE  LONG  LINES 

IYEND  =  40 
DO  20  J  =  1,  3 

CALL  SPL IN ( IXPOS,  IYPOS,  IXEND,  IYEND) 
IXPOS  =  IXPOS  +  400 
IXEND  =  IXPOS 
CONTINUE 

RESET  VARIABLES  AND  DO  THREE  MIDDLE  LINES 

IYEND  =  50 
IXPOS  =  202 


IX END.  I VEND) 


1  <  END  -  202 
DO  '30  J  =  1 .  3 

CALL  3 PL  I N  UX PC'S.  I  f  POS. 

IxPOS  =  IxPOS  +  4oo 
lx END  =  IxPOS 
CONTINUE 

SET  UP  VARIABLES  FOR  2c.  SMALL  LINES 

I  VEND  =  55 
I  r  POS  =  c-0 
IXPOS  =  2 
IX.END  =  2 
DO  40  .J  =  1  ■  26 

CALL  SPLIN  (IXPOS,  IV POS,  IXEND,  I YEND  > 
I x  POS  =  IXPOS  +  40 
IX.END  =  IxPOS 
CONTINUE 


GET  BACK  TO  PRINT  NUMBERS  UNDER  HISTOGRAM 


IXPOS  =  2 
IYPOS  =  20 
IXEND  =  2 
I VEND  =  20 

CALL  SPMOV  (IxPOS,  IYPOS/ 

WRITE  (1)  "0" 

PRINT  NUMBER  50 

IXPOS  =  192 
I YPOS  =  20 
IXEND  =  192 
I YEND  =  20 

CALL  SPMOV (IXPOS,  IYPOS) 

WRITE(l)  "50" 

DC*  LOOP  TO  GIVE  l1.1*.',  150,  200,  250 

I  XF'OS  —  385 
I YPOS  =  20 
IXEND  =  385 
I YEND  =  20 


DO  50  1=1,4 

CALL  SPMOV (IXPOS,  IYPOS) 

GO  TO  <51,  52,  53.  54),  I 


B-38 


1 


SUBROUTINE  TRACE(  IRON.  .JCOL,  10.  JO,  LTH,  IJTH,  F'ERIM,  AREA,  NSPL. 

IMIN.  I  MAX.  .JM  IN,  JMAX,  B,  NF'T) 

PREPARED  FOR  SCI  SY STEMS •  INC. 

BY  BILL  POPE,  TELCOM  DATA  CORPORATION 

ONCE  AN  AREA,  OR  HOLE  WITHIN  AN  AREA,  HAS  BEEN  LOCATED 
THIS  SUBROUTINE  TRACES  THE  EDGES  AND  MARKS  THEM  AS  PART 
OF  THE  AREA  ADDITIONALLY,  TRACE  BUILDS  A  TWO  DIMENSIONAL  ARRAY,  B, 
WHICH  CONTAINS  THE  ORDERED  I  AND  J  VALUES  OF  THE  EDGE  POINTS  OF 
THE  AREA  OF  HOLE  BEING  TRACED. 


SOME  OF  THE  VARIABLES  USED  ARE: 


IRON. JCOL 

IO,  JO 

IP.  JP 
LTH, UTH 

AREA 

NSPL 


*•* 


STARTING  COORDINATES,  WITHIN  ARRAY.  OF  THE  EDGE  TO 
BE  TRACED. 

NEXT  TO  LAST  POINT  FOUND  ON  EDGE  BEFORE  CURRENT  POINT. 
LAST  POINT  FOUND,  USED  AS  CENTER  OF  SEARCH  FOR  NEXT  Pt 
LOWER  AND  UPPER  THRESHOLDS  PASSED  IN  FROM  BOARD.  USED 
AS  DELIMITERS  FOR  THE  FEATURES  IN  THE  SCAN. 

NUMBER  USED  TO  MARK  POINTS  AS  BELONGING  TO  A  GROUP. 
NUMBER  USED  TO  MARK  POINTS,  WITHIN  AN  AREA,  AS 
BEING  PART  OF  A  SUBGROUP (  EDGE  OF  HOLE). 

BOTH  AREA  AND  NSPL  ARE  POWERS  OF  2  SO  THEY  BASICLLY 
JUST  SET  A  BIT. 


IMIN, IMAX 
JMIN, JMAx 


THE  MIN  AND  MAX  COLUMNS  OF  THE  EDGE. 
THE  MIN  AND  MAX  ROWS  OF  THE  EDGE. 


PERIM 


MEASURED  PERIMETER  OF  THE  AREA  TRACED  AROUND. 


B 

NF'T 


AN  ARRAY  CONTAINING  THE  I,J  VALUES  OF  THE  EDGE. 
NUMBER  OF  POINTS  IN  B  WHICH  ARE  THE  CURRENT  EDGE. 


COMMON  ARRAY. SUM, LP,  IMN,  IMX, JMN, JMX, NCNT , NAREA, AMAX (5, 20) 
COMMON  /PLT/ISTARTX, I ST ARTY, MIDX, MIDY, NxSIZE, NYSIZE 
INTEGER  ARRAY  (  100,  1  VO  >,  SUM (256.) 

INTEGER  AREA,  UTH, UTHF 
INTEGER  B (2, 1500) 

VERIFY  STARTING  ROW  AND  COLUMN  ARE  VALID. 


IF  (  I  ROW.  GT 
IF-  I  ROW.  LT. 
IF  (.JCOL  GT 
IF  (JCOL.  LT 
I = I ROW 
•J=JCOL 


1 00 ) GO  TO  905 
1 >G0  TO  905 
1 00 ) GO  TO  905 
1 ) GO  TO  905 
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SET  STARTING  POINT  AS  FIRST  IN  EDGE 
NPT  =  1 

E(l- NPT ) = I ROW 
E  (  2,  NPT  =-JCOL 

C  INITIALIZE  PERIMETER  AND  MIN  MAX 


PERIM=0.  0 
IMAX= IRON 
I M I N= I ROW 
JMAX=JCOL 
•JM I  N= JCOL 
I  XE*0 


ONCE  A  POINT  HAS  BEEN  MARKED  BY  ADDING  AREA  TO  IT,  THEN 
YOU  NEED  A  DIFFERENT  THRESHOLD  VALUE.  SET  LOWER  AND  UPPER 
THRESHOLDS  SO  WE  CAN  STILL  IDENTIFY  POINTS  FROM  THE  CURRENT 
EDGE  AFTER  THEY  HAVE  BEEN  " MARKED 

LTHF=LTH+AREA 

UTHF=UTH+AREA+NSPL 

DI.  D-J,  DIAG  ARE  THE  VALUES  USED  IN  CALCULATING  THE  PERIMETER 
THEY  VARY  DEPENDENT  ON  DIRECTION  OF  MOVEMENT  FROM  ONE  POINT 
TO  THE  NEXT.  THE  DIFFERENCE  IN  DI,  DJ  IS  DUE  TO  THE  CAMERA 
NOT  HAVING  A  SQUARE  SCAN. 

D I = 1  25 
D-J=  1  0 

DIAG=1  60078 

INITIALIZE  VALUES  FOR  TEKTRONIX  PLOTTING. 

I XP*ISTARTX+I*NXSIZE 

I  YP=  I  STARTY- -J*NYS  I Z  E 

MS=2 

GO  TO  21 

CONTINUE 

HAVE  WE  GONE  FULL  CIRCLE 

IF  (I.  EQ.  I  ROW.  AND.  J.  EQ.  JCOL)  GO  TO  77 


PUT  AREA  FLAG  ON  DATA  POINT 


ARRAY <  I ,  J ) =AREA+NSPL+ ARRAY (I , J ) 


If  HEW  POINT 


MIN  OP  Mh'  VALUE 


IP  ■  I  L  T.  ININ  .  IMIN=  I 
IF*  I  Of  IMAx  )  I M1-1X  =  I 
I F  (  J.  LT  .JMIN).JMIM=J 
I F  >  0.  OT  JMh  x  )  .JMh  <  =.j 

IF  MARKING  NEU  EDGE.  PLOT  IT 


IF  (AREA  EG.  0)G0  TO  1<.> 

NS  IS  A  VALUE.  1  TO  8.  WHICH  IS  THE  DIRECTION  THE  EDGE  IS  MOVING. 
AS  LONG  AS  POINTS  CONTINUE  IN  SAME  DIRECTION,  SKIP  PLOT  AND  THEN 
DRAW  ONE  LONG  VECTOR  INSTEAD 

IF  (MS.  EQ  NS)  GO  TO  1 «.«  4ft 

CALL  SPLINvIxR.  IVP-  I xE,  I yE > 

IxP=IXE 

IYF-IYE 

MS=NS 

CONTINUE 

ADD  CORRECT  VALUE  TO  PERIMETER. 

IF (  I.  EQ.  IP) GO  TO  1 1 

IF(.J.  EQ.  JP ) GO  TO  12 

PER I M=PER I M+D I AG 

GO  TO  13 

PER I M=PER I M+OJ 

GO  TO  13 

PEP I M=PER I M+D I 

CONTINUE 

PLACE  NEW  POINT  IN  ARRAy  B 
NPT=NPT+1 

IFfNF'T.  GT.  1500)  GO  TO  900 
B< 1, NPT ) = I 
B  (  2 ,  NPT)=.J 
CONTINUE 

SAVE  CURRENT  POINT  FOR  PLOTTING. 

I XE= ISTARTX+ I *NxS I  i  E 

iye=istapty-.j*nysize 

I  p=  I 
JP=.J 


FOUND  AND  PLOTTED  A  POINT  SEARCH  FOR  NEXT  POINT  BEGINS. 
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C 

c 

c 

c 

c 

c 

c 

115 

116 

117 

ns 

120 

125 

126 

127 

129 


ASSUME  THE  SEARCH  OF  ADJACENT  POINTS- THE  9  POINTS  TO 
BE  EXAMINED  ARE  NUMBERED  1  TO  9  STARTING  WITH  THE  ONE 
AT  THE  I-l-J  POSITION  AS  1  AND  MOVING 
CLOCKWISE.  FOR  THE  POINT  WE  ARE  CURRENTLY  AT,  THERE 
WILL  BE  A  PREVIOUS  POINT  AS  ONE  OF  THE  9  OUR 
EXAMINATION  WILL  START  WITH  THE  NEXT  CLOCKWISE  POINT 
AND  PROCEED. 

IF  (  J— JO )  125,  120,  115 

IF(I-IO)  118,  117,  11/;- 

NS =3 

GO  TO  19 
NS»4 

GO  TO  19 
NS=5 

GO  TO  18 
NS=2 

IF (  I  LT.  1 0 ) NS=6 
GO  TO  18 

IF(I-IO)  126,  127,  129 

NS=7 

GO  TO  19 
NS=9 

GO  TO  18 
NS=1 


THROUGH  WITH  OLD  POINTS.  SAVE  CURRENT  POINTS  FOR  USE  NEXT  TIME 

8  CONTINUE 

10=1 
,J0=J 

GO  TO  (  25,  30,  35,  40,  45,  50.  55,  60 ) ,  NS 

NOW  WE  KNOW  WHERE  TO  START,  CIRCLE  CLOCKWISE  LOOKING  FOR 
NEXT  POINT  IN  AREA  BOUNDARY. 

5  CONTINUE 

TO  PREVENT  LOOPING  FOR  SPECIAL  CASE  OF  ISOLATED  POINT 
ON  ROW  1  (1=1),  NEED  FOLLOWING  CHECK. 

IF  (IP  EG.  1.  AND.  NS.  EQ.  2)  GO  TO  65 
C 

I=IP-1 

IF <  I  LT  1  ) GO  TO  34 
J=JP 

I V=ARRAY ( I. J) 

IF(  IV.  GE.  LTH.  AND.  IV.  LE  UTH)GO  TO  20 
IF (  IV.  GE.  LTHF  AND.  IV  LE.  UTHF ) GO  TO  76 
IF  (NS.  EQ.  2)  GO  TO  65 
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V-v"  "w 


■-  '  ’  T' 


*11  '  '  V  K~". 


cunt i nue 
i=ip-i 

I F  <  I  LT  1  >  GO  TO  35 
■J= JP- 1 

I F  * J.  LT.  1  )  GO  TO  45 
I V= ARRAY ( I, J) 

IF*  IV.  GE  LTH  AND.  IV  LE  UTH)GO  TO  20 
IF(  IV  GE.  LTHF.  AND.  IV.  LE.  UTHF )  GO  TO  76 
I F  (  NS.  EQ.  3 )  GO  TO  65 


CONTINUE 

DITTO  COMMENT  AT  25 

IF  *  IP.  EQ  1.  AND.  NS.  EQ.  2)00  TO  65 

CONTINUE 
J=.JF- 1 

IF'  J.  LT.  1  ) GO  TO  45 
I  =  IP 

I V=ARRAV ( I , J ) 

IF (IV.  GE.  LTH.  AND.  IV.  LE.  UTH)GO  TO  20 
IF<  IV.  GE.  LTHF.  AND.  IV.  LE.  UTHF) GO  TO  76 
IF  (NS.  EQ.  4)  GO  TO  65 


I*IP+1 

IF*  I.  GT.  100) GO  TO  55 
■J=JP-1 

IF*  J  LT.  1  )  GO  TO  45 
I  V= ARRAY  *  I ,  .J) 

IF*  IV.  GE.  LTH.  AND.  IV.  LE.  UTH)GO  TO  20 
IF*  IV.  GE.  LTHF.  AND.  IV  LE.  UTHF ) GO  TO  76 
IF  (NS.  EQ.  5)  GO  TO  65 


I=IP+1 

IF >'  I.  GT.  1 00 ) GO  TO  55 
U=JP 

I V=ARRA V  (  I,  .J) 

IF*  IV  GE  LTH  AND.  IV.  LE  UTH  i  GO  TO  20 
IF* IV  GE  LTHF  AND  IV.  LE  UTHF ) GO  TO  76 
I F  *  N  .  E  Q .  U  Ci  X  Ci  5 
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1= IP-*- 1 

I F ( I  GT  1 00 ) GO  TO  55 
J=.JP+ 1 

IF(.J  GT.  100) GO  TO  25 
IV=ARRAY ( I  - J) 

IF(IV,  GE.  LTH.  AND.  IV.  LE.  UTH)GO  TO  20 
IF  (  IV.  GE.  LTHF.  AND.  IV.  LE.  IJTHF ) GO  TO  76 
IF  (NS.  EG!.  7  )  GO  TO  65 


J=.JP+1 

IF ( J.  GT.  100) GO  TO  25 

I  =  IP 

I V=ARRAY  (  I,  .J) 

IF(  IV.  GE.  LTH.  AND.  IV.  LE.  UTH)GO  TO  20 
IF (IV.  GE.  LTHF.  AND.  IV.  LE.  UTHF ) GO  TO  76 
IF  (NS.  EQ.  3)  GO  TO  65 


I  =  IP-1 

IF  (I.  LT.  1  )G0  TO  34 
J=JP+1 

I F  ( J.  GT.  100)  GO  TO  25 
IV=ARRAY(  I,  .J) 

IF ( IV  GE.  LTH.  AND.  IV.  LE.  UTH>GO  TO  20 
IF (IV.  GE.  LTHF.  AND.  IV.  LE.  UTHF) GO  TO  76 
IF  (NS.  EQ.  1  )  GO  TO  65 
GO  TO  25 


COMPLETED  CIRCLE  WITH  NO  FIND.  SEE  IF  WE  ARE  STILL  ON 
FIRST  POINT  IDENTIFIED. 

CONTINUE 

IF (  IP.  EQ.  IROW.  AND.  JP.  EQ.  JCODGO  TO  66 
TYPE  " BOY  ARE  WE  LOST ! M 
GO  TO  78 


STILL  ON  FIRST  POINT  SO  CONSIDER  IT  AS  ISOLATED  BY 

SETTING  VALUE  TO  ZERO 

ARRAY ( IROW,  JCOL ) =0 

PER  I  M=0.  0 

GO  TO  78 


175 

176 
1^7 


C 

c 

78 

900 


FND  OF  LOOP  FLAG  STARTING  POINT  AND  CONTINUE 
CONTINUE 

IF,  I  NE  IRON  OP  J  NE  JCOL  >  GO  TO  10 

ARRAY  t  I  ROW.  JCOL  >  =AREA+NSPL+ARRAY  <  IRON-  JCUL  > 

NPT=NPT+1 

IF  (  NPT.  GT,  1500)  GO  TO  900 

Ev 1. NPT ) = I 

B  v  2. NPT > =J 

IF  (I  EO.  IP)  GO  TO  175 

I F  <  J.  EG.  UP  * GO  TO  17c- 

PER  I  fi=PER  I N+D I  AG 

GO  TO  177 

PER I M=PER I M+D  J 

GO  TO  177 

PEP  I  ri=PER  I  M+D  I 

CONTINUE 

I XE= I STARTX+ I R0N*NXSI ZE 
I YE=ISTARTY— JCOL*NYSIZE 

IFtAREA.  NE.  0)CALL  SPLlNtlxP.  IYP.  IXE-  I  YE ) 


CONTINUE 

RETURN 

TYPE  "  TRACE  FOUND 
STOP 

TYPE  "  BAD  TRACE-" 

STOP 

END 


1 5u0  POINTS  IN  ONE  AREA. 


•905 


IRON,  JCOL 


o  o  o  o  n)  o  o  o  o  o  o  o  o  o  o  o  onor 


SUBROUTINE  TRACK  < IROW, .JCOL,  IQ, JO, LT1, UT1,  ZMIN, MTOUCH) 

PREPARED  FOR  SCI  SYSTEMS,  INC. 

C  BY  EILL  POPE,  TELCOM  DATA  CORPORATION. 

C 

C  THIS  SUBROUTINES  DRAWS  A  BORDER  OF  WIDTH  ZMIN  AROUND  A  SEGMENT 

C  BOUNDARY  IF  ANY  OTHER  SEGMENTS  ARE  ENCOUNTERED  WITHIN  THE  BORDER 

C  THEN  THE  FLAG,  MTOUCH,  IS  RETURNED  INDICATING  CRITICAL  CLOSENESS 

C  TRACK  USES  THE  SAME  ALGORITHM  AS  TRACE  FOR  FINDING  THE  EDGE  POINT 

OF  A  SEGMENT  WHOSE  POINTS  ARE  EETWEEN  THE  VALUES  LT1 ,  UT1. 


COMMON  ARRAY, SUM, LP, IMN, I MX, JMN , JMX, NCNT , NAREA, AMAX 
COMMON  /PLT/ISTARTX,  ISTARTY,  MIDX,  MIDY, NXSIZE, NYSIZE 
I NTEGER  ARRAY ( 1 00 , 1 0O > ,  SUM ( 256 ) 

INTEGER  AREA,  UTH, UT 1 , UTHF 
INTEGER  AMAX (5, 20) 

I = IROW 

J=JCOL 

JP=J0 

IP=IO 

MTQIJCH=0 

NPT=1 

ZD=SC!RT  (  ZMIN-»*2y  2  ) 

NORMALLY,  EACH  SIDE  OF  THE  45  DEGREE  TRIANGLE  WOULD  BE  EQUAL. 
WITH  THE  CAMERA  BIAS,  WE  HAVE  DIFFERENT  SIZE  DELTAS. 

JD=IF I  X  <  ZD/ 1 .  +  5) 

ID=IFIX<  ZD/1.  25+.  5) 


SET  DELTAS  FOR  THE  PERPENDICULARS  IN  THE  BOX 

JE=IFIX<  ZMIN/1  +  5) 

IE=IFIX<  ZMIN/1.  25+.  5) 


LTH=LT 1-256 

UTH=LTH+NAREA 

LTHF=LTH 

uthf*uth 

GO  TO  21 
0  CONTINUE 


IF  ON  EDGE  OF  SCAN  FORGET  DISTANCE  CHECK 
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I  r  ' 

iT'U 

EO 

1 ) GO  70 

IF 

EO 

i  ’  00  T  0 

l  V  ■ 

I  Ml 

EO 

1  *u  ■  >  oo 

i  c 

,jM  j 

E-  0 

1  *.»v  1 00 

m':  WE  i  f'M1. E  EjGE  OF  SEGMENT •  WE  WANT  TO  PROJECT  OUTWARDS  FROM  THE 
CURRENT  POINT  SINCE  WE  ARE  MOVING  IN  THE  SAME  DIRECTION  AROUND  EDGE. 
THIS  PAN  BE  THOUGHT  OF  A  PROJECTING  TO  THE  LEFT.  FOR  THIS  PROGRAM,  THE 
DIRECTION  TO  PROuECT  HAS  FROM  THE  CURRENT  POINT  1,0  HAS  BEEN  DETERMINED 
E:  THE  ANGLE  RELATION  BETWEEN  IT  AND  THE  TWO  PREVIOUS  POINTS  I-l.J-l  AND 
1-2. J- 2  CALL  THEM  POINTS  P.  PI.  AND  P2.  IF  THE  POINTS  FORMED  A 
1  Bn  DEGREE  ANGLE  THE  A  SINGLE  PROJECTION  AT  90  DEGREE  FROM  P  WAS  ALL 
ThhT  WhS  PEGU I  RED  SIMILARL  Y  FOR  ANY  ANGLE  LESS  THAN  90  DEGREES. 

FOR  ANGLE'S  GREATER  THAN  9>.>  DEGREES,  YOU  WOULD  BE  TURNING  AN  OUTSIDE 
CORNER  AND  THUS  WOULD  NEED  ADDITIONAL  PROJECTIONS  TO  PROVIDE  GOOD 
COVERAGE  OF  POINT'S  AROUND  THE  OUTSIDE  CORNER. 

POP  INSTANCE.  WITH  THE  OUTSIDE  ANGLE  OF  P, P1.P2  BEING  225  DEGREES, 
y  0 U  W 0  U L  D  HAVE  S  P F' 0  J EC  T I ONS : 

90  DEGREEE'S  TO  P-Pl  STARTING  AT  P. 

DEGREES  TO  F-P1  STARTING  AT  PI. 

’  1 ?5  DEGREES  TO  P-P1-P2  STARTING  AT  PI *  SPLIT  THE  OUTSIDE  ANGLE) 

By  USING  THIS  TECHNIQUE.  WITH  EVEN  MORE  PROJECTIONS  FOR  LARGER  ANGLES, 
yOU  ARE  mBLE  TO  FORM  A  FAIRLY  COHERENT  BORDER  WITH  A  MINIMUM  OF 
P  RO  JEC  T  I  ON'S 


FOUND  NEW  POINT  ALONG  EDGE  FIND  WHICH  DIRECTION* S )  TO  EXTEND 
MINIMUM  VALUE  AND  SEE  hBGUT  CRITICAL  DISTANCES. 


FOP:  ANGLE'S  MOVING  DIAGONAL  TO  THE  SCAN,  THE  DISTANCE  TO  PROJECT  WILL 
BE  SOME  COMBINATION  OF  THE  "DELTA'S"  ID  AND  JD.  FOR  PROJECTIONS 
MOVING  HOP  I ZONALL i  THE  DELTA  IS  +0R-  IE.  FOR  PROJECTIONS 
MOVING  VERTICALLY  THE  DELTA  IS  +0R-  IE. 

WHEN  A  DIRECTION  OF  PROJECTION  FROM  A  POINT  IS  DECIDED,  CRF'T  IS 
CALLED  TO  FIND  IF  THE  PROJECTION  FINDS  A  CRITICAL  DISTANCE. 

THE  POINTS  INVOLVED  ARE  I.J  10, JO 
I F ' J- JO i 5- 10  is 
I F  a  -  I  0  >  .  7,s 
1 1  =  -  I D 
J1  =  JD 

C  mLL  O  F  T  •  I  .  J.  I  1  •  Jl  .  MTOUC'H.  L  I  H.  UTH  > 

I57'  I  M2  OT  10 'GO  r0  l':' 

IP’  I  EO  I  M2  AND  J  NE  UM2  ’GO  TO  2'r-' 

CALL  C  1  FT  ’  10.  JO-  I  1  ■  J 1  ■  MTOUC H,  L TH .  LIT H  > 


I  M2,  JM2 


CALL  edge  to  BREAK  ARRAY  into  areas  within  lth,  uth 


CALL  AEDGE'LTL, LTU) 


CAlL  OVLOD (5, OVERS. 0. I ERR) 
IF( I ERR  NE  1 > GO  TO  9So 


NOW  TRY  TO  FIND  DISTANCES  BETWEEN  AREAS 


CALL  DISTANCE 

IF  (NANS.  EG.  "N "  )  GO  TO  40 


SET  VALUES  FOR  AUTO  MOVE  TO  NEXT  SCAN 

LEFTX=LEFTX+NXSI Z 

IF ( LEFTX.  LT.  NXSTOFuGO  TO  10 

MOVE  UP  ON  SCREEN  AND  MAKE  NEXT  PASS 

LEFTX=NXST 

LOWERy=LOWERY+NYSIZ 

IF (LOWERY.  LT  NYSTOF ) GO  TO  lv 


END  OF  SCREEN 
CONT I NUE 

CALL  CLOSE (5, I ERR) 

CALL  RESET 
STOP 

TYPE  "  ERROR  LOADING  OVERLAY  ",  I ERR 

STOP 

END 


B-t.2 


continue 

INITIALIZE  DISK  FILE  AND  UNPhDK  EiATA 

CALL  OVLOD  v  OVER1.'.  '-•>  I  ERR  t 
IF (I ERR  ME  1 ' GO  TO  *30 

CALL  UNPACK 

CALL  OVLOD '  5>  0VER1  ■  I  ERR  - 

IF'  I  ERR  NE.  1  )  CO  TO  9S>.' 

SMOOTH  DATA  IF  REQUIRED 

CALL  SMOOTH 

CONTINUE 

CLEAR  TEKTRONIX  SCREEN 
WRITE! 1 • 245) 

FORMAT ( "  <33>< 1 4>  " > 

DETERMINE  THRESHOLD 

CALL  OVLOD! 5. 0VER6.  u,  IERR) 

IF ! IERR  NE.  1)00  TO  980 

LTU=0 

CALL  CONVAL ( ARRAY . LP>  LTL>  LTU ) 


FIND  L'U '  ABU'UT  LIVE  ‘E'Lr-tN 


rvpE  "  LIVE  SCAN0  y,u" 
READ  11  1  >  NAN';. 

IF'NhNS.  EC'  "N"  >0Q  TO  2'.* 
I F < NANO  NE  " Y "  ■  GO  TO  1 2 


SET  OF'  EOF  LIVE  SCANS 


FORMAT'.  SI  > 

INITIALIZE  SCAN  VALUES  TO  START  WITH  LOWER  LEFT  OF  VISIBLE  SCREEN. 

NX ST  ® 1 5v 
NYST= 1 50 
LEFTX=NXST 
LOWERY=NYST 

DETERMINE  FOR  A  GIVEN  STEP  SIZE  HOW  MANY  SCANS  ACROSS  AND  DOWN 
THE  SCREEN.  SET  UP  VALUES  AND  ENDING  CONDITIONS. 

MS=400-2*LEFTX 

STEP=1 

NUM=MS  >  (, STEP* 1 OU >  + 1 
NXSIZ=MS/NUM 

NXST0F'=4C|C'-LEFT X-NXSI  Z/2 
MS=500-2*L0WERY 
NUM=MS./  (  1 00* STEP  .»  +  1 
NYS I Z  =MS / NUM 

N  Y  S  T  0  P = 5  0  O  -  L  0  W  E  R  Y  -  3  *  N  Y  S I  Z  /  2 
NUMSCN=0 

R I GHT X  =99*STEP+LEF T  X 
UPPER  Y =99*STEF’+L0WER  Y 
N  U  M  S  C  N = N I J  M  S  C  N  +  1 
WR I TE ( LP, S3 7 ) NUMSCN 

FORMATt TSO. "*  *  *  SCAN  NUMBER  "• IS- "  *  *  *" > 


GO  SCAN  SCENE 

CALL  TSPOT \ LEFTx, LOWERY, ARRAY' 


LL '  ^'MhF.D 


F'PEP'ARED  FuF'  SCI  SySTEMS,  INC. 

B<  BILL  POPE,  TELCOM  DATA  CORPORATION. 

THIS  PROGRAM  IS  THE  DRIVER  ( MA I N )  ROUTINE  FOP  THE  LIVE  BOARD 
SCANNING  SYSTEM.  IT  PERFORMS  A  NUMBER  OF  SCANS,  EACH  OF  100X100  POINTS, 
TO  COVER  THE  AREA  SEEN  By  THE  CAMERA.  THEN  SUBROUTINES  ARE 
CALLED  TO  DETERMINE  THE  FEATURES  IN  THE  SCAN. 

ALTERNATELY  BOARD  CAN  EXAMINE  A  SINGLE  SCAN  PREVIOUSLY  STORED 
ON  DISC  BY  THE  PROGRAM  SCAN2. 

NOTE  SINCE  SUM  DATA  ON  DISC  WAS  STORED  BEFORE  THE  HARDWARE  FILTER 
WAS  PLACED  ON  THE  DIGITIZER,  IN  THIS  MODE  THERE  IS  AN  OPTIONAL 
DATA  SMOOTHING  SUBROUTINE. 


S80 

C 

c 

c 

c 

c 

c 


COMMON  ARRAy . SUM. LP,  IMIN,  I MAX,  JMIN,  UMAX,  NCNT • NAREA, AMAA, 
LEFT  X . RIGHT  X , UPPER y. LOWERY, STEP 
INTEGER  RIGHTx, UPPERY,  STEP 

INTEGER  ARRAY(  lO1.1,  I'.".'),  SUM  (236),  AMAX  <  5-  20  ) 

EXTERNAL  OVERv, OVERl, OVERS, OVERS 


INITIALIZE  OVERLAYS 

CALL  OVOPNf  3,  “BOARD.  OL".  I ERR) 

I F  (  I  ERR.  EG.  1  )  GO  TO  SSu 

TYPE  •*  OVERLAY  OPENING  ERROR  ",  I  ERR 

STOP 

CONT I NUE 


OPEN  TEKTRONIX  TERMINAL  FOR  PLOTTING 

CALL  0PEN(1,  "tTTOl  "  . IER,  128) 

LTU=0 

LP=10 

TYPE  "  PRINT  OR  TYPE  RESULTS?  P, T  " 
READ ( 1 1, 100 ) NANS 
IF  (NANS.  NE  "P“  )  GO  TO  7 
LP=  1 2 

CALL  FOF'EN  >  12,  "$LPT"  ) 

CONTINUE 


i 

i 


I 

■x 


Es-59 


READ  BLOCK  INTO  ARRAY  CALLED  PACK 


Ml  =  1 
1  =  1 

IBLOCK-U 
I END-0 
I ELK-0 
I COUNT =0 

READ  BINARY '2.'  PACK 
I BLOCK- I BLOCK* 1 
I ELK  =  PACK ( 1 ) 

IF  ^ I BLOCK.  NE  I ELK) TYPE  "  BLOCK  COUNTS  NOT  EQUAL  -  BUG 
IF ( I ELK.  EQ  20)  IEND=1 

UNPACK  DATA 


DO  6  J=7. 256 

ARRAY < MC ,  I )  =  IRIGHT ( PACK  <  J ) . 8 > 

I  HOLD- 1 LEFT ( ARRAY ( MC.  1)8) 

ARRA y ( MC ■  I +2 ) -PACK  UJ ) - 1 HOLD 
1  =  1+4 

IF l I  LT  lOO >00  TO  6 
1=1-99 

IF  (  I.  EQ.  2)00  TO  6 

1  =  1 

MC-MC+1 

CONTINUE 

IFlIEND.  EQ  0)00  TO  5 
RETURN 

TYPE  "  NOT  10.  OoO  POINT  SCAN  FILE” 

STOP 

END 


E 


'  *7V1V  w  T  V  V  ^  ’  -•  »  V—  J-J-7W  -  w  - 


DO  1  J= 1 .  1 00 
DO  11=1,  lOO 
ARRAY  (1,0)  =  0 

1  CONTINUE 
C 

C  ZERO  OUT  SEVERAL  ARRAYS  FOR  INITIALIZATION  PURPOSES 

C 

DO  2  I  =  1.  256 
PACK ( I )  =0 
SUM ( I )  =0 
TOT  < I )  =  0 

2  CONTINUE 
C 

C  I TOT  SHOULD  BE  THE  TOTAL  NUMBER  OF  POINTS  IN  SUM 

C  I COUNT  IS  COUNTER  INCREMENTED  EACH  TIME  A  DOT  IS  DISPLAYED 

C  I BLOCK  IS  CURRENT  BLOCK  BEING  PROCESSED  AND  IS  COMPARED 

C  AGAINST  I  ELK  WHICH  IS  EQUAL  TO  PACK',  1)  OF  THE  BLOCK  JUST 

C  READ.  THIS  COMPARISON  IS  DONE  FOR  VALIDATION  PURPOSES. 

C  I END  IS  FLAG  SET  AFTER  BLOCK  20  IS  READ  TO  PREVENT  READING 

C  PAST  VALID  DATA.  I END  IS  SET  TO  1  AFTER  BLOCK  20  IS  READ. 

I TOT  =  0 
I COUNT  =  0 
I BLOCK  =  0 
I END  =  0 
C 

C  THIS  PROGRAM  SHOULD  BE  RUN  FROM  THE  DOC  TERMINAL. 

C 

CALL  0PEN(1,  ”$TT01"»  O,  IER,  123 > 

TYPE  *'  ENTER  FILENAME  UP  TO  10  CHARACTERS  .  *' 

READdl,  100)  INAMEd) 

100  FORMAT  (S10) 

WRITE  <LP,  1011  INAMEd) 

101  FORMAT  <"  FILENAME  IS  :  ",  S10) 

CALL  FOPEN ( 2,  I NAME,  512) 

TYPE  "  FILE  IS  OPEN" 

C 

C  LOCATE  TO  RELATIVE  BLOCK  0  FOR  HISTOGRAM  INFO 

C 

CALL  FSEEK (2,  O' 

C 

C  READ  RELATIVE  BLOCK  o  INTO  SUM 

C 

READ  BINARY <2 )  SUM 
C 

C  CHECK  TO  SEE  IF  1 0000  POINTS  IN  BLOCK  O 

C 

DO  3  I  =  1,  256 
I  TOT  =  I  TOT  +  SUMU) 

3  CONTINUE 

IF  <  ITOT.  NE.  10000)  GO  TO  990 
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OVERLAY  overo 
SUBROUTINE  UNPACK 

THIS  PROGRAM  UNPACKS  A  FILE  OF  SCANNED  DATA  PRODUCED 
BY  SUBROUTINE  SCAN2  AND  PLACES  IT  IN  A  100  BY  100  ARRAY. 
CHECKS  ARE  MADE  DURING  PROCESSING  TO 

MAKE  SURE  THE  FILE  IS  VALID  AND  THE  THRESHOLDS  ARE  CORRECT. 
THE  NUMBER  OF  POINTS  WHICH  WILL  BE  PLOTTED  IS  GIVEN  ON  THE 
CONSOLE  AFTER  THE  THRESHOLDS  ARE  TYPED  IN.  THIS 
TOTAL  IS  THEN  PROGRAMATICaLLY  CALCULATED  AND  LATER  CHECKED 
AGAINST  THE  FIRST  NUMBER  TO  VALIDATE  PROCESSING. 

SUBROUTINE  LOGIC  IS  USED  TO  UNPACK  DATA  FROM  BLOCKS  2-21. 
SUBROUTINE  SPDOT  IS  USED  TO  PLOT  "  ON  THE  SCREEN  OF  THE 
TEKTRONIX  4006-1  TERMINAL. 


ARRAY  WILL  CONTAIN  THE  SCANNED  VALUES  OF  EACH  BLOCK  WHICH 
WILL  BE  UNPACKED  BY  SUBROUTINE  LOGIC. 

PACK  WILL  CONTAIN  THE  RELATIVE  BLOCK  1  THRU  RELATIVE 
BLOCK  20  DATA  WHICH  WILL  BE  READ  FROM  THE  FILENAME. 

I HOLD  WILL  BE  USED  IN  THE  UNPACKING  DO  LOOP 

SUM  WILL  HOLD  THE  HISTOGRAM  DATA  AND  IS  USED  TO  VALIDATE 

THE  FILENAME. 

TOT  WILL  BE  USED  TO  ACCUMULATE  TOTALS  TO  CHECK  AGAINST 
THE  VALUES  IN  SUM  TO  MAKE  SURE  EVERYTHING  IS  OK 

SIZE  IS  THE  MULTIPLE  OF  A  100  X  100  SCAN  FILE 
WHICH  WILL  BE  DISPLAYED  ON  THE  SCREEN.  THUS  A 
SIZE  OF  4  WILL  GIVE  A  400  X  400  CARTOON  ON  THE  SCREEN. 

UTH  IS  THE  UPPER  THRESHOLD  TAKEN  FROM  THE  CONSOLE 
STEP  IS  THE  WAY  THE  IMAGE  WAS  ORIGINALLY  SCANNED  BY  SCAN2 
AND  IS  NEEDED  IN  ORDER  TO  UNPACK  THE  DATA  INTO  ITS  RELATIVE 
SCANNED  POSITION. 

I NAME  IS  THE  NAME  OF  THE  FILE  TO  BE  PROCESSED. 


COMMON  ARRAY, SUM, LP, IMIN, IMAX, UMIN, UMAX, NCNT , NAREA, AMAX 
I NTEGER  ARRAY ( 1 00,  1 00 ) ,  PACK ( 256  ) ,  SUM ( 256 ) ,  TOT ( 256 ) 
INTEGER  AMAX (5, 20) 

INTEGER  SIZE,  UTH,  UTHLOC,  STEP 
DIMENSION  I NAME <  6 ) 

EQUIVALENCE  (AMAX, I NAME) 


INITIALIZE  ARRAY 
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IF  (IV  GE.  LTHF  AND.  IV.  LE  IJTHF  )  GO  TO  76 
IF  (NS.  EQ.  S  )  GO  TO  65 


1  =  1  F—  1 

IF  (  I.  LT.  1  )  GO  TO  34 
J=  IP+1 

I F ( J.  GT  ICO) GO  TO  25 
I V= ARRAY  < I, J) 

IF (IV.  GE.  LTH.  AND.  IV.  LE.  UTH)GO  TO  20 
IF  (  IV.  GE.  LTHF.  AND  IV  LE.  UTHF )  GO  TO  76 
IF  (NS.  EG.  1  >  GO  TO  65 
GO  TO  25 


COMPLETED  CIRCLE  WITH  NO  FIND.  SEE  IF  WE  ARE  STILL  ON 
FIRST  POINT  IDENTIFIED. 

CONTINUE 

IF  (IP.  EQ.  I  ROW.  AND.  JP.  EQ.  JCOL )  GO  TO  66 
TYPE  "BOY  ARE  WE  LOST!" 

GO  TO  77 


STILL  ON  FIRST  POINT  SO  CONSIDER  IT  AS  ISOLATED  BY 

SETTING  VALUE  TO  ZERO 

ARRAY  (  I  ROW,  JCOL)=V 

PERIM=0.  0 

GO  TO  77 

END  OF  LOOP  FLAG  STARTING  POINT  AND  CONTINUE 
CONTINUE 

IF<  I.  NE.  IROW.  OR.  J.  NE.  JCOL)  GO  TO  10 


CONTINUE 
MT0UCH=0 
RETURN 
MTOUCH= 1 


!  t--  ■  IF'  EO  1  AND  NS.  £U  2  > GO  TO  c-5 

CONTINUE 
._!=•  IP-  1 

IFv-J  LT.  1  >00  TO  45 
1  =  IP 

I  V® ARRAY  (  I .  .J) 

IF '  IV  OE.  LTH.  AND.  IV.  L E.  UTH'GO  TO  20 
IF  I  IV.  OE  LTHF  AND  IV  LE  UTHF ) GO  TO  76 
IF  (NS.  £Q.  4  .■  GO  TO  65 


I=IP+1 

IF (I  GT  100, ’GO  TO  55 
J=JP-1 

I F  <  J.  LT  1  >  GO  TO  45 
IV-ARRAy ( I . J) 

IF  (IV.  GE.  LTH.  AND.  IV.  LE.  UTH)GO  TO  20 
IF(  IV.  GE  LTHF.  AND.  IV.  LE.  IJTHF  ) GO  TO  76 
IF  (NS.  EO.  5>  GO  TO  65 


I  =  IP+l 

IF(  I.  GT  100) GO  TO  55 
U=.JP 

IV® ARRAY ( I,  J) 

IF  (IV  GE.  LTH  AND.  IV  LE.  UTH)GO  TO  20 
IF<  IV  GE.  LTHF.  AND.  IV.  LE.  IJTHF) GO  TO  76 

IF  (  NS.  EQ.  6  )  GO  TO  6 5 


I  =  I P+ 1 

IF  (I.  GT.  100)  GO  TO  55 
J= JP+ 1 

IF  (  J.  GT.  100) GO  TO  25 
IV® ARRAY ( I,  J) 

IF  (IV.  GE.  LTH.  AND.  IV.  LE.  UTH.'GO  TO  20 
I F  (  IV.  GE.  LTHF.  AND  IV.  LE.  IJTHF  >G0  TO  76 
IF  (NS.  EQ.  7) GO  TO  65 


J=JP+1 

I F  (  U  GT.  1  uO )  GO  TO  25 
I  =  IP 

I V® ARRAY ( I,  J) 

I F  (  IV.  GE.  LTH.  AND.  IV.  LE  UTH.’GO  TO  20 


GO  TO  IS 

125  IF ( I- IOi  126,  127,  123 

126  NS=7 

GO  TO  13 

127  NS =3 

GO  T 0  1 3 
123  NS= 1 

C 

C  THROUGH  WITH  OLD  POINTS.  SAVE  CURRENT  POINTS  FOR  USE  NEXT  Ti 

C 

18  CONTINUE 

IO=I 
■JO=.J 

GO  TO  <  25,  30,  35,  40,  45,  50.  55,  60  > ,  NS 
C 

C  NOW  WE  KNOW  WHERE  TO  START,  CIRCLE  CLOCKWISE  LOOKING  FOR 

C  NEXT  POINT  IN  AREA  BOUNDARY. 

C 

25  CONTINUE 

C 

C  TO  PREVENT  LOOPING  FOR  SPECIAL  CASE  OF  ISOLATED  POINT 

C  ON  ROW  1  (1=1).  NEED  FOLLOWING  CHECK. 

IF  (IP.  EQ.  1.  AND.  NS.  EG  2. 'GO  TO  65 
C 

I=IP-1 

IF <  I.  LT  1  )G0  TO  34 
J=.JP 

I V= ARRAY! I, J) 

IF <  IV.  GE.  LTH.  AND.  IV.  LE.  UTHjGO  TO  20 
IF(  IV.  GE.  LTHF.  AND  IV  LE.  IJTHF )  GO  TO  76 
IF  (NS.  EQ.  2) GO  TO  65 
C 
C 
C 

30  CONTINUE 

I=IP-1 

IF (  I.  LT.  1  )  GO  TO  35 
J=JP-1 

I F  (  J.  LT  1  )  GO  TO  45 
IV=ARRAY( I, J) 

IF ( IV.  GE.  LTH.  AND.  IV.  LE  UTH)GO  TO  20 
IF (  IV.  GE.  LTHF.  AND.  IV.  LE.  UTHF ) GO  TO  76 

IF  (NS.  EQ.  3)  GO  TO  65 
C 
C 
C 

34  CONTINUE 

C 

C  DITTO  COMMENT  AT  25 
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CALL  CKF'T  (  I .  J.  1 1 .  J 1 .  MTOUCH,  LTH.  UTH ) 

IF  (  10.  GT  I  M2  >  GO  TO  29 

IF<  I  EQ.  I  M2.  AND.  J.  NE.  JM2>G0  TO  29 

CALL  CKF'T  (10.  JO.  II.  J 1 .  MTOUCH,  LTH,  UTH ) 

IF  (  J  NE.  JM2  )  GO  TO  29 
11=0 
J 1 =-JE 

CALL  CKF'T  (  10,  JO,  1 1 ,  J 1 ,  MTOUCH,  LTH,  UTH  ) 

IF  (  I  M2.  LT.  10)  GO  TO  29 

I 1=— ID 

J1=-JD 

CALL  CKPT ( 10, JO.  1 1 ,  J 1 ,  MTOUCH,  LTH, UTH ) 

IF  (  I.  NE.  I  M2 )  GO  TO  29 
J 1  =0 
I 1=-I£ 

CALL  CKF'T  (  10,  JO,  II,  J 1 ,  MTOUCH,  LTH,  UTH  ) 

GO  TO  29 

29  IF  (MTOUCH.  NE.  0)G0  TO  80 

C 

C  HAVE  WE  BEEN  FULL  CIRCLE 

C 

IF  (I  EQ.  I  ROW.  AND  J.  EQ.  JCODGO  TO  77 
C 

C 

21  CONTINUE 

IM2=IF' 

JM2= JP 
I  F'=  I 
JF'=J 
C 

C  FOUND  AND  PLOTTED  A  POINT.  SEARCH  FOR  NEXT  POINT  BEGINS. 

C  ASSUME  THE  SEARCH  OF  ADJACENT  POINTS, THE  8  POINTS  TO 

C  BE  EXAMINED  ARE  NUMBERED  1  TO  8  STARTING  WITH  THE  ONE 

C  AT  THE  I - 1 , J  POSITION  AS  1  AND  MOVING 

C  CLOCKWISE.  FOR  THE  POINT  WE  ARE  CURRENTLY  AT,  THERE 

C  WILL  BE  A  PREVIOUS  POINT  AS  ONE  OF  THE  8.  OUR 

C  EXAMINATION  WILL  START  WITH  THE  NEXT  CLOCKWISE  POINT 

C  AND  PROCEED. 

C 

IF(J-JO)  125,  120,  115 

115  IF(I-IO)  118,  117,  116 

116  NS =3 

uO  TO  l y 

117  NS=4 

uO  Tu  18 

1 1 8  NS =5 

GO  TO  18 

1 20  NS* 2 

IF  (  I  LT.  1 0 ) NS* 6 
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GO  TO  29 


IF ( 1-10)16. 17, 19 


1 1  =  I D 
J 1  =  JD 

CALL  CKPT ( I , J,  1 1 , J1 , MTOUCH, LTH, UTH ) 

IF ( 10  LT  I M2 ) GO  TO  29 
IF  ( J.  EQ.  JM2.  AND.  I.  NE  I  M2)  GO  TO  29 
CALL  CKPT ( 10, JO, II, J 1 , MTOUCH, LTH, UTH ) 
IF  (  I.  NE.  I  M2 )  GO  TO  29 
J 1  =0 
1 1  =  IE 

CALL  CKPT ( 10, JO,  II,  J 1 ,  MTOUCH. LTH, UTH ) 
IF  (JO.  GT.  JM2 )  GO  TO  29 

11  =  10 
J1=-JD 

CALL  CKPT < 10, JO,  II, J1 , MTOUCH, LTH, UTH ) 
IF (  10.  NE.  I  )GG  TO  29 
11=0 
J1=-JE 

CALL  CKPT ( 10, JO,  1 1 ,  J 1 ,  MTOUCH,  LTH, UTH  > 
GO  TO  29 


Jl-0 
1 1  =  IE 

CALL  CKPT ( I,  J,  1 1,  Jl,  MTOUCH,  LTH,  UTH  ) 

IF(  10.  LT.  I  M2  J  GO  TO  29 

IF(  I.  EQ.  I  M2.  AND.  J.  NE.  JM2)G0  TO  29 

J1=0 

CALL  CKPT ( 10, JO, 11, Jl, MTOUCH, LTH, UTH) 
IF ( 10  GT.  I M2 ) GO  TO  29 
1 1  =  ID 
Jl  =-JD 

CALL  CKPT ( 10,  JO,  II, Jl, MTOUCH, LTH, UTH) 
I F  <  J.  NE.  JM2 )  GO  TO  29 
11=0 
Jl =-JE 

CALL  CKPT ( 10, JO.  1 1 , Jl , MTOUCH,  LTH,  UTH ) 
IF  (  I.  NE.  I  M2  )G0  TO  29 
1 1 =- 1 D 
Jl =-JD 

CALL  CKPT ( 10, JO, II, Jl, MTOUCH, LTH, UTH) 
GO  TO  29 


1 1  =  ID 
J1=-JD 
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CALL  CKPT ( 10. JO.  1 1 . J 1 , MTOUCH, LTH. UTH ) 
IF  <  .JO  EC!.  JM2  >  00  TO  29 
11=0 
J1=JE 

CALL  CKPT ( 10. JO,  1 1 , Jl, MTOUCH,  LTH,  UTH ) 
GO  TO  29 


IF ( I.  LT.  10)00  TO  11 
11=0 
J 1 =- JE 

CALL  CKPT < I , J,  1 1 ,  Jl ,  MTOUCH, LTH, UTH  > 

IF v  JM2.  LT.  J ) GO  TO  29 
IF  ( J.  EQ.  JM2.  AND.  I.  NE.  I  M2)  GO  TO  29 
11=0 

CALL  CKPT ( 10, JO,  1 1 , Jl , MTOUCH, LTH, UTH ) 
IF  (I  M2.  LT.  10)  GO  TO  29 
I 1=-ID 
J1»-JD 

CALL  CKPT ( 10, JO,  1 1 ,  Jl ,  MTOUCH,  LTH,  UTH ) 
IF <  IM2.  LT.  I  ) GO  TO  29 
1 1=-IE 
J 1  =0 

CALL  CKPT < 10, JO, II, Jl, MTOUCH, LTH, UTH) 
IF ( JM2.  GT.  J ) GO  TO  29 
I 1=— ID 
J1=JD 

CALL  CKPT ( 10, JO, II, Jl, MTOUCH, LTH, UTH) 
GO  TO  29 


11=0 

J1=JE 

CALL  CKPT ( I , J,  1 1 , Jl , MTOUCH, LTH, UTH ) 

IF ( J.  LT.  JM2 ) GO  TO  29 
IF  ( J.  EQ.  JM2.  AND.  I.  NE.  I  M2)  GO  TO  29 
11=0 

CALL  CKPT < 10, JO,  1 1 , Jl , MTOUCH, LTH, UTH ) 
IF(  I  M2  GT.  10)  GO  TO  29 
1 1  =  I D 
J1=JD 

CALL  CKPT < 10, JO,  1 1 , Jl , MTOUCH, LTH, UTH ) 
IF  (  I  M2.  GT.  I  )  GO  TO  29 
J 1  =0 
1 1  =  IE 

CALL  CKPT < 10, JO,  1 1 , Jl , MTOUCH, LTH, UTH) 
IF< I  NE  I M2) GO  TO  29 
1 1  =  I D 
J1=-JD 

CALL  CKPT ( 10, JO,  1 1 , Jl , MTOUCH, LTH, UTH ) 


IF  (JO  EG!.  JM2  )  GO  TO  29 
1 1  =0 
J 1  =  JE 

CALL  CKPT ( 10, JO, II. J 1 , MTOUCH, LTH, UTH ) 
IF  (10.  LT.  I  M2)  GO  TO  29 
1 1  =  ID 
J 1  =  JD 

CALL  CKPT  (  10,  JO,  1 1 ,  ,J  1 ,  MTOUCH,  LTH,  UTH  > 

IF  ( I.  NE.  I  M2)  GO  TO  29 

•J1=0 

1 1  =  I E 

CALL  CKPT ( 10, JO,  1 1 ,  J 1 ,  MTOUCH,  LTH,  UTH  > 
GO  TO  29 


J1=0 
I 1=— IE 

CALL  CKPT ( I, J, II, J 1 , MTOUCH, LTH, UTH ) 

IF  (  I  M2.  LT.  I  )  GO  TO  29 

IF  (I.  EQ.  I  M2.  AND.  JO.  LT.  JM2)G0  TO  29 

J 1  =0 

CALL  CKPT (10, JO. II, Jl, MTOUCH, LTH, UTH) 
IF(  JM2.  GT.  10) GO  TO  29 
1 1  =—  I D 
•J 1  = -JD 

CALL  CKPT < 10, JO,  1 1, Jl, MTOUCH, LTH, UTH) 

IF  ( JM2.  NE.  J ) GO  TO  29 

11=0 

Jl-JE 

CALL  CKPT ( 10, JO,  1 1 , Jl, MTOUCH, LTH, UTH) 
IF  (  I.  NE.  I  M2 )  GO  TO  29 
1 1  =  ID 
J1=JD 

CALL  CKPT ( 10, JO,  1 1 , Jl , MTOUCH, LTH, UTH ) 
GO  TO  29 


I 1=- ID 
Jl  =-,JD 

CALL  CKPT ( I , J,  1 1 , Jl , MTOUCH, LTH, UTH ) 

IF  (  10.  GT  I  M2)  GO  TO  29 

IF ( J.  EG!.  JM2  AND.  I.  NE  IM2)G0  TO  29 

CALL  CKPT ( 10, JO,  1 1 , Jl , MTOUCH, LTH, UTH ) 

IF  <  I  GT.  I  M2)  GO  TO  29 

I 1=-IE 

J1=0 

CALL  CKPT ( 10, JO,  1 1 , Jl , MTOUCH, LTH, UTH ) 
IF  (JO.  LT.  JM2 )  GO  TO  29 


TITLE  TSPOT 


PREPARED  FOR  SC  I  St STEMS.  INC. 

BY  EILL  POPE,  TELCOM  DATA  CORFORTION. 


THIS  SUBROUTINE  PROVIDES  FO  A  RAPID  100*100  SCAN 
AUTOMATICALLY'  PLACING  THE  VALUES  INTO  THE  STORAGE  ARRAY 
SO  ONLY  ONE  CALL  TO  THE  SUBROUTINE  IS  NEEDED. 

.  ENT  TSPOT 

.  EXTN  .  U I E X 

.  EXTU 
.  NREL 

INTEGER  ARRAY  v  lOV,  1*.'VJ 

CALL  TSPOT ( I START X, I ST ARTY. ARRAY) 


ISX 

ISY 

ARRAY 

IX 

IY 

COUNT 
COLCNT  = 
YCNT 
LOC 

FS. 

VIDEO 


-167 

ISx  +  1 

ISY+1 

ARRAY+1 

IX+1 

IY+1 

COUNT  + 1 

COLCNT + 1 

YCNT+1 

LOC-ISX+l 

70 
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SUBROUTINE  HOLE  <  IN,  ON-  10.  00.  LTH,  UTH,  PERIM,  NBIAS,  B,  NBP ) 

FREF'AFED  F OR  SCI  SYSTEMS •  INC. 

By  BILL  POPE.  TELCOM  DaTh  CORPORATION. 

WHEN  WORr  I  NO  WITH  AN  AREm  hNO  FILLING  IN  THE  CENTER  , 
CALLED  TO  TRy  AND  TRACE  THE  EDGES  OF  THE  "HOLE",  FIND  THE 
PER I MEIER  AND  AREA  OF  THE  HOLE,  AND  FIND  THE  MINIMUM 
SURROUNDING  DISTANCE. 


VARIABLES  INCLUDE 


I  N  ■  ON 

10.  00 

LTH, UTH 
PERIM 
NBIAS 
B< .2, NBP ) 


THE  CURRENT  POINT  WITHIN  THE  SEGMENT  WHICH  SHOULD  BE 
THE  FIRST  POINT  ON  THE  EDGE  OF  THE  HOLE. 

THE  POINT  FOUND  WHICH  WAS  OUTSIDE  THE  THRESHOLDS  AND 
MADE  US  EXPECT  TO  FIND  A  HOLE. 

THRESHOLDS  OF  THE  CURRENT  SEGMENT. 

PERIMETER  OF  THE  HOLE,  TO  BE  RETURNED. 

CURRENT  SEGMENT  NUMBER,  I.  E.  VALUE  USED  TO  MARK  POINT 
ARRAY  CONTAINING  COORDINATES  OF  CURRENT  SEGMENT. 


COMMON  ARRAY,  SUM,  LP,  IMIN,  IMAX,  .JMIN,  JMAX,  NCNT ,  NAREA,  AMAX 
COMMON  /PL T/ 1 START X,  I ST ARTY,  MIDx,  MIDY, NXSIZE, NYSIZE 
I NT EGER  ARRA  Y  < 1  oO .  1 00 > , SUM <  256 ) , AMA  X ( 5 , 20 ) 

INTEGER  B(2,  1500) ■ Cl  2,  500) 

INTEGER  IJTH.UTHF 


FIRST,  CALL  TRACE  TO  FIND  EDGE  OF  HOLES. 

PI  =3.  1415926 
I  M=  I N 
JM=JN 
I  MO*  10 
JMO* JO 

SET  NSF'L  TO  MARK  EDGE  OF  THIS  HOLE.  MARK  IT  HIGHER  THAN  NORMAL 
SO  THAT  WE  CAN  DISTINGUISH  IT  FROM  OTHER  HOLES  IN  AREA. 

NSPL=76S 

CALL  TRACE  (  IN.  UN,  10,  JO,  LTH,  UTH,  PERIM,  NBIAS,  NSF'L,  IMN,  IMX,  JMN,  JMX, 
C, NCR > 

IF  l  PERIM.  LE.  20.  )  RETURN 

QUIT  IF  SPOT  LESS  THAN  20  AROUND 
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WE  HAVE  A  HOLE ■  AND  FROM  SUBROUTINE  TRACE  KNOW  THE 
MINI  MUM  AND  MAXIMUM  VALUES'  FOR  THE  ROWS  AND  COLUMNS  THAT 
IT  FALLS  BETWEEN.  NOW.  FOP  EACH  ROW  AND  EACH  COLUMN,  FIND 
THE  MINIMUM  AND  Ma  < I MUM  HOLE  BOUNDARY  POINT  Br  TAKING 
THE  AVERAGE  OF  THE  ROWS  AND  THE  AVERAGE  OF  THE  COLUMNS, 
FIND  THE  CENTER  OF  THE  HOLE . PRETEND  I NG  IT  IS  A  CIRCLE). 


IC=0 
NC  =0 

FIND  THE  MARKED  THRESHOLD  LOWER  AND  UPPER  LIMITS  OF  THE  HOLE 

MTL=LTH+NE I AS+NSPL 
MTU=UTH+NB I AS+NSPL 
DO  6  JMX 

SET  THE  LEFT  uND  RIGHT  LIMITS  TO  ZERO 

IVS=-1 

IWS=-1 

I u= I M  <  +  1 

DO  4  I  ~  I  Mi;  ■  I  M  e 

AS  WE  SEARCH  FROM  l.EF  i  TO  RIGHT  USING  INDEX  I,  WE  WILL  SEARCH 
FROM  RIGHT  TO  LEFT  US IN  INDEX  lu 

IU=!U-l 

IS  LEFTMOST  HOLE  BjJNOaRy  IDENTIFIED 

Ici  I  VS  NE.  -  i  :  GO  TO  1 
I  V=ARPa  r  v  I,  J  ;• 

IF  IV  LT  MTL  'GC  T$  : 

IF'  IV  GT  MTU  GO  TO  1. 

FOUND  LEFTMOST  BOUNDARY  VALUE  IN  IV), 

SAVE  LOCATION  IN  I VS. 


LOOK  FOR  RIGHTMOST  BOUNDARY 

IF  IWS.  NE  -1  ) GO  TO  2 
I W= ARRAY < IU, J) 

I F  <  I W  L T  M T L ) GO  TO  2 
IF  '  IW  GT.  MTU  1  GO  TO  2 
IWS=IU 

IF'  I  VS  NE  -1  AND.  IWS  NE.  - 1  .00  TO  5 
CONTINUE 


•  I  LED  TO  FIND  HOLE  dOUNDaF  i . 
0T0P  HOLERF: 


CALL ULhTE  -'Ul'l  OF  THE  AVERAGE  LOWER  ANE'  UPPER  EDGE  F'OINTS. 

I  C  =  1 G+  (  I  \'S+  I WS  ) /  Z 
NC=NC+ 1 
CONTINUE 

HAVE  SUM  OF  EDGE  POINTS  AND  HOW  MANY  FOUND.  GET  AVERAGE 
I C= I C/NC 

IC  NOW  HAS  X  COORDINATE  OF  THE  HOLE  CENTER.  WILL  REPEAT 
PROCESS  ALONG  COLUMNS  TO  GET  V  COORDINATE 


uc=o 

NC=0 

DO  11  I=IMN. IMX 

•JVS=- 1 

JWS=- 1 

I IJ— JM  X  + 1 

DO  9  J=.JMN»  JMX 

IU=IU-1 

IF  (  UVS.  NE.  -1  ) GO  TO  7 
I V= ARRAY',  I,  J) 

IF  v  IV.  LT.  MTL  )  GO  TO  7 
IF<  IV.  GT.  MTU > GO  TO  7 
JVS=.J 

IF  (  JWS.  NE  -1  l  GO  TO  8 
I W=ARRAY  ;  I .  I  IJ  > 

IF  (  IW.  LT  MTU  GO  TO  8 
IF IIW.  GT.  MTU) GO  TO  8 
JWS=IU 
CONTINUE 

IF ( JVS.  NE.  -1.  AND.  JWS.  NE.  -l)GO  TO  10 

CONTINUE 

STOP  H0LER2 

.JC= JC+  <  JVS+ JWS ) /2 

NC=NC+ 1 

CONTINUE 

JC=JC/NC 


THE  POINT  IC. JC  REPRESENTS  THE  CENTER  OF  OUR  HOLE.  IF  IT  WAS 
ROUND. 


NOW  WE  WANT  TO  FIND  THE  "RADIUS"  FROM  IC. JC  TO  EACH  POINT  ON 
THE  CIRCUMFERENCE  OF  THE  HOLE.  ALSO  FIND  THE  MINIMUM  AND  MAXIMUM 
RADIUS  DEFINED  MOST  IMPORTANTLY .  FIND  THE  MINIMUM  THICKNESS 
OF  THE  PAD*  AROUND*  THE  HOLE. 

FIRST.  SET  THRESHOLD  LIMITS  FOR  THE  SEGMENT  CONTAINING  THE  HOLE. 
SECOND,  SET  LIMITS  FOR  THE  BOUNDARY  OF  THE  SEGMENT  CONTAINING  HOLE. 

MTL=LTH+NB I AS 
MTIJ=MTL+NAREA 
MEL = NB I  AS +256 
MEU=MEL+256 


XMIN*  175. 

RADTOT=v. 

RADL=0 

RADS-400. 

FOR  EACH  POINT  ON  CIRCUMFERENCE,  FIND  RADIUS. 

DO  15  NC= 1 , NCR 
IP=C(  1,  NO 
JP=C<  2  ,  NC  > 

RI-1  25*FLOAT t IF- 1C > 

RJ=FLOAT<  jP-JC.- 

RADI  US = SORT  (  RI**2+RJ**2  > 

WR I TE ( LP ,  247) IP, JP , RaD I US 

47  FORMAT  (  T20,  "PT=",  13,  ",  ",  13,  2X,  "RADIUS*",  F6.  2> 
RADTOT =RADTOT +RAD I US 
IF ( RADIUS.  GT.  RADL>RaDL*RAOIUS 
IF  (RADIUS.  LT.  RADS )  RADS=RAD  I  US 
IF  (RADIUS.  LT.  1.  U)GO  TO  15 


NOW  WE  HAVE  RADIUS.  TO  FIND  MINIMUM  THICKNESS.  PROJECT  ALONG 
THE  RADIUS  A  DISTANCE  OF  XMIN.  XMIN  IS  THE  PREVIOUSLY  FOUND 
MINIMUM.  IF  STILL  IN  THE  SEGMENT,  THEN  THERE  IS  NOT  A  NEW 
MINIMUM.  IF  PROJECTION  FALLS  OUTSIDE  THE  MINIMUM,  THEN  WE 
NEED  TO  POINT  BY  POINT  MOVE  OUTWARD  ALONG  RADIUS  DIRECTION 
UNTIL  WE  GET  TO  EDGE  OF  SEGMENT.  THAT  POINT  SHOULD  BE  NEW 
MINIMUM. 


DEPENDENT  ON  WHETHER  I  OR  J  VECTOR  LARGER,  SET  ABSOLUTE  RATIO 
OF  ONE  TO  THE  OTHER. 

I F ( ABS ( R J ) .  GT  AES  >  R I >  » GO  TO  1 6 
XDI=RI/ABSt PI ) 
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XDJ=RJ.'  AE"S  v  R I  ' 

GO  TO  17 
XDI-RI / ABSvR.J 
XDJ=R J/ AES t  R J ' 

CONTINUE 

IF  THERE  IS  NO  MINIMUM  FOUND  VET,  USE  POINT  TO  POINT  MOVE. 

IF  (  XMI N  EO.  175.  *GO  TO  5«.- 
ID=IP+IFIX( XDI*XMIN+.  5) 

IF  i  ID.  LT.  1  « ID=1 
I F  <  ID.  GT,  1  VO )  I D=  1 
UD=.JP+ 1 F I  X  (  XDJ*  XM I N+.  5  > 

IF(JD.  LT  1  ) JD=  1 
I F '  JD.  GT.  1 00 )  JD=  1  '_>>.* 

ID,  JD  IS  THE  PROJECTION  OF  THE  PREVIOUSLY  FOUND  MINIMUM  VALUE 
ALONG  THE  DIRECTION  OF  THE  CURRENT  RADIUS. 

I X  =  I ST ART  X+ 1 D*NXS I ZE 
I Y= I STARTY- JD*NYS I ZE 
CALL  SPDOT (IX. I Y ) 

GET  VALUE  OF  POINT  TO  FIND  OUT  IF  IT  IS  STILL  IN  SEGMENT. 
IV=ARRAYv ID.  JD) 

IF (  IV.  GE.  LTH.  AND.  IV  LT.  UTH)GO  TO  15 
IF v  IV  GE.  MTL.  AND.  IV.  LT.  MTU) GO  TO  .15 


PROJECTION  IS  OUTSIDE  SEGMENT  SO  MUST  NOW  MOVE  ALONG  PROJECTION 
VECTOR  TIL  WE  FIND  NEW  MINIMUM. 

XDI  AND  XDJ  CONTAIN  THE  " DELTA V"  TO  MOVE  IN  THE  I  AND  J 
DIRECTION  ONE  OF  THEM  WILL  HAVE  A  VALUE  OF  1.  U  AND  THE 
OTHERS  VALUE  WILL  BE  LESS  THAN  1  0.  THAT  WAY,  WHEN  WE 
MAKE  A  "MOVE",  WE  WILL  ALWAYS  GO  ONE  DIRECTION  AT  LEAST 
1.  0  UNITS.  BY  SAVING  THE  "REMAINDER”  NOT  MOVED  AND  SUMMING 
IT,  WHEN  THE  OTHER  DIRECTIONS  VALUE  EXCEEDS  1.  u,  THEN 
WE  CAN  MOVE  IN  THAT  DIRECTION  ALSO.  THIS  WILL  GIVE  US 
AS  APPROXIMATE  DIRECTED  MOVEMENT. 


RMI=0  0 
RM  J=|.'  0 
I L  =  I P 
JL= JF 

10=  IF  I  X<  XDI+RMI  ).,  1 
RM I  =RM  I  -*•  X  D I 

IF  ABSiRMl/  LT  1  >00  TO  53 


B-t.S 


RM I  =RM I -FLOAT  ■;  ID  > 

GO  TO  52 

•JD=  IF  I  '  i  xDJ+RMJ  i  /  1 
RMJ=RMJ+XDJ 

i  F  ■ ■  ABS  <  RMJ  »  LT  1  >  00  TO  55 

F:MJ=F:MJ-FLOAT '  JD  ) 

GO  TO  54 
CONTINUE 

ID  AND  JD  ARE  THE  NEW  DELTAS.  MAKE  A  MOVE 

IL=IL+ID 

JL-JL+JD 


IF(  IL.  EQ.  1  ) GO  TO  15 

IF ( JL.  EQ.  1  ) GO  TO  15 

IF( IL  EQ.  100)00  TO  15 

I F  ( JL.  EQ.  100)  GO  TO  15 

PLOT  SPOT  JUST  MOVED  TO,  IF  DESIRED. 

I X= I START X+ IL*NXSI ZE 
IY=ISTARTY-JL*NYSIZE 
CALL  SPDOT < I  X ,  I Y ? 

IS  NEW  POINT  ON  EDGE  YET 

I V=ARRAY ( IL, JL) 

IF  (  IV.  GT.  MEU ) GO  TO  51 
IF(  IV  GE.  MEL) GO  TO  56 
IF (  IV.  GE.  LTH.  AND.  IV.  LT  UTH)GO  TO  51 

DID  WE  INADVERTENTLY  PASS  THROUGH  EDGE  ON  DIAGONAL 

IF < IV  GE.  NAREA ) GO  TO  51 

YES,  LOOK  AT  ADJACENT  POINTS  TO  TRY  AND  FILD  EDGE. 

1 1  =  IL 
J1=JL 


IL= 1 1 - 1 

IF ( IL  LT  1 ) GO  TO  151 
IV=ARRAY< IL, JL  > 

IF <  I V  GE  MEL  AND  IV  LT.  MEU ) GO  TO  56 

JL-Jl-1 

IL=  I  1 

IF(JL  LT  1 >G0  TO  152 
I V=ARRAY < IL, JL  > 

IF( IV  GE  MEL  AND  IV  LT  MEU)GO  TO  56 


IL=I 1+1 
JL  =  J 1 

IF  IL  GT  1  ;  GO  TO  1 50 

I GRF'i-i  f  IL  JL  ■■ 

IF'  10.  GE  MEL  AND  IV  LT  MEU ' GO  TO  56 
JL= J 1 + 1 
IL=  1 1 

I F  i  JL  GT  1  *_<«.■ .  GO  TO  15 
I V= ARRAY'S  IL.  JL) 

IF'  IV  LT  MEL  OR.  IV.  GE  MEG  .'GO  TO  15 


CONTINIJE 

CALCULATE  DISTANCE  AND  CHECK  AGAINST  MINIMUM. 

*01  =  1  25-»FL0rtT<  Il-IF'J 
XDU=FLOATv JL-JFN 
xD=SORT  v  xDl**l+xDu**2  ) 

IF <  XD.  GE.  XMIN  'GO  TO  15 

XMIN=*D 

I  E=  I P 

•JB= JP 

IK=IL 

JK»  JL 

CONT I NUE 


WRITE  S  LP,  1  '.'i.' )  NCNT ,  XMIN 

FORMATS  "'.'HOLE  IN  AREA  ",  IS/  1UX ,  "MINIMUM  THICKNESS  AROUND  HOLE  IS 
Ft.  2) 


NOW  PLOT  LINE  ALONG  MINIMUM  DISTANCE 
CALL  DOTL INS  IE, JB,  IK,  JK ) 


NOW  SCAN  THE  AREA  BETWEEN  THE  MIN  AND  MAX  VALUES  FOR  THE  HOLE. 

FOR  POINTS  WITHIN  THE  HOLE  BOUNDARY,  CALCULATE  AREA.  CALCULATION 
EASED  ON  LOOKING  AT  FOUR  CORNERS  OF  A  SQUPE  AND  COUNTING  HOW  MANY 
HAVE  VALUE  WITHIN  THE  RANGE  OF  THE  AREA  SCANNED.  SINCE  AREA  OF  HOLE 
DEFINED  BY  POINTS  "OUTSIDE"  RANGE,  EITHER  ALL  OR  HALF  A  SQUARE  IS 
WITHIN  HOLE  AREA  IF- 

MN=4  NO  PART  IN  HOLE 

mn=  ?■  half  a  square  in  hole 

MN  !  ?.  ALL  OF  SQUARE  IN  HOLE 
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C  L  THF  AND  IJTH F  SET  AS  LIMITS  OF  HOLE  BOUNDARY. 

LTHF=LTH+NBI  AS+NSF'l 
U  THF-Ll  TH+NB I  AS-TISPL 
JM1=.JMN+1 
RAREA=0.  >.» 

I F  ( -JM 1  G  T  .JM  X  >  00  T  0  32 
DO  31  J=JM1,JMX 
C 

C  INITIALIZE  FLAG  TO  INDICATE  SCAN  IS  STILL  OUTSIDE  HOLE. 

NFL=0 

DO  30  I=IMN, I MX 
IF (NFL.  GT  0)00  TO  25 
C 

C  LOOK  FOR  FIRST  BOUND  AF'1'  OF  HOLE 

IF! ARRAY!  I,  J).  GE.  LTHF  ANu  ARRAY !  I ,  U  L  L£.  UTHF  )  GO  TO  26 
GO  T 0  30 

25  IF  (  ARRAY  <  I,  J).  GE.  LTHF  AND.  ARRAY  (  I,  J).  LE.  UTHF  )  GO  TO  26 

C  PASSING  OUT  OF  HOLE  AREA,  SET  FLAG  BACK 

IFi  ARRAY  !  I,  J).  LT.  NBIaS/GO  TO  26 
NFL=0 
GO  TO  27 

26  NFL=1 

27  NH=0 
C 

C  COUNT  CORNERS,  WATCHING  OUT  FOR  ROW  AND  COLUMN  1 

IF<  I.  EQ.  1  >00  TO  3*.' 

IF(  J.  EQ.  1  >00  TO  3>J 
IF ! ARRAY ( I , J  J .  GE.  NEI AS i NM=NM+ 1 
I F ( ARRAY  < I - 1 , J ) .  GE.  NB I  AS ) NM=NM+ 1 
I F  «  ARRAY  (  I  - 1 ,  J- 1  > .  OE.  NB  I  AS  )  NM=NM+ 1 
IF  (  ARRAY  (  I,  J-l  ).  OE.  NB  I  AS  >  NM=NM+ 1 
IF ( NM  EQ  3 ) PAREA=RAREA+1.'.  *25 
IF  ( NM  LT  3  )  RAREA=F:aREh+  1  25 

30  CONTINUE 

31  CONTINUE 
C 


C 


32 

107 

CONTINUE 

WRITE (  LF\  107)  IC,  JC 

FORMAT ( lux , "HOLE  CENTER-", 13, ", " 

,  13) 

1 05 

RAV=RADTOT/FLOAT <  NCR ) 

WRITE'.LP,  105 )  RAV 

FORMAT ( iOX, "AVERAGE  RADIUS®", F7 

2 ) 

1 06 

WR I TE ( LP ,  1 06 ) RADS, RADL 

FORMAT! IOX, "RADIUS  RANGE  FROM-", 

F7.  2,  "  TO 

",  F7.  2) 

101 

WR I TE ( LP, 101 )PERIM 

FORMAT!"  CIRCUMFERENCE 

OF  HOLE  IS 

",  F7.  2) 
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WRITE  iLF\  1  *.>2 RaREA 

FORMAT*"  AREA  OF  HOLE  IS  ",F7.  2> 

PAT  1 0=4  -*P  I  *RAREA/  * PER  I  M**2 
WRIT£*LP  104 'RATIO 

FORMAT*"  4*PI*aREA/*C**2)=  "«F6.  2) 

WR ITEtLP.  lOSxIMN, I M  X . JMN • JM  X 

FORMAT*"  LOCATION  -  IMN-" .13."  IMX="«I3»" 

13. "  JMX= " i 13' 


JMN=" . 


NOW  THE  ABNORMAL  NSPL  VALUE  AT  FIRST  OF  SUBROUTINE  MUST  BE  BACKED  OUT 
This  WILL  LEAVE  THE  EDGE  OF  ALL  HOLES  BEING  MARKED  WITH  A  VALUE 
OF  NBIAS  +  512 


NSPL=-256 
DO  41  J=JMN,  JMX 
DO  40  I  =  I MN <  I MX 

IF  (ARRAY*  I,  J).  LT.  LTHF)GO  TO  40 

I F  (  ARRAY  (  I ,  J )  GT.  UTHF  )  GO  TO  40 

ARRAY  *  I .  J ) = ARRAY  *  I . J  > +NSPL 

CONTINUE 

CQNT I NUE 

RETURN 

END 


o o  n  o o  no  o  o  o o  o o  o o  o o  o o o o o o 


OVERLAY  07ERB 

SUBROUTINE  INTERNAL  S, NB, XMIN) 


C  PREPARED  FOR  SCI  SYSTEMS,  INC. 

C  BY  BILL  POPE,  TELCOM  DATA  CORPORATION. 

C  THIS  SUBROUTINE  DETERMINES  THE  MINIMUM  INTERNAL  THICKNESS  OF  A 

C  PREVIOUSLY  DEFINED  SEGMENT.  THE  POINTS  BELONGING  TO  THE  SEGMENT 

C  HAVE  BEEN  "MARKED  IN  THE  ARRAY  WITH  A  VALUE  UNIQUELY  IDENTIFYING 

C  THEM. 

C  TO  FIND  THE  MINIMUM,  EACH  POINT  ON  THE  EDGE  WILL  BE  EXAMINED 

C  TO  FIND  THE  NECESSARY  INTERNAL  DISTANCE  FROM  THAT  POINT. 

C  IF  THE  DOCUMENTED  EXPLANATION  OF  THE  SEARCH  PATTERN  IS  NOT 

C  SUFFICIENT,  RECOMPILE  THIS  SUBROUTINE  AFTER  " UN-COMMENT I NG " 

C  THE  STATEMENTS  CONCERNING  PLOTTING  OF  SPDOT.  THE  SEARCH  PATTERN 

C  WILL  THEN  BE  TRACED  IN  A  VISUAL  MANNER. 

VARIABLES  INCLUDE: 

B<2, NB)  ARRAY  CONTAINING  THE  COORDINATES  OF  ALL  POINTS  ON  EDGE 

XMIN  THE  MINIMUM  THICKNESS  FOUND. 


INITIALIZE  BY  SETTING  SEARCH  LIMITS  FOR  AREA 

COMMON  ARRAY,  SUM,  LP,  IMIN,  IMAX,  .JMIN,  UMAX,  NCNT ,  NAREA.  AMAX 
INTEGER  ARRAY  (  lUU,  lw  >  ,  SUM  <  256  ) ,  AMAX  (5,  20  ) 

COMMGN/PLT/ 1 START X ,  ISTARTY,  MIDX,  MIDY,  NXSIZE, NYSIZE 
INTEGER  B(2, 1500) , UTT, UT1 


MAKE  SURE  THERE  IS  ENOUGH  PERIMETER  TO  REALLY  CONSIDER 
THIS  AN  AREA.  SMALLER  RUNS  INTO  PROBLEMS  IN  THE 
METHOD  USED  FOR  CALCULATING  SLOPES 

IF ( NB.  LE.  20) GO  TO  9UO 


PI =3  14159 
PI 3*3  *PI/8 
PI 5=5  *PI/8. 

IP.JP  IS  FIRST  POINT  ON  EDGE.  NAREA  IS  THE  "BASE  VALUE"USED  TO  MARK, 
I  E  NUMBER,  EDGES.  CALCULATE  IN  NC  THE  VALUE  THAT  HAD  BEEN  USED  TO 
MARK  THIS  EDGE. 

IP=B<  1,1) 

JP*B(2, 1 ) 
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V  - 


B 


NC=APPA  I  'IF  -  JF  i  '  tinREn 

FROM  THAT-  CALCULATE  THE  OFFER  AND  LOWER  LIMITS  OF  THIS  SEGMENT. 


LTL=NC*NAREA 

IF  LTL  LT  1  o24  >  STOP'  GONG 

UTT=LTL+NAREA 

LTl=LTL  +  2Gc- 

UT1=LT 1+255 

LH1=UT 1+1 

LHU=LH 1+255 

'M INR- 145 

XMIN=20000. 

START  WITH  EACH  POINT  ON  EDGE  OF  AREA  FOR  SEARCH  PATTERN 

DO  75  I  B=  1  >  NB 

SAVE  SEARCH  START  POINT 


IF'=B<  1. IE) 

■JP=B ( 2-  IE) 

IF  POINT  IS  ON  EDGE.  OF  SCAN.  IGNORE  IT. 

IF < IP  EQ.  1 ) GO  TO  75 
IF  (IF'  EQ.  100)  GO  TO  75 
IF  ( JP  EQ.  1  >  GO  TO  75 
I F  ( -JF'.  EQ.  1  Oo )  GO  TO  75 

FIND  SLOPE  OF  PERIMETER  IN  AREA  OF  THIS  POINT  BY  LOOKING 
AT  THE  LINE  FROM  S  POINTS  BACK  TO  S  POINTS  AHEAD. 

KPS®  I E 
KM8= I E 
DO  31  M=1,S 

AS  WE  MOVE  BACKWARD  ALONG  POINTS  IN  B,  IF  WE  REACH  B(l,l>  THEN 
WILL  NEED  TO  LOOP  AROUND  TO  END  OF  B(l.NB)  TO  PICK  UP  NEXT  POINT. 
THIS  IS  BECAUSE  B'l-1)  AND  B(1,NB)  ARE  ADJACENT  POINTS  WITHIN  THE 
GEOMETRY  OF  THE  SCANNED  ARRAY. 

) M8=KMS- 1 

I F ( KM8  LT  1 )KM8=NB+KM8-1 


IF  EDGE  OF  SEGMENT  HIT:-  EDGE  OF  SCAN.  THEN  QUIT  MOVING  BACKWARD 
BECAUSE  SCAN  EDGE  POINTS  mPE  NOT  RElAVENT  TO  SHAPE  OF  THE  SEGMENT 

IF ( Bf 1  \  M3 i  EO  1  GO  TO  22 
I F ( B ‘ 2 • >  MS i  EO  1  'Go  TO  22 


B 


A 
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IF*  B*  1  ■  K MS  > .  EG'  1*  .»«.•  >GO  TO  32 

IF<B*2.  KMS>  EQ  1  /GO  TO  32 

CONTINUE 

CONTINUE 

DO  33  M= 1 • S 

MOVE  FORWARD  ALONG  ARR.Ay  B  WITH  SAME  CONSTRAINTS  OBSERVED  WHEN  ' 

MOV I NG  BACK WAR DS 

K  PS=K  PS+  1  ! 

I F  ( K  PS.  GT.  NB  >  KPS— KPS— NB  ] 

I F  ( B  (  1 ,  KPS ) .  EG!  1  > GO  TO  34  J 

IF(B(2,  KPS).  EG!.  l)GO  TO  34  Ij 

I F  <  B  (  1 ,  KPS  )  EQ.  1  *:■>.' )  GO  TO  34 
IF <  B<  2.  KPS >  EQ.  l'X'.'GO  TO  34 
CONTINUE 
CONTINUE 

THE  DIFFERENCE  IN  COORDINATES  OF  THE  TWO  POINTS  FOUND  GIVE  US 
A  DIRECTED  SLOPE  OF  THE  EDGE  IN  THE  AREA  OF  THE  POINT  I, J. 

I DI F=B ( 1. KPS ) — B 1 1- KMSi 
•JD I F=B  (  2,  KPS  )  -B  (.  2  •  KMS  > 

XDI=FLOAT( ID1F) 

XDJ=FLOAT<  JDIF,' 


IF  FIRST  TIME  THROUGH  LOOP,  GO  ON  TO  FIND  A  TEST  MINIMUM 
IF  (IB.  EC!.  1  )  GO  TO  49 


EASED  ON  THE  SLOPE  JUST  CALCULATED,  PROJECT  A  POINT  THE 
MINIMUM  THICKNESS  ALREADY  FOUND  TOWARD  THE  INSIDE  OF  THE 
AREA.  IF  THAT  POINT  IS  INSIDE  THE  AREA,  THEN  THERE  IS  NOT 
A  NEW  MINIMUM  SO  WE  CAN  SKIP  REST  OF  TEST  FOR  THAT  POINT. 


IF ( ABS ( XDI )  LT.  ABSl XDJ  >  >G0  TO  43 

FIND  ABSOLUTE  RATIO  OF  SMALLER  DELTA  TO  LARGER 

BY  REVERSING  THE  I  AND  J  VALUES,  AND  CHANGING  THE  SIGN  OF  THE 
NEW  I  VALUE-  WE  ARE  PROJECTING  AT  9u  DEGREES  TO  THE  SLOPE. 

R J=  *  D I 'ABS ( x  D I ) 

P I  =—  x  D  J '  ABS  (  x  E*  I  ) 

GO  TO  44 
R J= x  D I / ABS \ X  D J ) 

RI=— XDJ/ABS( XDJ / 
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I  i  =  iP+IFI  x .  ;-MINF:-<PI  i 
Jl=JF+IFIx  •;  XMINR*RJ 
IF ■ 1 1  LT  l • 11=1 
IF'  Jl.  i_T  1  >  Jl  =  l 

IF  ' ii  gt  i >  11  =  1 

IF  v  .Jl  GT  luu  i  J 1  =  1 
PLOT  SPOT 

[  xF'=  I  ST  ART  X  + 1  l*NxSIZE 
I V  F'=  I  START Y-J 1  #N'V  SIZE 
CALL  SF'DOT  (  I  XF,  I  y  F) 

I V=ARRA Y  <  1 1 ■  Jl  ■ 

IS  POINT  INSIDE  AREA.  IF  SO  SKIP  TO  75. 
I F  v  IV  GE  LTL  AND.  IV.  LT.  UTT>GO  TO  75 
CONTINUE 


FIND  THE  ANGLE  ASSOCIATED  WITH  THE  LINE. 

I F  <  ADS  (  XD I  ) .  GT.  1 .  E-4  >  GO  TO  8 
SL0PE=3000  *FLOAT  v  JD I F  i 
PH  I  =P  1/2. 

IFfJDIF  LT  1 >PHI=PI+PHI 

GO  TO  9 

CONTINUE 

SLOF'E= xDJ/x  DI 

F'H I  =AT AN2 XDJ.  >DI ) 

CONTINUE 

SET  APPROPRIATE  SEARCH  ANGLE 


F:F'HI  =P  1 3+PHI 
QF'H  I  =P  1 5+F'H  I 


FIND  SLOPES  OF  SEARCH  LINES 

F:T  =  TAN(  PPHI  > 

OT  =  T  AN  '.  OF  HI  ) 


BASED  ON  SLOPES  AND  WHETHEF  THE  DELTA  I  OR  DELTA  J  SHOULD  BE 
LARGER ■  GENERATE  <_ OF'REC T  DELTae 

IF  aDS <  FT  ■'  L’E  1  1  oJ  i  u  1 

FDI=-1  O 
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SUBROUTINE  SFTi'J'v  I  ‘POS  I  i  Fij:-.  ■ 

'HI  :•  SUSP1  HJTINE  PETS  THE  POSITION  OF  THE  NE  *  T  CHARACTER  TO  EE 
F'F'  I  r.‘  TED  ON  THE  TERM  I  NHL  TO  I  « POS.  I  YPOS. 

the  FUNCTION  IS  SIMILAR  TO  SF'LiN,  BUT  ONLY  THE  FIRST  v  DARK ! 

VECTOR  i £  PLOTTED 

THE  TEKTRONIX  4..*06  MUST  HAVE  BEEN  ASSIGNED  TO  CHANNEL  1 
CALL  FOF'EN  •  1,  "STTOl"  >  "  .  BEFORE  CALLING  THE  SUBROUTINE. 

I  YH= 239+  <  I  YPOS/ 32  > 

COMPUTE  Y  START  HIGH  BYTE 

I  y  '_=  352+MOD  (  I  y  POS  ■  32 

COMPUTE  Y  START  LOW  BYTE 

I  *  H=28S+ \ I  X  POS >  32 ) 

COMPUTE  X  START  HIGH  BYTE 

I  x L  =  320+M0D ( I xPOS> 32 > 

COMPUTE  x  START  LOW  BYTE 

FORMAT  < S2.  Z > 

OUTPUT  1  CHARACTER. INHIBIT  CARRIAGE  RETURN 

FORMAT  >  S2. S2. S2. S2,  l ) 

OUTPUT  4  CHARACTERS.  INHIBIT  CARRIAGE  RETURN 

IGF:aP=2S5 

GRAPHICS  MODE  CHARACTER 

WRITE  ( 1  i  "  " 

OUTPUT  A  SPACE  FOLLOWED  BY  A  CARRIAGE  RETURN 

WRITE  (1,5)  I GRAF 

SHIFT  TO  GRAPHICS 
WRITE  (1.6)  I  YH.  IYL-  IXH,  I  XL 

MOVE  CURSOR  TO  START  POSITION 

I3F-23S 

NUMERICAL  REPRESENTATION  OF  A  SPACE 

I US =2 87 

NUMERICAL  REPRESENTATION  OF  ALPHA  MODE  CHARACTER 
WRITE  >1.6)  I SP ■  I SP ,  I  SR  I  US 

OUTPUT  SPACES  TO  ALLOW  TIME  FOR  VECTOR  AND  THEN 

GO  ALFHA  MODE 

RETURN 

END 


SUBROUTINE  SFLIN  .  I  ■  FOS.  IiPOS-  IxEND-  I  y  END  > 
p-,1;,  SUBROUTINE  CONVERTS  >  POSITION  AND  Y  POSITION  TO 

the  ppopef  cons tan r.-  ro  produce  a  plot  on  a  Tektronix  400&. 

FHE  PLOT  WILL  EC  FROM  I <POS-  I rPOS  TO  I  *  END  IyEND 

THE  T6>  TF'UNlX  4>. MOOT  HAVE  SEEN  ASSIGNED  TO  CHANNEL  1 

.  ••CALL  F OPEN  1  1,  '  T  T To  1  '  >  “  ■  BEFORE  CALLING  THE  SUBROUTINE. 

I i EH=CSS+ > I i END ■ SC i  ; COMPUTE  V  END  HIGH  BYTE 

I T h=2SS+ '  I YPOS/ S2 >  ; COMPUTE  V  START  HIGH  BYTE 

I  tEL=?'?2+M0D'  IyEND.  SC  ■  ;  COMPUTE  Y  END  LOW  BYTE 

I  YL  =  S52+M0D  I  (FOS.  S2  1  ;  COMPUTE  Y  START  LOW  BYTE 

I <EH=2SS+> I xEND' S2;  .COMPUTE  x  END  HIGH  ByTE 

I  x H-CSS+  I  XPOSx  32  )  ;  COMPUTE  X  START  HIGH  BYTE 

I yEL=32o+M0lM  I XEND*  32 >  > COMPUTE  X  END  LOW  BYTE 

I  »L= 320+MOD <  IXF'OS.  32  >  ,  COMPUTE  X  START  LOW  BYTE 

FORMAT  < S2  C'  ;  OUTPUT  1  CHARACTER.  INHIBIT  CARRIAGE  RETURN 

FORMAT  ( SC.  SC.  S2- S2.  L >  ■ OUTPUT  4  CHARACTERS,  INHIBIT  CARRIAGE  RETURN 

IGRAP=CS5  .GRAPHICS  MODE  CHARACTER 

WRITE  . 1 <  "  "  , OUTPUT  A  SPACE  FOLLOWED  BY  A  CARRIAGE  RETURN 

WRITE  ( 1 . 5  >  I GRAF  SHIFT  TO  GRAPHICS 

WRITE  '’l.ii  IYH.  IYL.  IXH.  IXL  ;  MOVE  CURSOR  TO  START  POSITION 

WRITE  ( 1 . S )  IYEH. IYEL- IxEH. IxEL  DRAW  LINE  TO  END  POSITION 
ISF'=2SS  ;  NUMERICAL  REPRESENTATION  OF  A  SPACE 

I US=2S7  .NUMERICAL  REPRESENTATION  OF  ALPHA  MODE  CHARACTER 

WRITE  <l.o '  ISP.  ISP.  ISP.  I US  ; OUTPUT  SPACES  TO  ALLOW  TIME  FOR  VECTOR 

AND  then  GO  ALPHA  MODE 

RETURN 

END 


NTOp=hRF:h  v  .  I .  j) 
GO  TO  22' 

NBOT  =hPF'h'C  <!-._•  i 
00  TO  _ 21.1 


AT  LABEL’;.  230  AND  25>.',  WE  HAVE  BEEN  MOVING  BETWEEN  LEVELS.  WHEN  THE 
VALUE  OF  THE  SLOPE  FALLS  BELOW  THE  LIMIT  NSTP,  THINGS  ARE  “LEVELING" 
OUT  IF  THERE  WAS  A  PREVIOUS  LEVEL  (NOT  STARTING  A  SIDE  OF  SCAN), 
THEN  GET  AVERAGE  BRIGHTNESS  AND  SAVE. 

IF (NS  LT  -NSTP ) GO  TO  224 

IF  <  NTOF'  LT  i.> )  GO  TO  231 

X  A=  (  FLOAT  >  NTOP+ARRAY  (  I  ■  J  )  )  i  ,/2. 

AVG=AVG+  X A 
CNT =CNT  + 1 . 

NTOF'=- 1 
NBOT =- 1 
GO  TO  220 

IF  (NS.  GT  NSTF  )  GO  TO  224 

IF  (NBOT  LT.  0 )  GO  TO  231 

XA=  ( FLOAT  <  ARRAY  (  I  •  U )  +NBOT  >  >  /  2. 

AVG=AVG+  X A 
CNT=CNT+1. 

GO  TO  231 
SLOF'E=NS 

IF  (  I  AES  <NS).  LT.  NSTP)  SLOPE-O 

CONTINUE 

CONTINUE 


THATS  ALL  THE  POINTS.  FIND  AVERAGE  FOR  UPPER  LIMITS. 

IF  (CNT.  LE.  1.  ) GO  TO  933 
LT= I F I  X ( AVG/CNT  > 

WRITE(LF',  100) MIN,  LT 

FORMAT ( "OCOMB I NED  CONTOUR  THRESHOLDS  BETWEEN  ",  14,"  AND  ",I4) 

RETURN 

LT =U 

M I N— 0 

RETURN 

END 


20 

24 

25 

C 

c 

c 

I* 

c 

c 

c 

c 

c 

c 


c 

c 

r 

C 

c 

r 

C 

c 

c 

c 

c 

c 

c 

c 

c 

24*.. 


IF  •  Nr.  GT  NSTP  'GO  TO  24 

IF-NBGT  LT.  0  GO  TO  31 

X  A=  *  FLOAT  '  ARRAY  *  I  •  U  )  +N3(U  T  .  .  /  2 

AVG=AVG+Xh 

CNT=CNT+1 

GO  TO  3 1 

SLOPE=NS 

IF  <  I  AES .  NS  >  LT  NSTP  .  SLOPE*-.- 


CONTINUE 

CONTINUE 


DO  225  1  =  10.  lOu,  l'.« 

AS  WE  EXAMINE  EACH  FOW.  NS  REPRESENTS  THE  NEW  SLOPE  FOR  THE  NEXT 
POINT  TO  EXAMINE  WHILE  SLOPE  IS  THE  GENERAL  SLOPE  IN  THE  REGION 
WE  ARE  SEARCHING  A  SLOPE  OF  ZERO  INDICATES  WE  APE  ON  THE  BOARD, 

A  NEGATIVE  SLOPE  INDICATES  MOVEMENT  TOWARDS  A  BRIGHTER  LEVEL,  AND 
A  POSITIVE  SLOPE  TOWARD  A  DARKER  LEVEL. 


NS = AFP  Ax  1  I  .  1  i  —  APF,-t  x  I  •  S  Ti  f  1  i 
SLOPE=NS 

IF',  I  AES*  NS  >  lT  NS  TP  SLOPE*-.- 

NT0F=-1 

NE0T=- 1 

IB=*.» 

I T =0 

DO  22 4  .J=1LPEND 

FIRST  CHEO  NEW  POINT  FOP  MAXIMUM  BRIGHTNESS 

IF  < ARRAY  *  I, UK  LT  MIN )MIN=aRRAY *  I, J) 

NS=ARRAY ( I , J ) -ARRAY ( I . J+STP ' 

JUMP  DEPENDENT  ON  WHICH  REGION  WE  h«VE  BEEN  IN. 

IF  i  SLOPE)  230.  24",  25-.' 


HAVE  BEEN  SEARCHING  ALONG  LEVEL  v SLOPE  = 

IF  THE  ABSOLUTE  VALUE  OF  THE  SLOPE  EXCEEDS  THE  LIMIT  NSTP, 
THEN  CONSIDER  STARTING  MOVE  TO  NEXT  LEVEL. 


SAVE  THE  PRESENT  VALUE  If)  NTOP  OR  NBOT  DEPENDENT  ON  WHICH 
DIRECTION  WE  ARE  STARTING  TO  MOVE 

I F  *  I ABS  *  NS  )  L T  NSTP  >  GO  TO  224 
IF  (NS.  GT  O.'GO  TO  245 
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NS=ARRa  y  <  1  •  J  )  -ARRAY  (  STP+  1  ■  .J  ? 

SLOPE=NS 

IF  (  IADS  >  MS  i  L  T.  NS  TP  >  SLOPE**.' 

ntop*— i 

NEC'T=-1 

i  b=o 

I  T=0 

DO  24  1=1. LPEND 

FIRST  CHECK  NEW  POINT  FOR  MAXIMUM  BRIGHTNESS 

IFiARRAy  i  I,  J !  GT.  MAX  )  MAX=ARRAY  <  I ,  J) 

IF (ARRAY (I. J).  LT  MIN ) MI N=ARRAY ( I , J ) 

NS* ARRAY v I . J > -ARRAY ( I +S  TP . J ) 

JUMP  DEPENDENT  ON  WHICH  REGION  WE  HAVE  BEEN  IN. 

IF ( SLOPE) 30. 40. 50 

HAVE  BEEN  SEARCHING  ALONG  LEVEL  (SLOPE  =  0). 

IF  THE  ABSOLUTE  VALUE  OF  THE  SLOPE  EXCEEDS  THE  LIMIT  NSTP, 
THEN  CONSIDER  STARTING  MOVE  TO  NEXT  LEVEL. 


SAVE  THE  PRESENT  VALUE  IN  NTOP  OR  NBOT  DEPENDENT  ON  WHICH 
DIRECTION  WE  ARE  STARTING  TO  MOVE. 

IF<  I  AES  (NS).  LT.  NSTP )  GO  TO  24 
IF  (NS.  GT.  0 >  GO  TO  45 
NTOP=ARRAY ( I , J ) 

GO  TO  20 
NBOT*ARRAY ( I , J  > 

GO  TO  20 


AT  LABELS  30  AND  50,  WE  HAVE  BEEN  MOVING  BETWEEN  LEVELS.  WHEN  THE 
VALUE  OF  THE  SLOPE  FALLS  BELOW  THE  LIMIT  NSTP,  THINGS  ARE  "LEVELING" 
OUT.  IF  THERE  WAS  A  PREVIOUS  LEVEL  (NOT  STARTING  A  SIDE  OF  SCAN), 
THEN  GET  AVERAGE  BRIGHTNESS  AND  SAVE. 

IF <NS.  LT.  -NSTP) GO  TO  24 

IF  (NTOP.  LT.  0 )  GO  TO  31 

X A= ( FLOAT ( NTOP+ ARRAY ( I , J ) ) > /2. 

AVG=AVG+XA 

CNT*CNT+1. 

NTOP*— 1 
NBOT =- 1 
GO  TO  20 


o  o  o  o  o  o  o  o  o  o  o  o  *  >  o  o  o  o  *.  >  o  o  o  r>  o  o  o  o  o  o  o  r»  o  c 


OVEFLA i  OVER* 

SUBROU  r  I NE  CONVAL  .  »R R«  r  •  I..P  M I N,  L  T  > 

C  PREPARED  FOP  -'O.I  3  f  S  TEM  r.  •  I  NC 

c  by  sill  pope,  telcom  d«ta  corporation. 

THIS  SUBROUTINE  13  BEING  MAINTAINED  IN  FILE  SVCONVAL  AND  IS 
MOVED  TO  CONVAL  FOR  COMPILATION  WHEN  DESIRED  THERE  IS  AN 
ALTERNATE  VERSION  OF  CONVAL  KEPT  IN  FILE  CONSAVE. 


THIS  SUBROUTINE  TRIES  A  DIFFERENT  METHOD  FOR  THRESHOLDS. 

ASSUMMING  THAT  THE  DATA  ARRAY  RECEIVED  FROM  A  SCAN  BASICALLY 
FORMS  A  CONTOUR  OF  THE  BRIGHTNESSES,  CONHIST  SEARCHES  ONE  ROW 
AT  A  TIME  LOOKING  FOR  THE  CHANGES  IN  CONTOUR.  THE  CHANGING  CONTOUR 
CAN  THEN  BEE  SEEN  AS  LEVELS  REPRESENTING  THE  BOARD.  THE  PUNS,  AND 
THE  "BOTTOM"  OF  HOLES.  BY  FINDING  THE  AVERAGE  MIDPOINT  BETWEEN  LEVELS 
EACH  TIME  THE  CONTOUR  MOVES,  THESE  MIDPOINTS  CAN  THEN  BE 
AVERAGED  TO  FIND  THE  "UPPER”,  OR  DARKER,  THRESHOLD.  AT  THE 
SAME  TIME,  THE  BRIGHTEST  SPOT  ON  THE  BOARD  CAN  BE  FOUND  TO  USE 
AS  THE  LOWER  THRESHOLD. 


INTEGER  ARRAY  (  lOU,  1  <.'<  «)  •  STR.  SLOPE 


THE  METHOD  FOR  DETECTING  MOVEMENT  BETWEEN  LEVELS  IS  EXAMINATION  OF  THE 
SLOPE  OF  THE  CONTOUR.  SINCE  THE  "X"  COORDINATE  WILL  BE  CONSTANT  FOR 
EACH  COMPARISON,  ONLY  THE  "Y"  COORDINATE  *  DIFFERENCE  IN  BRIGHTNESS ) 
WILL  BE  EXAMINED. 


INITIALIZE  VALUES 

NSTP= 10 
STP*4 

LPEND= 1 OO-STP 
MAX*0 
MIN-255 
CNT =0. 

AVG=U. 

DO  25  0*10, 100, lo 

AS  WE  EXAMINE  EACH  ROW,  NS  REPRESENTS  THE  NEW  SLOPE  FOR  THE  NEXT 
POINT  TO  EXAMINE  WHILE  SLOPE  IS  THE  GENERAL  SLOPE  IN  THE  REGION 
WE  ARE  SEARCHING.  A  SLOPE  OF  ZERO  INDICATES  WE  ARE  ON  THE  BOARD, 

A  NEGATIVE  SLOPE  INDICATES  MOVEMENT  TOWARDS  A  BRIGHTER  LEVEL,  AND 
A  POSITIVE  SLOPE  TOWARD  A  DARKER  LEVEL. 
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IF  POINT  IS  IN  SOME  OTHER  mP.Eh  GO  SET  FLAG. 


IF\  IV  LT.  LTH'GO  TO  lV 
I F  •  I V.  GT  IJTH  /  GO  T 0  1  V 
RETURN 


MTOUCH= 1 

FOUND  CRITICAL  DISTANCE-  DRAW  VECTOR  TO  SHOW  AND  MARK  IT  WITH  "C" 

IX=ISTARTX+I*NXSIZE 
I  y'=  I STARTY-  J*NYS  I ZE 
I  X S= I START X+ 1 0*NX  S I ZE 
I YS= I START Y -JO*NYS I Z E 
CALL  SPL IN  < IXS.  I YS>  I*-  IV ) 

I XE=I XS- > I XS-IX  >/2 
IYE=IYS-( IYS-IY)/2 
CALL  SPMOVCIXE.  I  YE  > 

WRITER  1-  1 00 ) 

FORMAT  (  ’*  C"  ) 

RETURN 

END 


O  O  O  O  O  O  O  O  u  u  O  U  M  ;  O  O  i  \  o  o  o  o 


SUBROUTINE  CKPT  i  10,  JO,  1 1  •  •  J 1 ,  MTOUCH,  LTH,  IJTH  > 

PREPARED  FOP  SCI  SYSTEMS.  INC 
B'i  BILL  POPE .  TELCOM  Da TA  CuRF’URATI UN 


SUBROUTINE  CkPT  PREPARES  A  "PPROJECT ION”  FROM  POINT  10, JO  ALONG 
A  DELTA  OF  II, Jl.  IF  THE  NEW  POINT  IS  WITHIN  SOME  AREA  OTHER  THAN 
THE  CURRENT  ONE  < CURRENT  AREA  HAS  THRESHOLDS  OF  LTH, UTH),  THEN 
THE  FLAG  MTOUCH  IS  SET  INDICATING  CRITICAL  DISTANCE. 


C  IF  PROJECTION  LOGIC  FOUND  IN  TRACK  IS  UNCLEAR,  TRY  PLOTTING  POINTS 

C  FROM  THIS  SUBROUTINE  TO  SEE  DIRECTIONS  PROJECTED. 


COMMON  ARRAY,  SUM.  LP,  IMN.  IMX,  .JMN .  JMX.  NCNT ,  NAREA.  AM  AX 
INTEGER  ARRAY  v  i«X»,  1  •.»•.*> .  SUM <  25ft, ) ,  AMAX  (5,  2*.») 

INTEGER  UTH 


COMMON  /PL T/ 1 START X,  ISTARTY,  MIDX,  MIDY, NxSIZE. NtSIZE 


FIND  NEW  PROJECTED  POINT. 

1=10+11 

J-J0+J1 


IF  PROJECTION  OFF  SCAN,  BACK  UP  TO  EDGE. 

IF X  I.  LT  1  )  1  =  1 
IF  ( J.  LT  1  )  J=1 
IF (I  GT.  100)I=1U0 
I F  ( J.  GT  1 00 )  J=  1 OU 


IF  DESIRED,  PLOT  PROJECTED  POINT. 

IX=ISTARTX+I*NXSIZE 
I Y=ISTARTY-J*NYSI ZE 
CALL  SPDOT (IX, It) 


IV=APRAr ( I. J > 

IF  POINT  ON  BOARD  (NOT  MARKED  FOR  ANY  AREA > ,  RETURN. 


IF<  IV.  LT.  NAREA; RETURN 


B-8 


RETURN 

WR I  TE  *  LF'»  903  1 
FORMAT \  "  AREA 
RETURN 
END 


THICKNE 


o  o  o  o  o  o  ^  o  n  n 


NDPT=MF'NT 
MPNT=NSTRT 
00  TC  70 


i 


OUT'S  I  Cl  OF  AREA,  NOW  WE  RE  IN  TROUBLE 

SEE  IF  AN  ADJACENT  POINT  IS  ON  BOUNDARY.  IF  SO  , 

•70  IF (  IN.  EQ.  100)00  TO  91 

IU=ARRAY  (  IN+1,  .JN) 

IF ( IU.  LT.  LT1.  OR.  IU.  OT  UTl/GO  TO  91 

IN=IN+1 

00  TO  95 

91  IF' JN.  EQ.  1  •  00  TO  92 
IU=ARRAY( IN. JN— 1 ) 

IF  i  IU.  LT.  LT1.  OR.  IU.  GT.  UT1-G0  TO  92 
JN= JN- 1 
00  TO  95 

92  IF ( JN.  EQ.  100)00  TO  93 
IU=ARPAY( IN, JN+1 ) 

IF (  IU.  LT  LT1.  OR.  IU.  OT.  UT1)00  TO  93 
JN= JN+ 1 
GO  TO  95 

93  IF <  IN.  EQ.  1  >00  TO  94 
IU=ARRAYv IN-1, JN) 

IF <  IU  LT.  LT1.  OR.  IU.  GT.  UT1  >00  TO  94 
IN=IN— 1 

95  NBFLA0=1 

GO  TO  64 

94  CONTINUE 

IF  NOT  ON  EDGE  OF  SCAN  JUST  KEEP  TRUCK IN 

IF (  IN.  EQ.  1  ) GO  TO  75 
IF IJN.  EQ.  1  )G0  TO  75 
IF (  IN  EQ  1 00 ) GO  TO  75 
IF ( JN  EG  100)00  TO  75 
IF ( NBFLAG.  EQ.  0)00  TO  60 
5  CONTINUE 

DID  WE  FIND  A  MINI NUN 

IF<  XMIN.  EQ.  20O00.  )G0  TO  900 

YES,  PLOT  IT 


USE  IT 


X  =  XMIN 
XMIN=SQRT ( X ) 

CALL  D0TLINU3I,  I3J,  ICI,  ICJ) 


I C  .J=UN 
NBFL AG= i 
CONT  I  NUE 

IF  NOT  TRACT.  1  NO  Ai_ONO  EDGE.  GO  MOVE  TO  NEW  POINT 
IF  •;  NBFLAG  EQ.  '.‘/GO  To  c1-1 

HAVE  WE  DEEM  UN  EDGE  BEFORE 
IF  (MPNT.  ME.  -.••GO  TO  7v 

SEARCH  THROUGH  BOUNDARY  ARRAt  TO  FINE*  WHERE  WE  ARE  LOCATED 

DO  66  M 1  =  1  ■ NB 

IF  (  IN  NE  B-  1,  Ml  •  •  GO  TO  c-6 

IF  (.JN  NE  B  ‘  2 •  Ml  )  )  GO  TO 

MF'NT=M1 

GO  TO  80 

CONT I NUE 

GO  TO  75 


MOVE  1  POINT  ALONG  EDGE 
MPNT =MF'NT  + 1 

IF  (  MF'NT  EG.  NO  FT  )  GO  TO  75 
IF  (MPNT  GT.  NE  ;  GO  TO  7  5 
IN=E( 1  MPNT ) 

JN=E ( 2 • MPNT) 

00  7  0  81 


SEARCH  OUTWARD  FROM  POINT  HAS  FOUND  AN  EDGE.  FIND  IF  IT  IS  BEGINNING  OR 
END  POINT  OF  EDGE 

IF ( NSTRT.  NE.  0>G0  TO  S5 
IF (MPNT  LT.  IB) GO  TO  75 

SAVE  START  POINT  OF  BOUNDARY  SEARCH  AND  SETUP  TO  GO  FIND  END  POINT. 

NSTRT =MF'NT 
RDI=QDI 
F:DJ=QDJ 
GO  TO  59 


WE  HmVE  THE  ENDPOINT  ON  BOUNDaPi  CAVE  IT  THEN  GU  SEARCH 
BOUNDARY  FROM  NSTRT  TO  NDF  F. 


O  O  O  O  O  O  O  o  o 


IF  NOT.  LOOP  BACK  FOR  NEXT  POINT 

IF  •  IV  LT  L  T 1  >00  TO  60 
IF'  IV  GT  UT1 >00  TO  *0 


C  IF  SEARCH  POINT  IS  WITHIN  9  POINTS  OF  ORIGINAL  ALONG  BOUNDARY  OF 

C  AREA,  DON'T  CONSIDER  IT. 

C 

MS* IB-9 

IF  <  MS.  LT.  1  )MS=MS+NB 
DO  62  Ml  =  l.  IS 
MS=MS+ 1 

IF  v  MS.  GT  NB  > MS=MS-NB 
I F  (  I N.  NE.  B  1 ,  MS  >  >  GO  TO  62 
I  FUN.  NE.  B  v  2,  MS )  )  GO  TO  62 

C 

C  TOO  CLOSE 

IF  < NSTRT.  NE.  0)00  TO  64 
C 

C  IF  JUST  RUNNING  ALONG  EDGE,  DROP  POINT. 

C 

MOVE=MOVE+ 1 
IF  (MOVE.  GT.  5) GO  TO  75 
C 

GO  TO  64 
62  CONTINUE 

NBFLAG=1 


IFvMPNT.  EC*.  0 ) GO  TO  67 

MAKE  SURE  NOT  CLOSING  ON  EACH  OTHER. 

IF < ( MPNT-IB ) .  LT.  15; GO  TO  70 

IF<  IB.  LT.  8.  AND.  ( NB-MPNT ) .  LT.  8) GO  TO  75 


FOUND  A  POINT,  CHECK  ITS  DISTANCE 

67  XDI»1.  25*FL0AT( IN-IP, 

XDJ*FLOAT ( JN-JP > 

D I  ST = X  D I  **2+xDJ**2 

IF  (  DIST.  GT.  XMIN ) GO  TO  64 

C 

C  FOUND  NEW  MINIMUM 

XMIN=DIST 
XM I NR=SQRT (DIST ) 

IBI=IP 
I B J= JP 
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RMJ-PM J+RDJ 

I F  >  ABS  ( RMJ  .  LT.  1 .  vi  GO  TO  58 
PM  _I=RMJ— FLOAT  JD  » 

00  TO*  57 
C ON  T I NUE 


SET  FLAG  FOR  NOT  SEARCHING  ALONG  EDGE  CURRENTLY 

NEFLAG=0 
MR NT =0 

MAKE  A  MOVE 
IN=IN+ ID 
.JN= JN+ JD 
CONT I NUE 

PLOT  THE  POINT  BEING  EXAMINED 


I  X  P= I START  X  + 1 N*N X S 1 2E 
I YP» I START Y-JN*NYS I L E 
CALL  SF'DOT (  IXP.  IYP) 

IF  STARTING  POINT.  STOP  LOOP 

IF  (  IN.  EG.  IP.  AND.  JN.  EG.  JP>60  TO  75 

FIND  VALUE  OF  POINT  BEING  LOOKED  AT 
I V=ARRAY (.  IN,  JN) 

MAKE  SURE  POINT  IS  IN  AREA 

IF(  IV.  LT.  LTL ) GO  TO  90 
IF':  IV.  GT.  UTT  ) GO  TO  90 

IF  NEW  POINT  ON  EDGE  OF  SCAN/  IGNORE  IT 

I F {  IN.  EG.  1  ) GO  TO  6 S'. 

IF  (  JN.  EG.  1  )  GO  TO  63 
IF  (IN.  EG.  100)  GO  TO  63 
I F  <  JN  EG.  100)  GO  TO  63 


IF  WE  HAVE  HIT  A  HOLE  IN  THIS  AREA,  IGNORE  THIS  THICKNESS 
IF  v  IV  GE  LH1  AND.  IV  LE.  LHU ) GO  TO  75 


FIND  OUT  IF  IT  IS  ON  BOUNDARY 
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RDU=-PT 

IF'  UD IF  OE  "  GO  Tu  11 


PD  1  =  1  V 
RDJ=-PDJ 
GO  TO  1 1 
PD 1=1  / RT 
RDJ=1  o 

IF  (ID  IF.  OE.  >0  ;  GO  TO  11 
RDJ=- 1  O 


RDI=-PDI 

I F  <  AE(S  •'  QT  >  OE.  1 .  u  >  OU  T 0 
OD I =- 1 .  O 
QOU= -OT 

IF  v  JDIF  OE.  V  >  GU  Tu  10' 

OD  1  =  1  O 
QDJ=-GDU 
GO  TO  10 
QD  1  =  1.  /QT 
OD  0=  1 .  O 

IF(  I  DIF.  OE.  0)00  TO  13 
000= - 1  0 
OD I =—00 I 
CONTINUE 
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INITIALIZE 

NSTRT=0 

NDPT=U 

MPNT =0 

M0VE=0 

CONTINUE 


NOW  WE  APE  READY'  TO  BEGIN  SEARCH.  SET  CURRENT  POINTS  (IN,  JN)  AND 
REMAINDER  FROM  LAST  MOVE 

I  N=  I P 
ON=.JP 
RMI=0.  0 
RMJ=0.  O 

I D= IF  I  X  < RD I +RM I )  '1 
PMI =RMI +RD I 

IF  <  AES'  RMI  •  LT  1  u  1  GO  TO  5c. 

RM I =RMI -FLOAT i ID ) 

GO  TO  55 

•JD=  IF  I  X  ( RDO+RMO  /  /  1 
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SUBROUTINE  DOTL IN>  MX, My, NX, NY) 

PREPARED  FOP  SCI  SYSTEMS.  INC. 

BY  BILL  POPE.  TELCOM  DAT A  CORPORATION 

THIS  SUBROUTINE  DRAWS  A  DOTTED  LINE  BETWEEN  2  POINTS 
GIVEN  BY  MX, MY  AND  NX, Ny  MX, MY  AND  NX, NY  ARE  GIVEN 
AS  ROW,  COLUMN  VALUES  FROM  THE  lUOXlOO  SCAN  MATRIX. 

THEY  ARE  CONVERTED  TO  SCREEN  COORDINATES  BASED  ON  THE 
INITIALIZATION  VALUES  GIVEN  IN  THE  FOLLOWING  COMMON  BLOCK. 

COMMON  /PLT/ISX. ISY, MIDX, MIDY, NxSIZE, NYSIZE 


I X«ISX+MX*NXSIZE 

I Y= I SY—M  Y *NYS I ZE 

JX=ISX+NX*NXSIZE 

•JY=ISY— NY*NYSI  ZE 

IXS=IX 

IYS=IY 

MAXE=ISX+100*NXSI ZE 
M I NY= I SY- 1 00*NYS I ZE 


DASH=2.  AFLOAT'. NYSIZE) 


CALCULATE  DISTANCE  BETWEEN  POINTS 

XDI  =  1.  25*FL0AT ( MX-NX ) 

XDU=FLOAT ( MY-NY ) 

X X » X D I **2+ XD J**2 
XD=SGRT (  XX ) 

IF  t  XD.  LE.  1.  75)  GO  TO  20 
SCRN  X  D=  X  D*FLOAT ( NYSI ZE ) 

IF  DISTANCE  BETWEEN  POINTS  TOO  SMALL  FOR  DOTTED  LINE, 
REDUCE  SIZE  OF  EACH  DASH. 

IF ( SCRNXD.  LT.  <3  *DASH )  ) DASH=SCRNXD/3. 

IF  (DASH.  LE.  0  )GO  TO  20 


CALCULATE  RECTANGULAR  DELTAX  AND  DELTAY  FOR  EACH  DASH 
IN  THE  CONNECTTING  LINE. 

PTS=SCRNXD/DASH 
RDX=-FLOAT< Ix-JX )/PTS 
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RDY=-FLOAT ( I t - JY ) j RTS 
ND  X  =  I F I  x  ( F.D  x ) 

IFfND X.  NE  «.• )  GO  TO  5 
NDX  =  1 

IF  (  PDX.  LT.  0.  )  NDX=- 1 
QDX=FLGAT<NDX > 

NDY=IFI X ( RDY ) 

IF  ( NDY.  NE  0)00  TO  6 
NDY=  1 

IF  (RDY  LT.  0.  ) NDY=—  1 
QDY=FLOAT ( NDY  ) 

RMX=0. 

RMY=0. 


TAKE  OFF  AND  PLOT  LINE 

I XE= I  X— IF I  X ( RDX  ) 

IYE=IY-IFIX (RDY) 

NDASH=7* I F I X ( DASH ) /4 
I XS=I XE+IFI X ( RDX+RMX ) 

I XE=I XS+IFI X (RDX ) 

I YS= I YE+ 1 F I X ( RDY+RMY ) 

I YE=I YS+IFI X ( RDY ) 

RMX=RMX+RDX 

IF ( ABS ( RMX ) .  GT.  ABSvQDx) )RMx=RMX-QDx 
RMY=RMY+RDY 

IF ( ABS ( RMY ) .  GT.  ABS(QDY) )RMY=RMY— QDY 
KX=IABS<  IXE-.JX  ) 

KY=IABS( IYE-JY) 

I F  (  K X .  LT  NDASH.  AND.  KY.  LT.  NDASH  )  GO  TO  20 
K  X  =  I  ABS  (  IXS-.JX  ) 

KY= I ABS ( IYS-JY) 

I F  (  K X .  LT.  NDASH.  AND.  KY.  LT.  NDASH  >  GO  TO  20 
IF  (  I  XE.  LT.  ISX.  OR.  I XE.  GT.  MAXE )  GO  TO  21 
IF  ( I  YE  LT.  MI  NY.  OR.  I  YE.  GT.  ISY)GO  TO  21 
CALL  SF'L  I N  ( I X  S »  I  YS»  IXE,  I  YE  ) 

GO  TO  1 0 


REACHED  END  POINT.  CONNECT  TO  IT 

CALL  SPL I N ( I  X  S .  I YS - J X , OY ) 
CONTINUE 


WRITE  SIZE  BESIDE  LINE 

I XE=NXSI ZE+ 1  X— ( IX-JXJ/2 
I YE= I Y— < IY-JY>/2 


overlay  overs 

SUB ROUT I NE  D I STANCE 

PREPARED  FOR  SCI  SySTEMS-  INC. 

BY  BILL  POPE'  TELCOM  DAT A  CORPORATION. 


AFTER  ALL  SEGMENTS  WITHIN  THE  SCANNED  ARRAY  ARE  MARKED,  DISTANCE 
IS  CALLED  TO  INSPECT  FOR  THE  POSSIBLE  DEFECT  OF  TWO  SEGMENTS  BEING 
TOO  CLOSE  TOGETHER.  FOR  SOME  PRESET  MINIMUM  CRITICAL  DISTANCE  (10.  ), 
A  VECTOR  WILL  BE  PROJECTED  OUTWARD  FROM  EACH  POINT  ON  THE  EDGE  OF 
THE  BOUNDARY  OF  EACH  SEGMENT  IN  THE  SCAN.  IF  THE  PROJECTED  POINT 
HITS  ANOTHER  SEGMENT,  IT  IS  ASSUMED  THAT  THE  SEGMENTS  ARE  TOO  CLOSE. 

ALL  REQUIRED  INPUT  IS  THROUGH  COMMON. 


COMMON  ARRAY, SUM. LP,  IMIN.  IMAX,  JMIN, UMAX, NCNT, NAREA, AMAX 
I NTEGER  ARRAY ( 1 OO,  1 oO  > . SUM ( 256 > . AMAX ( 5, 20 ) 


I  B=0 
IC=0 

NCNT  IS  NUMBER  OF  SEGMENTS  PREVIOUSLY  IDENTIFIED 

IF  (NCNT.  EQ.  1  )  GO  TO  990 
WR I TE ( LP, 200 ) NCNT 

200  FORMAT ( "0*  *  *  THERE  WERE  ",13,  "  AREAS  LOCATED.  " > 


WANT  TO  COMPARE  EACH  SEGMENT  TO  ALL  OTHER  SEGMENTS  SO  NEED 
NESTED  LOOP. 

ZMIN=10. 

ZMIN2=ZMIN**2 


DO  81  NB=1 , NCNT 
DO  80  NC- 1 , NCNT 
IF  SAME  AREA,  SKIP 


IF  ( NB.  EQ.  NO  GO  TO  80 

FOR  THE  NEW  AREAS  TO  BE  COMPARED,  GET  THE  MINIMUM  AND  MAXIMUM 
AREA  COORDINATES  SAVED  FROM  SUBROUTINE  EDGE. 

IMINB=AMAX ( 1, NB) 

IMAXB=AMAX (2, NB ) 
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JMINB=AMAX  <  3.  NB) 
JMAXE-AMAX ( 4. NB > 
IMINC=AMAx ( 1 , NC » 
IMAXC=AMAX <  2,  NC • 
JM I NC=AMA  X  ( 3 > NC / 
JMAXC=AMAX(4,  NO 


FIRST  .  ASSUME  A  EOx  GROUND  EACH  AREA.  IF  THE  BOXES  DON'T 
OVERLAP,  AND  THE  PERPENDICULARS  BETWEEN  THE  BOXES  IS  GREATER 
THAN  THE  MINIMUM  DISTANCE  FOUND  SO  FAR,  THEN  THE  BOX 
BEING  INSPECTED  MAY  BE  ELIMINATED. 


SET  THE  DIRECTION  WE  CAN  SEE  TO  U,  I.  E.  BOXES  OVERLAP 
DIRECTED 

IF (  IMAXB.  GE.  IMINOGO  TO  2 
X BUM= 1 .  25*FL0AT( IABS< IMAxB-IMINO ) 

IF  (  XDUM.  GT.  ZMIN  )  GO  TO  80 
D I RECTO 
GO  TO  8 

2  I F  <  JM  A  X  B.  GE.  JM  1  NC  >  GO  TO  4 

XDUM*FLOAT< I ABS ( JMA  X  B-JM I NO ) 

IF  (XDUM.  GT.  ZMIN)  GO  TO  SO 
DIRECT® 5 
GO  TO  8 

IF <  JMINB.  LE.  JMAXOGO  TO  6 
X  DOM®  FLO  AT  <  IABS(  JMINB-.JMAXC)  ) 

IF  (XDUM.  GT.  ZMIN)  GO  TO  80 

DIRECT=1 

GO  TO  8 

IF (  IMINB.  LE.  IMAXOGO  TO  8 

XDUM=1  25*FL0AT( IABS( IMINB-IMAXC) > 

IF  (XDUM.  GT.  ZMIN)  GO  TO  80 

DIRECT=7 

CONTINUE 


IF  BOTH  I  COORDINATES  ARE  ON  ONE  SIDE  OF  AREA  NB  AND 
ALSO  BOTH  J  COORDINATES  ARE  ON  ONE  SIDE  OF  AREA  NB,  THEN 
A  QUICK  EXAMINATION  OF  THE  DISTANCE  TO  THE  CORNERS  OF 
THE  AREA  MIGHT  ELIMINATE  IT. 

IF  <  IMAXB.  GT.  in  INC.  OR.  UMAXS.  GT.  UMINOGO  TO  12 
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,<  D I  =  1  25*FL0AT  (  I  MAX  B- 1 M I NC  > 

X  DJ=FLOAT  ( JMA.X  B-JM I  NC  ) 

xDUM=xDI**2+XD J**2 

IF'  XD'JM  GE  ZMIN2  >00  TO  8>.‘ 

DIRECT=4 
00  TO  18 

IF (  IMINB.  LT  IMAXC.  OR.  JMAxB.  OT.  JMINOGO  TO  14 
X  D I  =  1  25>FL0AT ( I M I MB- i MA X C  > 
xD.J=FLOAT  (  JMAX B-JM INC  > 

XDUM= xD I **2+ xDJ**2 
I F  ( X  BUM.  OE .  i  M 1  M2 )  GO  TO  SO 
DIRECT-6 
GO  TO  18 

IF  k  IMAXB.  OT.  IMINC.  OR.  JMINB.  LT.  JMAXOOO  TO  16 
X D I  =  1 .  25*FL0AT  ( I MA  XB- IMI NC ) 

X  DJ=FLOAT v  JM I NB- JMA X  C ) 

X  DUM= X  D I **2+X DJ*#2 

IF  (  XBUM.  GE.  ZMIN2  i  GO  TO  80 

DIPECT=2 

GO  TO  18 

IF(  IMINB.  LT.  IMAXC.  OR.  JMINB.  LT.  JMAxOGO  TO  18 
XDI=1.  25*FL0AT( IMINB-IMAXO 
XDJ=FLOAT ( JMINB— JMA XC ) 

XDUM=XDI**2+XDJ**2 

IF  (XBUM.  GE.  ZMIN2  >  GO  TO  80 

DIRECT=8 

CONT I  NlJE 


BOXES  TOO  CLOSE  FOR  AREAS  NB  AND  NC.  JUST  GO  CHECK  AREA  NB  THOROUGHLY. 


SET  THE  BOUNDARY  LIMITS  FOR  NB  SO  WE  CAN  IDENTIFY  AREA. 

L  T 1 =NB*NARE A+256 
UT1=LT 1+255 


FIND  THE  MAX  AND  MIN  COORDINATES  FOR  THIS  AREA  FROM  DATA  WE  SAVED 
IN  AEDGE. 

DO  10  I  =  I M I NB  >  IMAXB 

IF < ARRAY*  I,  JMINB).  LT.  LT1  >G0  TO  10 

IF  (ARRAY  (  I»  JMINB'.  LE.  UT1  >00  TO  11 


CONTINUE 

TV  PE  **  ERROR  IN  SUBROUTINE  DISTANCE" 
STOP 

CONTINUE 


WE  NOW  HAVE  A  PERIMETER  POINT  FROM  AREA  NB.  INITIALIZE 
SOME  VALUES  AND  CALL  TRACK  TO  CHECK  FOR  CRITICAL  DISTANCES. 
SUBROUTINE  TRACK  IS  A  SPINOFF  FROM  SUBROUTINE  TRACE.  IT  BASICALLY 
TRACES  THE  EDOES  OF  A  SEGMENT.  BUT  AS  EACH  POINT  IS  FOUND  MAKES 
THE  CRITICAL  DISTANCE  PROJECTION. 


10=1-1 
J0= JM INB 

CALL  TRACK < I. JM I NB •  10. JO. LT1 . UT1.  ZMIN, MTOUCH) 


DID  WE  FIND  A  CRITICAL  DISTANCE 
I F  ( MTOUCH.  EC!.  0)00  TO  81 
YES.  PRINT  MESSAGE 
WRITE  <LP.  1S1 /NB 

FORMAT ( "O***********  CRITICAL  DISTANCE  ADJACENT  TO  AREA  ".13) 

GO  TO  31 

CONTINUE 

CONTINUE 

RETURN 

WRITE <LP, 991 ) 

FORMAT ("  ONLY  ONE  AREA  SO  DISTANCES  UNAVAILABLE.  " ) 

RETURN 

END 


FILMED 


5-85 


DTIC 


